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Abstract 

The search for Higgs bosons in both the standai^d model and its extensions 
is well under way at the Tevatron. As the integrated luminosity collected in- 
creases into the multiple inverse femptobai^n range, these searches are becom- 
ing very interesting indeed. Meanwhile, the construction of the Large Hadron 
Collider (LHC) and its associated experiments at CERN are nearing comple- 
tion. In this TeV4LHC workshop, it was realized that any experience at the 
Tevatron with respect to backgrounds, experimental techniques and theoreti- 
cal calculations that can be verified at the Tevatron which have relevance for 
future measurements at the LHC were important. Studies and contributions 
to these efforts were made in three broad categories: theoretical calculations 
of Higgs production and decay mechanisms; theoretical calculations and dis- 
cussions pertaining to non-standard model Higgs bosons; and experimental 
reviews, analyses and developments at both the Tevatron and the upcoming 
LHC experiments. All of these contributions represent real progress towards 
the elucidation of the mechanism of electroweak symmetry breaking. 
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1 Introduction 

Contributed by: S. Willenbrock, A. Dominguez, I. lashvilli 

The Fermilab Tevatron, which has been coUiding protons and antiprotons for over twenty years, 
was not designed to search for the Higgs boson. However, the advent of high-efficiency b tagging, 
developed in the context of the search for the top quark, made it possible to consider searching for the 
Higgs boson, produced in association with a weak boson, via the decay h —>■ 66 [1]. It was realized 
that this would require very high luminosity, and that other discovery modes, such as /i — > W~^W~ — > 
£~^£~i'r', might also become viable with sufficient integrated luminosity [2]. The strategy for the Standard 
Model Higgs search was developed in the TeV2000 workshop [3], and was further refined, along with 
the case of the supersymmetric Higgs, in the SUSY/Higgs workshop [4]. 

The search for a Higgs boson, both standard and supersymmetric, is in full swing at the Tevatron, 
and is becoming increasingly interesting as the integrated luminosity mounts. Meanwhile, the con- 
struction of the CERN Large Hadron Collider (LHC) is nearing completion. At this workshop, dubbed 
TeV4LHC, the Higgs working group used the first meeting to decide what "TeV4LHC" means in the 
context of the Higgs boson. We decided that anything having to do with the Higgs at the Tevatron was 
relevant to the workshop, since this experience will surely be valuable at the LHC. Any experience at the 
Tevatron with backgrounds to Higgs searches is also relevant to the workshop. Finally, any experimental 
techniques being developed for the Higgs search at the Tevatron or the LHC should also be included in 
the workshop. 

The proceedings of the Higgs working group comprises a large number of contributions on a wide 
variety of topics. Roughly speaking, the contributions fall into one of three categories. 

The first category is theoretical calculations of Higgs production and decay processes, including 
higher-order corrections and resummation to all orders. There is an overview of Higgs total cross sec- 
tions, both in the Standard Model and with supersymmetiy. There is a review of calculations of Higgs 
production in association with heavy quarks, either bottom or top. In the case of Higgs production in 
association with bottom quarks, there is a discussion of the Higgs transverse momentum distribution, 
including the resummation of soft gluons, for both inclusive Higgs production as well as production in 
association with a high- Pt 6 jet. These calculations make use of the 6 distribution function in the proton, 
and there is a contribution regarding sets of parton distribution functions with no heavy quarks, with only 
c quai^ks, or with both c and 6 quaifc, at next-to-next-to-leading order in QCD. Finally, there is a calcu- 
lation of the electroweak corrections to Higgs production via gg — > h, which is the dominant production 
mechanism. 

The second category is non-standard Higgs bosons, either with or without supersymmetry. There 
is a discussion of the impact of radiative corrections on the search for supersymmetric Higgs bosons at the 
Tevatron and the LHC. There is an analysis of the search for a Higgs decaying via /i — > aa — > bbT~^T~ 
at the Tevatron, where a is also a Higgs scalar (or pseudoscalar). There is a discussion on how to use the 
processes 66 — > /i, /i — > r+r^, and /i — > 77 to disentangle the nature of electroweak symmetry breaking. 
Methods to search for a Higgs boson that decays invisibly are proposed. Finally, there is a discussion of 
the search for charged Higgs bosons at hadron colliders. 

The third category is experimental reviews, analyses, and developments. There are reviews from 
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both CDF and DO on the status and prospects for Higgs searches at the Tevatron. There ai^e studies 
on b jets, one on Z ^ 66 and the other on improving the 6-jet resolution. There are studies on h ^ 
W~^W~ i'^i'vu and h — > r+r^ at the LHC. There is a discussion of the diphoton background at 
the Tevatron, which is relevant to the search for the Higgs via /i ^ 77 at the LHC. 

All of these contributions represent real progress towards the elucidation of the mechanism of elec- 
troweak symmetry breaking. It will require the best efforts of us all to extract the maximal information 
from the data coming from the Tevatron and the LHC. 

Acknowledgment 

This material is based upon work supported by the National Science Foundation under Grant Number 
0547780. 
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2 SM and MSSM Higgs Boson Production Cross Sections 
Contributed by: T. Hahn, S. Heinemeyer, F. Maltoni, S. Willenbrock 

We present the SM and MSSM Higgs-boson production cross sections at the Tevatron and the 
LHC. The SM cross sections are a compilation of state-of-the-art theoretical predictions. The MSSM 
cross sections are obtained from the SM ones by means of an effective coupling approximation, as 
implemented in FeynHiggs. Numerical results have been obtained in four benchmark scenarios for two 
values of tan /3, tan /? = 5, 40. 

2.1 Introduction 

Deciphering the mechanism of electroweak symmetry breaking (EWSB) is one of the main quests of 
the high energy physics community. Electroweak precision data in combination with the direct top- 
quark mass measurement at the Tevatron have strongly constrained the range of possible scenarios and 
hinted to the existence of a light scalar particle [5]. Both in the standard model (SM) and in its minimal 
supersymmetric extensions (MSSM), the W and Z bosons and fermions acquire masses by coupling to 
the vacuum expectation value(s) of scalar SU(2) doublet(s), via the so-called Higgs mechanism. The 
common prediction of such models is the existence of at least one scalar state, the Higgs boson. Within 
the SM, LEP has put a lower bound on the Higgs mass, nih > 114 GeV [6], and has contributed to the 
indirect evidence that the Higgs boson should be relatively light with a 95% probability for its mass to 
be below 186 GeV [5]. In the MSSM the experimental lower bound for the mass of the lightest state is 
somewhat weaker, and internal consistency of the theory predicts an upper bound of 135 GeV [7, 8, 9]. 

If the Higgs sector is realized as implemented in the SM or the MSSM, at least one Higgs bo- 
son should be discovered at the Tevatron and/or at the LHC. Depending on the mass, there are various 
channels available where Higgs searches can be performed. The power of each signature depends on 
the production cross section, a, and the Higgs branching ratio into final state particles, such as lep- 
tons or 6-jets, the total yield of events being proportional to a- BR. In some golden channels, such as 
gg ^ h ^ Z^*^Z — > 4/x, a discovery will be straightfoward and mostly independent from our ability 
to predict signal and/or backgrounds. On the other hand, for coupling measurements or for searches in 
more difficult channels, such as tth ttbb associated production, precise predictions for both signal 
and backgrounds are mandatory. Within the MSSM such precise predictions for signal and backgrounds 
are necessary in order to relate the experimental results to the underlying SUSY parameters. 

The aim of this note is to collect up-to-date predictions for the most relevant signal cross sections, 
for both the SM and the MSSM. In Section 2.2 we collect the results of state-of-the-art calculations for 
the SM cross sections as a function of the Higgs mass. In Section 2.3 we present the MSSM cross sections 
for the neutral Higgs-bosons in four benchmark scenarios. These results are obtained by rescaling the 
SM cross sections presented in the previous sections, using an effective coupling approximation. 

2.2 SM Higgs production cross sections 

In this section we collect the predictions for the most important SM Higgs production processes at the 
Tevatron and at the LHC. The relevant cross sections are presented in Figs. 2.2.1 and 2.2.2 as function of 
the Higgs mass. The results refer to fully inclusive cross sections. No acceptance cuts or branching ratios 
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SM Higgs production 




Fig. 2.2.1: Higgs-boson production cross sections (fb) at the Tevatron (y/s — 1.96 TeV) for the most relevant production 
mechanims as a function of the Higgs-boson mass. Results for gg —* h, qq Vh, bb are at NNLO in the QCD 

expansion. Weak boson fusion {qq — > qqh) and ti associated production are at NLO accuracy. 

are applied^ . We do not consider here diffractive Higgs production, ^ ff©p [10, 11, 12, 13, 14]. 
For the discussion of this channel in the MSSM we refer to Ref. [15]. 

We do not aim here at a detailed discussion of the importance of each signature at the Tevatron or 
the LHC, but only at providing the most accurate and up-to-date theoretical predictions. To gauge the 
progress made in the last years, it is interesting to compare the accuracy of the results available in the 
year 2000, at the time of the Tevatron Higgs Working Group [4], with those shown here. All relevant 
cross sections ai^e now known at least one order better in the strong-couphng expansion, and in some 
cases also electroweak corrections are available. 

• gg ^ h + X: gluon fusion 

This process is known at NNLO in QCD [16, 17, 18] (in the large top-mass limit) and at NLO 
in QCD for a quark of an arbitrary mass circulating in the loop [19, 20]. Some N'^LO results 
have recently been obtained in Refs. [21, 22]. The NNLO results plotted here are from Ref. [23] 
and include soft-gluon resummation effects at NNLL. MRST2002 at NNLO has been used [24], 

'More details and data files can be found at maltoni.web.cem.ch/maltoni/TeV4LHC . 
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Fig. 2.2.2: Higgs-boson production cross sections (fb) at the LHC (^/s — 14 TeV) for tiie most relevant production mechanims 
as a function of the Higgs-boson mass. Results for gg h, qq ^ Vh, bb ^ h are at NNLO in the QCD expansion. Weak 
boson fusion (qq qqh) and ti associated production are at NLO accuracy. Single-top associated production (qb qth) is at 
LO. 

with the renormalization and factorization scales set equal to the Higgs-boson mass. The overall 
residual theoretical uncertainty is estimated to be around 10%. The uncertainties due to the large 
top mass limit approximation (beyond Higgs masses of 2 x mt) are difficult to estimate but expected 
to be relatively small. Differential results at NNLO are also available [25]. NLO (two-loop) EW 
corrections are known for Higgs masses below 2my/, [26, 27], and range between 5% and 8% of 
the lowest order term. These EW corrections, however, are not included in Figs. 2.2.1, 2.2.2, and 
they are also omitted in the MSSM evaluations below. The same holds for the recent corrections 
obtained in Refs. [21, 22]. 

• qq qqh + X: vector boson fusion 

This process is known at NLO in QCD [28, 29, 30]. Results plotted here have been obtained 
with MCFM[31]. Leading EW corrections are taken into account by using a{Mz) as the (square 
of the) electromagnetic coupling. The PDF used is CTEQ6M [32] and the renormalization and 
factorization scales ai^e set equal to the Higgs-boson mass. The theoretical uncertainty is rather 
small, less than 10%. 
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• qq ^ Vh + X\ W,Z associated production 

These processes are known at NNLO in the QCD expansion [33] and at NLO in the electroweak 
expansion [34]. The results plotted here have been obtained by the LH2003 Higgs working group 
by combining NNLO QCD and NLO EW corrections [35]. The PDF used is MRST2001 and the 
renormalization and factorization scales are set equal to the Higgs-vector-boson invariant mass. 
The residual theoretical uncertainty is rather small, less than 5%. 

• bb ^ h + X: bottom fusion 

This process is known at NNLO in QCD in the five-flavor scheme [36]. The cross section in the 
four-flavor scheme is known at NLO [37, 38]. Results obtained in the two schemes have been 
shown to be consistent [35, 39, 40]. The results plotted here are from Ref. [36]. MRST2002 at 
NNLO has been used, with the renormalization scale set equal to mh and the factorization scale 
set equal to mh/4:. For results with one final-state fe-quai^k at high-p^ we refer to Ref. [41, 39]. For 
results with two final-state 6-quarks at high-py we refer to Ref. [37, 38]. 

• qq^ 99 tih + X: tt associated production 

This process is known at NLO in QCD [42, 43, 44]. The results plotted here are from Ref. [44]. The 
PDF used is CTEQ6M and the renormalization and factorization scales are set equal to mt+mf, /2. 

• qb ^ qth : single-top associated production 

This process is known at LO in QCD [45]. The results plotted here (t-channel production, LHC 
only) are from Ref. [46]. The PDF used is CTEQ5L and the renormalization and factorization 
scales are set equal to the Higgs-boson mass. 

2.3 MSSM Higgs production cross sections 

The MSSM requires two Higgs doublets, resulting in five physical Higgs boson degrees of freedom. 
These are the light and heavy CP-even Higgs bosons, h and H, the CP-odd Higgs boson. A, and the 
charged Higgs boson, H^. The Higgs sector of the MSSM can be specified at lowest order in terms of 
Mz, Ma, and tan/3 = V2/V1, the ratio of the two Higgs vacuum expectation values. The masses of the 
CP-even neutral Higgs bosons and the charged Higgs boson can be calculated, including higher-order 
corrections, in terms of the other MSSM parameters. 

After the termination of LEP in the year 2000 (the final LEP results can be found in Refs. [6, 47]), 
the Higgs boson search has shifted to the Tevatron and will later be continued at the LHC. For these anal- 
yses and investigations a precise prediction of the Higgs boson masses, branching ratios and production 
cross sections in the various channels is necessary. 

Due to the large number of free parameters, a complete scan of the MSSM parameter space is too 
involved. Therefore the search results at LEP [47] and the Tevatron [48, 49, 50], as well as studies for 
the LHC [51] have been performed in several benchmark scenarios [52, 53, 54]. 

The code FeynHiggs [55, 7, 8] provides a precise calculation of the Higgs boson mass spectrum. 



10 



couplings and the decay widths^. This has now been supplemented by the evaluation of all relevant 
neutral Higgs boson production cross sections at the Tevatron and the LHC (and the corresponding three 
SM cross sections for both colliders with Mf^ = mh,TnH,TnA)- They are calculated by using the 
effective coupling approach, rescaling the SM result^ . 

In this section we will briefly describe the benchmark scenarios with their respective features. The 
effective coupling approach, used to obtain the production cross sections within FeynHiggs, is discussed. 
Results for the neutral Higgs production cross sections at the Tevatron and the LHC are presented within 
the benchmark scenarios for two values of tan [5, tan /? = 5, 40. 



2.4 The benchmark scenarios 

We start by recalling the four benchmark scenarios [53] suitable for the MSSM Higgs boson search at 
hadron colliders^. In these scenarios the values of the parameters of the i and h sector as well as the 
gaugino masses are fixed, while tan f3 and Ma are the parameters that are varied. Here we fix tan /3 to a 
low and a high value, tan /3 = 5, 40, but vary Ma- This also yields a variation of M^ and niH- 

In order to fix our notations, we list the conventions for the inputs from the scalar top and scalar 
bottom sector of the MSSM: the mass matrices in the basis of the current eigenstates ti, t/j and ^l, &ij 
are given by 

_ / M^+m2 + cos2/3(i-|4)M| m,X, \ 

1^ mtXt M^+m2 + |cos2/?4M| J ' ^^"^-'^ 

_ ( Ml+ml + cos2P{-\ + \sl)Ml m,X, \ 

' [ m,X, M?^ +ml-l cos 2/34m| J ' ^^"^-^^ 



where 



mtXt = mt{At — fj, cot (3), rribXi, = mi, {Ah — fi tan P). (2.4.3) 



Here At denotes the trilinear Higgs-stop coupling, Ai, denotes the Higgs-sbottom coupling, and jj, is the 
higgsino mass parameter. SU(2) gauge invariance leads to the relation 

Mi,=M-,^. (2.4.4) 

For the numerical evaluation, a convenient choice is 

Mj^ = M^^ = = M^^ =: MsusY. (2.4.5) 

The parameters in the i/b sector are defined here as on-shell parameters, see Ref. [56] for a discussion 
and a translation to DR parameters. The top-quark mass is taken to be mt = rn^^^ = 172.7 GeV [57]. 

• The m™^^ scenario: 

This scenario had been designed to obtain conservative tan/5 exclusion bounds [58]. The pa- 
rameters are chosen such that the maximum possible Higgs-boson mass as a function of tan /? 



^The code can be obtained from www.feynhiggs.de . 

'The inclusion of the charged Higgs production cross sections is planned for the near future. 

''in the course of this workshop they have been refined to cover wider parts of the MSSM parameter space relevant especially 
for heavy MSSM Higgs boson production [54]. 
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is obtained (for fixed AfsusY and mt, and Ma set to its maximal value, Ma = 1 TeV). The 
pai^ameters are^: 

^SUSY = 1 TeV, fj, = 200 GeV, M2 = 200 GeV, 

Xt = 2 MsusY At = At,mg = 0.8 MgusY • (2.4.6) 

• The no-mixing scenario: 

This benchmark scenario is associated with vanishing mixing in the i sector and with a higher 
SUSY mass scale as compared to the m™^^ scenario to increase the parameter space that avoids 
the LEP Higgs bounds: 

MsusY = 2 TeV, /i = 200 GeV, M2 = 200 GeV, 

Xt = 2 MsusY Ai, = At,mg = 0.8 MgusY • (2.4.7) 

• The gluophobic Higgs scenario: 

In this scenario the main production cross section for the light Higgs boson at the LHC, gg h, 
can strongly suppressed for a wide range of the Ma — tan /3-plane. This happens due to a can- 
cellation between the top quark and the stop quark loops in the production vertex (see Ref. [59]). 
This cancellation is more effective for small i masses and for relatively large values of the i mixing 
parameter, Xt. The partial width of the most relevant decay mode, T{h — > 77), is affected much 
less, since it is dominated by the W boson loop. The parameters are: 

MsusY = 350 GeV, fi = 300 GeV, M2 = 300 GeV, 

Xt = -750 GeV Ab = At,mg = 500 GeV . (2.4.8) 

• The small Ocff scenario: 

Besides the channel (7(7 — > /i — > 77 at the LHC, the other channels for light Higgs searches at 
the Tevatron and at the LHC mostly rely on the decays /i — > 66 and h — > r+r". Including 
Higgs-propagator corrections the couplings of the lightest Higgs boson to down-type fermions is ~ 
sin aes, where Ucs is the loop coiTccted mixing angle in the neutral CP-even Higgs sector. Thus, 
if Oes is small, the two main decay channels can be heavily suppressed in the MSSM compared to 
the SM case. Such a suppression occui's for large tan /3 and not too lai^ge Ma- The pai^ameters of 
this scenario are: 

MsusY = 800 GeV, fi = 2.5MsusY, M2 = 500 GeV, 

Xt = -1100 GeV, Ab = At, m-g = 500 GeV . (2.4.9) 

'As mentioned above, no external constraints are taken into account. In the minimal flavor violation scenario, better agree- 
ment with BR(6 s'y) constraints would be obtained for the other sign of Xt (called the "constrained m™^"" scenario [53]). 
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2.5 The effective coupling approximation 

We consider the following neutral Higgs production cross sections at the Tevatron and the LHC {(p 
denotes all neutral MSSM Higgs bosons, 4> = h,H, A): 

gg ^ </. + X, (2.5.10) 

qq ^ qqcp + X, (2.5.11) 

qq W/Z(j) + X, (2.5.12) 

bh cp + X, (2.5.13) 

gg.qq ^ tt> . (2.5.14) 

The MSSM cross sections have been obtained by rescaling the corresponding SM cross sections of 
Section 2.2 either with ratio of the corresponding MSSM decay with (of the inverse process) over the 
SM decay width, or with the square of the ratio of the corresponding couplings. More precisely, we apply 
the following factors: 

• 99 ^ (f' + X: 

^{(p 99)mssm ^2 5 15) 

r(</' 9g)sM 

We include the full one-loop result with SM QCD corrections. MSSM two-loop corrections [60] 
have been neglected. 

• qq ^ qq(t> + X: 

|ff0yy,MSSMp y^y^z. (2.5.16) 



We include the full set of Higgs propagator corrections in the effective couplings. 

qq W/Z(t) + X: 

' ' ",MSSM| 



\9<i>vv,i 



\94>vv. 



SM 



2 



V = W,Z. (2.5.17) 



We include the full set of Higgs propagator corrections in the effective couplings. 
• hb^cp + X: 

r((/) 66)mssm 



(2.5.18) 



T{4> bb)sM 

We include here one-loop SM QCD and SUSY QCD corrections, as well as the resummation of 
all terms of O {{as tan /?)"). 
• 99, qq tt(t): 

|g(/.M"MSSMl C9S1Q^ 
I _ 12 ' yZ.J.iy) 

\94)tt,SM\ 

where 9ff,tt,MSSM ^^d grptt.SM composed of a left- and a right-handed part. We include the full 
set of Higgs propagator corrections in the effective couplings. 

In the effective couplings introduced in eqs. (2.5.15)-(2.5.19) we have used the proper normaliza- 
tion of the external (on-shell) Higgs bosons as discussed in Ref. [61]. 

It should be noted that the effective coupling approximation as described above does not take into 
account the MSSM-specific dynamics of the production processes. The theoretical uncertainty in the 
predictions for the cross sections will therefore in general be somewhat larger than for the decay widths. 
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2.6 Results 

Results for the neutral Higgs production cross sections at the Tevatron and the LHC are presented within 
the four benchmark scenarios for two values of tan /?, tan /? = 5, 40, giving a total of eight plots for each 
collider. 

Figs. 2.6.3 and 2.6.4 show the results for the Tevatron, while Figs. 2.6.5 and 2.6.6 show the LHC 
results. In Fig. 2.6.3 (2.6.5) the Higgs production cross sections for the neutral MSSM Higgs bosons at 
the Tevatron (LHC) in the m™'*^ scenario (upper row) and the no-mixing scenario (lower row) can be 
found. Fig. 2.6.4 (2.6.6) depicts the same for the gluophobic Higgs scenario (upper row) and the small 
Oes scenario (lower row). 

For low Ma values the production cross section of the h and the A are. similar, while for large Ma 
the cross sections of H and A are very close. This effect is even more pronounced for large tan /3. 

The results presented in this paper have been obtained for the MSSM with real parameters, i.e. the 
CP-conserving case. They can can easily be extended via the effective coupling approximation to the 
case of non-vanishing complex phases (as implemented in FeynHiggs). 
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Fig. 2.6.3: Neutral Higgses production cross sections (fb) at the Tevatron, = 1.96 TeV for the most relevant production 
mechanims as a function of the Higgs-boson mass. Results are based on the SM cross sections and evaluated through an 
effective coupling approximation in the ra^^^ and no-mixing scenarios, for tan P — 5, 40. 
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Fig. 2.6.5: Neutral Higgses production cross sections (fb) at the LHC, = 14 TeV, for the most relevant production mechan- 
ims as a function of the Higgs-boson mass. Results are based on the SM cross sections and evaluated through an effective 
coupling approximation in the m™^" and no-mixing scenarios, for tan/3 = 5, 40. 
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Fig. 2.6.6: Same as Fig. 2.6.5, for the gluophobic Higgs and small a^s scenarios. 



18 



3 Towards understanding the nature of Electroweak Symmetry Breaking at the Tevatron and 
LHC 

Contributed by: A. Belyaev, A. Blum, S. Chivukula, E. H. Simmons 

PACS 14.80.Cp,11.30.Pb,11.15.Ex 

In this study we discuss how to extract information about physics beyond the Standard Model 
(SM) from searches for a light SM Higgs at Tevatron Run 11 and CERN LHC. We demonstrate that 
new (pseudo)scalar states predicted in both supersymmetric and dynamical models can have enhanced 
visibility in standard Higgs search channels, making them potentially discoverable at Tevatron Run 11 
and CERN LHC. We discuss the likely sizes of the enhancements in the various search channels for each 
model and identify the model features having the largest influence on the degree of enhancement. We 
compare the key signals for the non-standard scalars across models and also with expectations in the SM, 
to show how one could start to identify which state has actually been found. In particular, we suggest 
the likely mass reach of the Higgs search in pp/pp — > — > t~^t~ for each kind of non-standard scalar 
state and we demonstrate that pp / pp ^ H —>■ j'y may cleanly distinguish the scalars of supersymmetric 
models from those of dynamical models and shed the light on the pattern of Electroweak Symmetry 
Breaking. 

3.1 Introduction 

The origin of electroweak symmetry breaking remains unknown. While the Standai^d Model (SM) of 
particle physics is consistent with existing data, theoretical considerations suggest that this theory is only 
a low-energy effective theory and must be supplanted by a more complete description of the underlying 
physics at energies above those reached so far by experiment. 

If the Tevatron or LHC do find evidence for a new scalar state, it may not necessarily be the 
Standard Higgs. Many alternative models of electroweak symmetry breaking have spectra that include 
new scalar or pseudoscalar states whose masses could easily lie in the range to which Run 11 is sensitive. 
The new scalars tend to have cross-sections and branching fractions that differ from those of the SM 
Higgs. The potential exists for one of these scalars to be more visible in a standard search than the SM 
Higgs would be. 

Here we discuss how to extract information about non-Standard theories of electroweak symmetry 
breaking from searches for a light SM Higgs at Tevatron Run 11 and CERN LHC. Ref. [62] studied the 
potential of Tevatron Run 11 to augment its seaixh for the SM Higgs boson by considering the process 
99 hsM T^T~ . Authors determined what additional enhancement of scalar- production and 
branching rate, such as might be provided in a non-standai^d model like the MSSM, would enable a 
scalar to become visible in the r+r" channel alone at Tevatron Run 11. Similar work has been done for 
99 hussM T^T~ at the LHC [63] and for gg hsM 77 at the Tevatron [64] and LHC [65]. 

Our work builds on these results, considering an additional production mechanism (b-quark anni- 
hilation), more decay channels {bb, W~^W", ZZ, and 77), and a wider range of non-standard physics 
(supersymmetry and dynamical electroweak symmetry breaking) from which rate enhancement may de- 
rive. We discuss the possible sizes of the enhancements in the various search channels for each model 
and pinpoint the model features having the largest influence on the degree of enhancement. We suggest 
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the mass reach of the standai^d Higgs seaixhes for each kind of non-standard scalar state. We also com- 
pare the key signals for the non-standai^d scalai^s across models and also with expectations in the SM, 
to show how one could identify which state has actually been found. Analytic formulas for the decay 
widths of the SM Higgs boson are taken from [66], [67] and numerical values ai^e calculated using the 
HDECAY program [68]. 



3.2 Models of Electroweak Symmetry Breaking 

Supersymmetry 

One interesting possibility for addressing the hierarchy and triviality problems of the Standard Model is 
to introduce supersymmetry. 

In order to provide masses to both up-type and down-type quarks, and to ensure anomaly can- 
cellation, the minimal supersymmetric Standard Model (MSSM) contains two Higgs complex-doublet 
superfields: = and = {^t, ^u) which aquire two vacuum expectation values vi and V2 

respectively. Out of the original 8 degrees of freedom, 3 serve as Goldstone bosons, absorbed into longi- 
tudinal components of the and Z, making them massive. The other 5 degrees of freedom remain in 
the spectrum as distinct scalar states, namely two neutral CP-even states(/i, H), one neutral, CP-odd state 
(A) and a charged pair (H^). It is conventional to choose tan/3 = vi/v2 and Ma = — 
to define the SUSY Higgs sector. There are foloowing relations between Higgs masses which will be 
useful for determining when Higgs boson interactions with fermions are enhanced: 



{Ml + M|) ^ J {Ml + M|)2 - 4M2m| cos2 2/3 



2,^ , Ml{Ml-Ml) 
;cos^(/3-a)- ''^ ^ 



Ml{Ml-Miy 
(3.2.20) 

where a is the mixing angle of CP-even Higgs bosons. The Yukawa interactions of the Higgs fields with 
the quarks and leptons can be written as: ^ 

Y^,u/Y^i' = cos a/ sin /3 ^^mf/^^f ' = sin a/ sin [3 YAttlYm = cot /3 

yhb-b/ynbi = - sin a/ cos (3 Yh,i/Y^^' = cos a/ cos /3 Yj,,j,/Y^^f = tan ^B.2.21) 

relative to the Yukawa couplings of the Standard Model (Y^^j 'j = mj /v). Once again, the same pattern 
holds for the tau lepton's Yukawa couplings as for those of the b quark. There are several circumstances 
under which various Yukawa couplings ai^e enhanced relative to Standard Model values. For high tan /? 
(small cos /?), eqns. (3.2.21) show that the interactions of all neutral Higgs bosons with the down-type 
fermions are enhanced by a factor of 1/ cos /3. In the decoupling limit, where Ma — > oo, applying eq. 
(3.2.20) to eqns. (3.2.21) shows that the H and A Yukawa couplings to down-type fermions are enhanced 
by a factor of ~ tan /3. Conversely, for low niA — rn^, one can check that Y^^^^i/Y^^^^ = Yy^^^^jY^^J ~ 
tan /? that h and A Yukawas are enhanced instead. 



Technicolor 

Another intriguing class of theories, dynamical electroweak symmetry breaking (DEWSB), supposes 
that the scalar states involved in electroweak symmetry breaking could be manifestly composite at scales 

*Note that the interactions of the A are pseudoscalar, i.e. it couples to if'"f5'>lj- 
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not much above the electroweak scale v 250 GeV. In these theories, a new asymptotically free strong 
gauge interaction (technicolor [69, 70, 71]) breaks the chiral symmetries of massless fermions / at a 
scale A 1 TeV. If the fermions caiTy appropriate electroweak quantum numbers (e.g. left-hand (LH) 
weak doublets and right-hand (RH) weak singlets), the resulting condensate (/l/r) 7^ breaks the 
electroweak symmetry as desired. Thi^ee of the Nambu-Goldstone Bosons (technipions) of the chiral 
symmetry breaking become the longitudinal modes of the W and Z. The logarithmic running of the 
strong gauge coupling renders the low value of the electroweak scale natural. The absence of fundamental 
scalars obviates concerns about triviality. 

Many models of DEWSB have additional light neutral pseudo Nambu-Goldstone bosons which 
could potentially be accessible to a standard Higgs search; these are called "technipions" in technicolor 
models. Our analysis will assume, for simplicity, that the lightest PNGB state is significantly lighter than 
other neutral (pseudo) scalar technipions, so as to heighten the comparison to the SM Higgs boson. 

The specific models we examine are: 1) the traditional one-family model [72] with a full fam- 
ily of techniquarks and technileptons, 2) a variant on the one-family model [73] in which the lightest 
technipion contains only down-type technifermions and is significantly lighter than the other pseudo 
Nambu-Goldstone bosons, 3) a multiscale walking technicolor model [74] designed to reduce flavor- 
changing neutral currents, and 4) a low-scale technciolor model (the Technicolor Straw Man model) [75] 
with many weak doublets of technifermions, in which the second-lightest technipion P' is the state rele- 
vant for our study (the lightest, being composed of technileptons, lacks the anomalous coupling to gluons 
required for gg ^ P production). For simplicity the lightest relevant neutral technipion of each model 
will be generically denoted P; where a specific model is meant, a superscript will be used. 

One of the key differences among these models is the value of the technipion decay constant Fp, 
which is related to the number Nd of weak doublets of technifermions that contribute to electroweak 
symmetry breaking. We refer reader to [76] for details. 

3.3 Results For Each Model 

Supersymmetry 

Let us consider how the signal of a light Higgs boson could be changed in the MSSM, compared to 
expectations in the SM. There are several important sources of alterations in the predicted signal, some 
of which are interconnected. 

First, the MSSM includes three neutral Higgs bosons H = {h, H, A) states. The apparent signal 
of a single light Higgs could be enhanced if two or three neutral Higgs species are nearly degenerate, 
and we take advantage of this near-degeneracy by combining the signals of the different neutral Higgs 
bosons when their masses are closer than the experimental resolution. 

Second, the alterations of the couplings between Higgs bosons and ordinary fermions in the MSSM 
can change the Higgs decay widths and branching ratios relative to those in the SM. Radiative effects on 
the masses and couplings can substantially alter decay branching fractions in a non-universal way. For 
instance, B{h r + r~) could be enhanced by up to an order of magnitude due to the suppression 
of B{h — > bb) in certain regions of parameter space [77, 78]. However, this gain in branching fraction 
would be offset to some degree by a reduction in Higgs production through channels involving l^^^j [62]. 

Third, a large value of tan /? enhances the bottom-Higgs coupling (eqns. (3.2.21) ), making gluon 
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Fig. 3.3.7: NLO cross sections for Higgs production via the bb ^ Ti. and gg ^ TC processes (as well as their sum) at the 
Tevatron for the SM Higgs (a) the Supersymmetric axial Higgs boson with tan /3 = 30 (b). 



fusion through a 6-quark loop significant, and possibly even dominant over the top-quark loop contribu- 
tion. 

Fourth, the presence of superpartners in the MSSM gives rise to new squark-loop contributions to 
Higgs boson production through gluon fusion. Light squarks with masses of order 100 GeV have been 
argued to lead to a considerable universal enhancement (as much as a factor of five) [79, 80, 81, 82] for 
MSSM Higgs production compared to the SM. 

Finally, enhancement of the Y-^^f^i coupling at moderate to large tan /3 makes bb ^ H a. significant 
means of Higgs production in the MSSM - in contrast to the SM where it is negligible. To include both 
production channels when looking for a Higgs decaying a.sH ^ xx, we define a combined enhancement 
factor 



total I XX 



(^{99 — > W — > xx) + a{ 



7i — > xx) 



^'^gg/xx + ''^'bb/xx^bb:gg]/[^ + Rbb.gg] 



a{gg hsM ^ xx) + a{bb hsu ^ xx) 

(3.3.22) 

Here Rbb-.gg is the ratio of bb and gg initiated Higgs boson production in the Standard Model, which can 
be calculated using HDECAY. 

Figure 3.3.7 presents NLO cross sections at the Tevatron. For bb ^ H we are using the code 
of Ref. [83], ' while for gg ^ H we use HIGLU [84] and HDECAY [68] One can see that in the 
MSSM the contribution from bb ^ H becomes important even for moderate values of tan (3 10. For 
Mt-i < 110 — 115 GeV the contribution from gg ^ H process is a bit bigger than that from bb H, 



^Note that bb ^ H has been recently calculated at NNLO in [36]. 



Specifically, we use the HIGLU package to calculate the gg hern cross section. We then use the ratio of the Higgs decay 
widths from HDECAY (which includes a more complete set of one-loop MSSM corrections than HIGLU) to get the MSSM 



gg^H cross section: a^"^^^' = a^^' x r(H ^ gg)/r{hsM gg). 



22 



(a) gg+bb^A+H+h, tanp=30, Tevatron/LHC 



(b) gg+bb^A+H+h, tan^=50, Tevatron/LHC 
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Fig. 3.3.8: Enhancement factor n^f/^^ for final states xx = bb, t'^t~ , WW, ZZ, 77 when both gg —^ Ti. and bb ^ Ti. are 
included and the signals of all three MSSM Higgs states are combined. Frames (a) and (b) correspond to tan/3 — 30 and 50, 
respectively, at the Tevatron (solid lines) and at the LHC (dashed lines). 



while for M?^ > 115 GeV 6-quark-initiated production begins to outweigh gluon-initiated production. 
Results for LHC are qualitatively similar", except the rate, which is about two orders of magnitude higher 
compared to that at the Tevatron. 

Using the Higgs branching fractions with these NLO cross sections for gg ^ H and bb ^ H 
allows us to derive as presented in Fig. 3.3.8 for the Tevatron and LHC. There are several 

"physical" kinks and peaks in the enhancement factor for various Higgs boson final states related to 
WW, ZZ and top-quark thresholds which can be seen for the respective values of Ma- At very large 
values of tan f3 the top-quark threshold effect for the 77 enhancement factor is almost gone because the 
b-quark contribution dominates in the loop. One can see from Fig. 3.3.8 that the enhancement factors at 
the Tevatron and LHC aie very similar. On the other hand, the values of the total rates at the LHC are 
about two orders of magnitude higher than the corresponding rates at the Tevatron. In contrast to strongly 
enhanced bb and rf signatures, the 77 signature is always strongly suppressed! This particular feature of 
SUSY models, as we will see below, may be important for distinguishing supersymmetric models from 
models with dynamical symmetry breaking. 

It is important to note that combining the signal from the neutral Higgs bosons /i, A, H in the 
MSSM turns out to make our results more broadly applicable across SUSY parameter space. Combining 
the signals from A, h, H has the virtue of making the enhancement factor independent of the degree of 
top squark mixing (for fixed Ma, ^^ and Ms and medium to high values of tan /?), which greatly reduces 
the parameter-dependence of our results. 
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Table 3.3.1: Calculated enhancement factors for production at the Tevatron and LHC of a 130 GeV technipion via gg alone, via 
bb alone, and combined. Note that the small enhancement in the bb process slightly reduces the total enhancement relative to 
that of gg alone. In all cases, Ntc ~ 4. 
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Technicolor 

Single production of a technipion can occur through the axial-vector anomaly which couples the techni- 
pion to pairs of gauge bosons. For an SU (Ntc) technicolor group with technipion decay constant Fp, 
the anomalous coupling between the technipion and a pair of gauge bosons is given, in direct analogy 
with the coupling of a QCD pion to photons, by [85, 86, 87]. Comparing a PNGB to a SM Higgs boson 
of the same mass, we find the enhancement in the gluon fusion production is 

_ T{P ^ gg) _9 . .2 v"^ n 'X -yx^ 

Kgg prod - Y{h~^~^ ~ 4^^TC-^gg^ 

The main factors influencing Kgg p^od for a fixed value of Ntc are the anomalous coupling to 
gluons and the technipion decay constant. The value of Kgg pj-od for e,a.ch model (taking Ntc = 4) is 
given in Table 3.3. 

The value of ki,i, pj-od (shown in Table 3.3) is controlled by the size of the technipion decay con- 
stant. 

We see from Table 3.3 that prod is at least one order of magnitude smaller than Hgg prod in each 
model. From the Kgg prod/ >^bb prod ratio which reads as 



Kgg prod _ 9 ,^2 a2 \-2 / i 



2 \ 2 



we see that the larger size of Kgg prod is due to the factor of N^^, coming from the fact that gluons couple 
to a technipion via a techniquark loop. The extended technicolor (ETC) interactions coupling ft-quai^ks 
to a technipion have no such enhancement. With a smaller SM cross-section and a smaller enhancement 
factor, it is clear that technipion production via bb annihilation is essentially negligible at these hadron 
colliders. 

We now calculate the technipion branching ratios from the above information, taking Ntc = 4. 
The values are essentially independent of the size of Mp within the range 120 GeV - 160 GeV; the 
branching fractions for Mp = 130 GeV are shown in Table 3.3.2. The branching ratios for the SM Higgs 
at NLO are given for comparison; they were calculated using HDECAY [68]. Comparing the technicolor 
and SM branching ratios in Table 3.3.2, we see immediately that all decay enhancements. Model 2 is 
an exception; its unusual Yukawa couplings yield a decay enhancement in the t~^t~ channel of order 
the technipion's (low) production enhancement. In the 77 channel, the decay enhancement strongly 
depends on the group-theoretical sti'ucture of the model, through the anomaly factor. Our results for 
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Table 3.3.2: Branching ratios of Technipions/Higgs of mass 130 GeV 
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Table 3.3.3: Enhancement Factors for 130 GeV technipions produced at the Tevatron and LHC, compared to production and 
decay of a SM Higgs Boson of the same mass. The slight suppression of Kp^od due to the b-quark annihilation channel has 
been included. The rightmost column shows the cross-section (pb) for pp/pp P ^ xx at Tevatron Run II/LHC. 
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the Tevatron Run II and LHC production enhancements (including both gg fusion and bb annihilation), 
decay enhancements, and overall enhancements of each technicolor model relative to the SM are shown 
in Table 3.3.3 for a technipion or Higgs mass of 130 GeV. Multiplying k^^^^^^ by the cross-section for 
SM Higgs production via gluon fusion [84] yields an approximate technipion production cross-section, 
as shown in the right-most column of Table 3.3.3. 

In each technicolor model, the main enhancement of the possible technipion signal relative to that 
of an SM Higgs arises at production, making the size of the technipion decay constant the most critical 
factor in determining the degree of enhancement for fixed Ntc- 

3.4 Interpretation 

We are ready to put our results in context. The large QCD background for qq states of any flavor makes 
the tau-lepton-pair and di-photon final states the most promising for exclusion or discovery of the Higgs- 
like states of the MSSM or technicolor. We now illustrate how the size of the enhancement factors for 
these two final states vary over the parameter spaces of these theories at the Tevatron and LHC. We use 
this information to display the likely reach of each experiment in each of these standard Higgs search 
channels. Then, we compare the signatures of the MSSM Higgs bosons and the various technipions to 
see how one might tell these states apart from one another. 

In of Figure 3.4.9 we summarize the ability of Tevatron (left) and LHC (right) to explore the 
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gg+bb^ A+H+h -> zx, Tevatron, Vs = 1.96 TeV gg+bb^ A -> xx, LHC, = 14 TeV 




110 120 130 140 150 160 170 180 190 200 110 120 130 140 150 160 170 180 190 200 



(GeV) (GeV) 

Fig. 3.4.9: Results for gg-\-bb^h-\-H-\-A^ t^t~ at LHC. Left frame: Selected contours of given enhancement factor 
values kJ^j^,/^^ in the MSSM. Right frame: Predicted LHC reach, based on the hsM — > r+r^ studies of [63], in the MSSM 
parameter space. 

MSSM parameter space (in terms of both a 2cj exclusion curve and a 5cj discovery curve) using the 
process gg + bb^h + A + H^ t'^t". Translating the enhancement factors into this reach plot draws 
on the results of [62]. As the Ma mass increases up to about 140 GeV, the opening of the W~^W~ 
decay channel drives the t~^t~ branching fraction down, and increases the tan /3 value required to make 
Higgses visible in the r+r~ channel. At still larger M4, a very steep drop in the gluon luminosity 
(and the related 6-quark luminosity) at large x reduces the phase space for H production. Therefore 
for Ma >170 GeV, Higgs bosons would only be visible at very high values of tan (3. The pictures for 
tevaron and LHC are qualitatively similar, the main differences compared to the Tevatron are that the 
required value of tan /3 at the LHC is lower for a given Ma and it does not climb steeply for Ma >170 
GeV because there is much less phase space suppression. 

It is important to notice that both, Tevatron and LHC, could observe MSSM Higgs bosons in the 
r+T~ channel even for moderate values of tan /3 for Ma < 200 GeV, because of significant enhancement 
of this channel. However the 77 channel is so suppressed that even the LHC will not be able to observe 
it in any point of the Ma < 200 GeV parameter space studied in this paper! ^ 

The Figure 3.4. 10 presents the Tevatron and LHC potentials to observe technipions. For the Teva- 
tron, the observabiUty is presented in terms of enhancement factor, while for the LHC we present signal 
rate in term of u x Br{P tt/^^). At the Tevatron, the available enhancement is well above what is 
required to render the P of any of these models visible in the t~^t~ channel. Likewise, the right frame 
of that figure shows that in the 77 channel at the Tevatron the technipions of models 3 and 4 will be 
observable at the Scr level while model 2 is subject to exclusion at the 2a level. The situation at the LHC 

'in the decoupling limit with large values of Ma and low values of tan (3, the lightest MSSM Higgs could be dicovered in 
the 77 mode just like the SM model Higgs boson 
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is even more promising: all four models could be observable at the Scr level in both the r+r (left frame) 
and 77 (right frame) channels. 

Once a supposed light "Higgs boson" is observed in a collider experiment, an immediate important 
task will be to identify the new state more precisely, i.e. to discern "the meaning of Higgs" in this 
context. Comparison of the enhancement factors for different channels will aid in this task. Our study has 
shown that comparison of the t~^t~ and 77 channels can be particularly informative in distinguishing 
supersymmetric from dynamical models. In the case of supersymmetry, when the r+r" channel is 
enhanced, the 77 channel is suppressed, and this suppression is strong enough that even the LHC would 
not observe the 77 signature. In contrast, for the dynamical symmetry breaking models studied we expect 
simultaneous enhancement of both the r+r^ and 77 channels. The enhancement of the 77 channel is 
so significant, that even at the Tevatron we may observe technipions via this signature at the So" level 
for Models 3 and 4, while Model 2 could be excluded at 95% CL at the Tevatron. The LHC collider, 
which will have better sensitivity to the signatures under study, will be able to observe all four models 
of dynamical symmetry breaking studied here in the 77 channel, and can therefore distinguish more 
conclusively between the supersymmetric and dynamical models. 

3.5 Conclusions 

In this paper we have shown that searches for a light Standard Model Higgs boson at Tevatron Run II 
and CERN LHC have the power to provide significant information about important classes of physics 
beyond the Standard Model. We demonstrated that the new scalar and pseudo-scalar states predicted 
in both supersymmetric and dynamical models can have enhanced visibility in standard r+r" and 77 
search channels, making them potentially discoverable at both the Tevatron Run II and the CERN LHC. 
In comparing the key signals for the non-standai^d scalars across models we investigated the likely mass 
reach of the Higgs search in pp/pp H ^ t^t^ for each kind of non-standard scalar state, and we 
demonstrated that pp pp ^ H ^ may cleanly distinguish the scalars of supersymmetric models 
from those of dynamical models. 
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Fig. 3.4.10: Observability of technipions as a function of technipion mass and assuming the final state is a tau pair (left frame) 
or photon pair (right frame) at the Tveatron and LHC. Top raw: the 5cr discovery and 2a exclusion curves indicate the required 
enhancement factor for a Higgs-like particle at Tevatron Run II when the final state is r+r^ [62] (left frame) or 77 [64](right 
frame). Bottom raw: the lowest curve is the a x Br required to make a Higgs-like particle visible (5ct discovery) in t'^t~ [63] 
or in 77 [65] at LHC. 
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4 MSSM Higgs Boson Searches at the Tevatron and the LHC: Impact of Different Benchmark 
Scenarios 

Contributed by: M. Carena, S. Heinemeyer, C.E.M. Wagner, G. Weiglein 

The MSSM requires two Higgs doublets, resulting in five pliysical Higgs boson degrees of free- 
dom. These are the light and heavy CP-even Higgs bosons, h and H, the CP-odd Higgs boson. A, and 
the charged Higgs boson, H^. The Higgs sector of the MSSM can be specified at lowest order in terms 
of Mz, Ma, and tan/5 = V2/V1, the ratio of the two Higgs vacuum expectation values. The masses of 
the CP-even neutral Higgs bosons and the charged Higgs boson can be calculated, including higher-order 
corrections, in terms of the other MSSM parameters. 

After the termination of LEP in the year 2000 (the close-to-final LEP results can be found in 
Refs. [6, 47]), the Higgs boson search has shifted to the Tevatron and will later be continued at the 
LHC. Due to the large number of free parameters, a complete scan of the MSSM parameter space is too 
involved. Therefore the search results at LEP have been interpreted [47] in several benchmark scenar- 
ios [52, 53]. Current analyses at the Tevatron and investigations of the LHC [51] potential also have 
been performed in the scenarios proposed in Refs. [52, 53]. The m™**^ scenario has been used to obtain 
conservative bounds on tan (3 for fixed values of the top-quai^k mass and the scale of the supersymmetric 
particles [58]. These scenarios are, conceived to study particular cases of challenging and interesting 
phenomenology in the searches for the SM-like Higgs boson, i.e. mostly the light CP-even Higgs boson. 

The current searches at the Tevatron are not yet sensitive to a SM-like Higgs in the mass region 
allowed by the LEP exclusion bounds [6, 47]. On the other hand, scenarios with enhanced Higgs boson 
production cross sections can be probed already with the currently accumulated luminosity. Enhanced 
production cross sections can occur in particular for low Ma in combination with large tan /3 due to 
the enhanced couplings of the Higgs bosons to down-type fermions. The corresponding limits on the 
Higgs production cross section times branching ratio of the Higgs decay into down-type fermions can 
be interpreted in MSSM benchmark scenarios. Limits from Run II of the Tevatron have recently been 
published for the following channels [88, 89, 50] (here and in the following (j) denotes all three neutral 
MSSM Higgs bosons, 4> = h,H, A): 

bb(j), 4> ^ bb (with one additional tagged 6 jet), (4.0.25) 
pp ^ 4> ^ T^T~ (inclusive), (4.0.26) 
pp^tt^ H^W^ bb, TUr . (4.0.27) 

The obtained cross section limits have been interpreted in the m™'^^ and the no-mixing scenario with 
a value for the higgsino mass parameter of /i = —200 GeV [88] and = ±200 GeV [89]. In these 
scenarios for Ma ~ 100 GeV the limits on tan /3 are tan (3 < 50. 

Here we investigate the dependence of the CDF and DO exclusion bounds in the M^-tan /? plane 
on the parameters entering thr^ough the most relevant supersymmetric radiative corrections in the theoret- 
ical predictions for Higgs boson production and decay processes. We will show that the bounds obtained 
from the bbcp, (p —>■ bb channel depend very sensitively on the radiative connections affecting the relation 
between the bottom quark mass and the bottom Yukawa coupling. In the channels with r+T~ final states, 
on the other hand, compensations between large connections in the Higgs production and the Higgs decay 
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occur. In this context we investigate the impact of a lai^ge radiative correction in the gg ^ (p production 
process that had previously been omitted. 

In order to reflect the impact of the corrections to the bottom Yukawa coupling on the exclusion 
bounds we suggest to supplement the existing m™'^^ and no-mixing scenarios, mostly designed to search 
for the light CP-even MSSM Higgs boson, h, with additional values for the higgsino mass parameter //. 
In fact, varying the value and sign of //, while keeping fixed the values of the gluino mass and the common 
third generation squark mass parameter MsusY> demonstrates the effect of the radiative corrections on 
the production and decay processes. The scenarios discussed here are designed specifically to study the 
MSSM Higgs sector without assuming any particular soft supersymmetry-breaking scenario and taking 
into account constraints only from the Higgs boson sector itself. In particular, constraints from requiring 
the correct cold dark matter density, BR(6 — > 57) or {g — 2)^, which depend on other parameters of the 
theory, are not crucial in defining the Higgs boson sector, and may be avoided. 

We also study the non-standard MSSM Higgs boson search sensitivity at the LHC, focusing on 
the processes pp — > H/A + X, H/A t^t~ and pp —>■ tH^ + X, rvr, and stress the 

relevance of the proper inclusion of supersymmetric radiative corrections to the production cross sections 
and decay widths. We show the impact of these corrections by investigating the variation of the Higgs 
boson discovery reach in the benchmark scenarios for different values of /i. In particular, we discuss the 
resulting modification of the parameter region in which only the light CP-even MSSM Higgs boson can 
be detected at the LHC. 



4.1 Predictions for Higgs boson production and decay processes 

Notation and renormalization 

The tree-level values for the CP-even Higgs bosons of the MSSM, rrih and mu, are determined by tan j3, 
the CP-odd Higgs-boson mass Ma, and the Z boson mass Mz- The mass of the charged Higgs boson, 
m}j±, is given in terms of Ma and the W boson mass, Miy. Beyond the tree-level, the main correction 
to the Higgs boson masses stems from the t/t sector, and for large values of tan/3 also from the h/h 
sector. 

In order to fix our notations, we list the conventions for the inputs from the scalar top and scalar 
bottom sector of the MSSM: the mass matrices in the basis of the current eigenstates Il^Ir and 6^, bR 
are given by (modulo numerically small D-term contributions) 

* y mtXt + 

where 

mtXt = mt{At — fi cot (3), rribXi, = mi, {Ah — fj, tan P). (4.1.29) 

Here At denotes the trilinear Higgs-stop coupling, A^ denotes the Higgs-sbottom coupling, and fi is the 
higgsino mass parameter. SU(2) gauge invariance leads to the relation M^^ = M-^^ . For the numerical 
evaluation, a convenient choice is 

Mi^ = M-^^ = Mi^ = M^^ =: MsusY. (4.1.30) 



\ „ Ml +ml rribXb \ 

] , Ml= \ , . , (4.1.28) 

I ' \ mbXb Ml +ml r 



30 



The Higgs sector observables furthermore depend on the SU(2) gaugino mass pai^ameter, M2. The other 
gaugino mass pai-ameter, Mi, is usually fixed via the GUT relation Mi = f ^M2. At the two-loop level 
also the gluino mass, rrig, enters the predictions for the Higgs-boson masses. 

Corrections to the MSSM Higgs boson sector have been evaluated in several approaches. The 
status of the available corrections to the masses and mixing angles in the MSSM Higgs sector (with 
real parameters) can be summarized as follows. For the one-loop part, the complete result within the 
MSSM is known [90, 91]. The by far dominant one-loop contribution is the 0{at) term due to top and 
stop loops (at = /i^/(47r), ht being the top-quark Yukawa coupling). Concerning the two-loop effects, 
their computation is quite advanced and has now reached a stage such that all the presumably dominant 
contributions ai-e known [92, 93, 94, 95, 96, 7, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 
108, 109, 110]. The remaining theoretical uncertainty on the light CP-even Higgs boson mass has been 
estimated to be below ~ 3 GeV [8, 111]. The above calculations have been implemented into public 
codes. The program FeynHiggs [55, 112] is based on the results obtained in the Feynman-diagrammatic 
(FD) approach [7, 8, 1 10]. It includes all the above corrections. The code CPsuperH [1 13] is based on the 
renormalization group (RG) improved effective potential approach [93, 94, 95, 96, 114]. For the MSSM 
with real parameters the two codes can differ by up to ~ 4 GeV for the light CP-even Higgs boson mass, 
mostly due to formally subleading two-loop corrections that are included only in FeynHiggs. 

It should be noted in this context that the FD result has been obtained in the on-shell (OS) renor- 
malization scheme, whereas the RG result has been calculated using the MS scheme; see Refs. [1 14, 115] 
for a detailed comparison. Owing to the different schemes used in the FD and the RG approach for the 
renormalization in the scalar top sector, the parameters Xt and MgusY are also scheme-dependent in the 
two approaches. 



Leading effects from the bottom/sbottom sector 

The relation between the bottom-quark mass and the Yukawa coupling hh, which controls also the inter- 
action between the Higgs fields and the sbottom quarks, is affected at one-loop order by large radiative 
corrections [105, 106, 107, 108, 109]. The leading effects are included in the effective Lagrangian for- 
malism developed in Ref. [108]. Numerically this is by far the dominant part of the contributions from 
the sbottom sector (see also Refs. [103, 110, 104]). The effective Lagrangian is given by 



eMw 00 + 



tan p Aibj5b + y/2Vtb tan(3 H+tLbR (4.1.31) 



, sma ^ cosa\ - , /cos a . sma , , 

+ o - ^b—-T hhLbR - + A;,—— HbLbR 

^ cos p sm p J \ cos p sm p 



+ h.c. 



Here mj, denotes the running bottom quark mass including SM QCD corrections. In the numerical 
evaluations obtained with FeynHiggs below we choose m;, = mf,{mt) ~ 2.97 GeV. The prefactor 
1/(1 + A5) in Equation 4. 1.32 arises from the resummation of the leading corrections to all orders. The 
additional terms ~ Af, in the hbb and Hbb couplings arise from the mixing and coupling of the "other" 
Higgs boson, H and h, respectively, to the b quarks. 

As explained above, the function A;, consists of two main contributions, an 0{as) connection from 
a sbottom-gluino loop and an 0{at) connection from a stop-higgsino loop. The explicit form of A^ in 



31 



the limit of Ms > mt and tan /? » 1 reads [105, 106, 107] 

Ab = m-g tan (3 x /(m^^ , m^^, m^) + — A* ^ tan/3 x I{m^^,m^^, ^) . (4.1.32) 
The function / is given by 

liaAc) = (,2 _ _ ,2) {-'b' log ^ + ^^^^ log g + c^a^ log g)(4.1.33) 

1 

max(a2, b'^,c'^) 

The large 6 — ^ loops are resummed to all orders of (ostan/?)" via the inclusion of A;, [105, 106, 
107, 108, 109]. The leading electroweak contributions are taken into account via the second term in 
Equation 4.1.32. 

For large values of tan /3 and the ratios of M"^s/-^suSY ^^'^ f^^t/M^^^y, the Af, correction can 
become very important. Considering positve values of At and nig, the sign of the Af, term is governed 
by the sign of fi. Cancellations can occur if At and nig have opposite signs. For nig, At > the 
Af, connection is positive, leading to a suppression of the bottom Yukawa coupling. On the other hand, 
for negative values of A^, the bottom Yukawa couphng may be strongly enhanced and can even acquire 
non-perturbative values when A;, — > —1. 

Impact on Higgs production and decay at large tan P 

Higgs-boson production and decay processes at the Tevatron and the LHC can be affected by different 
kinds of lai^ge radiative connections. For large tan (3 the supersymmetric radiative corrections to the 
bottom Yukawa coupling described above become particularly important [78, 77]. Their main effect on 
the Higgs-boson production and decay processes can be understood from the way the leading contribution 
Af, enters. In the following we present simple analytic approximation formulae for the most relevant 
Higgs-boson production and decay processes. They are meant for illustration only so that the impact 
of the Afe conections can easily be traced. In our numerical analysis below, we use the full result from 
FeynHiggs rather than the simple formulae presented in this section. No relevant modification to these 
results would be obtained using CPsuperH. 

We begin with a simple approximate formula that represents well the MSSM parametric variation 
of the decay rate of the CP-odd Higgs boson in the large tan /? regime. One should recall, for that 
purpose, that in this regime the CP-odd Higgs boson decays mainly into r-leptons and bottom-quarks, 
and that the partial decay widths are proportional to the square of the Yukawa couplings evaluated at an 
energy scale of about the Higgs boson mass. Moreover, for Higgs boson masses of the order of 100 GeV, 
the approximate relations m;,(M^)^ ~ 9 GeV^, and mT-(Afyi)^ ~ 3 GeV^ hold. Hence, since the number 
of colors is Nc = 3, for heavy supersymmetric particles, with masses far above the Higgs boson mass 
scale, one has 

BR{A ^ 66) ~ ^ , BR(^ ^ r+r") ~ ^ . (4.1.34) 

(1 + A,)V9' (i + A,f + 9 

On the other hand, the production cross section for a CP-odd Higgs boson produced in association 
with a pair of bottom quarks is proportional to the square of the bottom Yukawa coupling and therefore 
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is proportional to tan^ [5/ (1 + Af,)^. Also in the gluon fusion channel, the dominant contribution in the 
large tan (3 regime is governed by the bottom quark loops, and therefore is also proportional to the square 
of the bottom Yukawa coupling. Hence, the total production rate of bottom quarks and r pairs mediated 
by the production of a CP-odd Higgs boson in the large tan /3 regime is approximately given by 

aihlA) X BR(A ^ hh) ~ <yihhA)su % x o , (4.1.35) 

(1 + Ab)2 (i + A,)2 + 9 

tan^ 3 

aigg,bb^A)xBR{A^T+T') ~ a{gg,bb ^ A)sm ^ X n ' (4.1.36) 

(1 + Ai) +9 

where a{bbA)sM and (T{gg, bb — > ^)sm denote the values of the corresponding SM Higgs boson pro- 
duction cross sections for a Higgs boson mass equal to Ma- 

As a consequence, the bb production rate depends sensitively on A;, because of the factor 1/(1 + 
Afo)^, while this leading dependence on A;, cancels out in the t^t^ production rate. There is still a 
subdominant parametric dependence in the r+r" production rate on A^ that may lead to variations of a 
few tens of percent of the r-pair production rate (compared to variations of the rate by up to factors of a 
few in the case of bottom-quark pair production). 

The formulae above apply, within a good approximation, also to the non-standard CP-even Higgs 
boson in the large tan P regime. Depending on Ma this can be either the h (for Ma ^120 GeV) or the 
H (for Ma ^ 120 GeV). This non-standard Higgs boson becomes degenerate in mass with the CP-odd 
Higgs scalar. Therefore, the production and decay rates of H (h) are governed by similar formulae as 
the ones presented above, leading to an approximate enhancement of a factor 2 of the production rates 
with respect to the ones that would be obtained in the case of the single production of the CP-odd Higgs 
boson as given in Equations 4.1.35, (4.1.36). 

We now turn to the production and decay processes of the charged Higgs boson. In the MSSM, 
the masses and couplings of the charged Higgs boson in the large tan /3 regime are closely related to 
the ones of the CP-odd Higgs boson. The tree-level relation rnjj^ = M\ + receives sizable 
connections for large values of tan/3, /i. At and A}). These corrections depend on the ratios ^^/Mgugy' 
- ^6At)VM|usY> {^t + Ab)VAf|usY [93, 94, 95, 96]. The couphng of the charged Higgs boson 
to a top and a bottom quark at lai^ge values of tan (5 is governed by the bottom Yukawa coupling and is 
therefore affected by the same Af, corrections that appear in the couplings of the non-standard neutral 
MSSM Higgs bosons [108]. 

The relevant channels for charged Higgs boson searches depend on its mass. For values of 
mj/i smaller than the top-quark mass, searches at hadron colliders concentrate on the possible emis- 
sion of the charged Higgs boson from top-quark decays. In this case, for large values of tan/3, the 
charged Higgs decays predominantly into a r lepton and a neutrino, i.e. one has to a good approximation 
BR(iJ^ — > tVt) ~ 1- The partial decay width of the top quark into a charged Higgs and a bottom quark 
is proportional to the square of the bottom Yukawa coupling and therefore scales with tan^ /?/(l + A;,)^, 
see e.g. Ref. [108]. 

For values of the charged Higgs mass larger than mt, instead, the most efficient production channel 
is the one of a charged Higgs associated with a top quark (mediated, for instance, by gluon-bottom 
fusion) [116]. In this case, the production cross section is proportional to the square of the bottom-quark 
Yukawa coupling. The branching ratio of the charged Higgs decay into a r lepton and a neutrino is. 
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apart from threshold coiTections, governed by a similar formula as the branching ratio of the decay of 
the CP-odd Higgs boson into r-pairs, namely 

where the factor (1 — r^)^ is associated with threshold corrections, and rj = mf/m|^±. 

As mentioned above, our numerical analysis will be based on the complete expressions for the 
Higgs couplings rather than on the simple approximation formulae given in this section. 



4.2 Interpretation of cross section limits in MSSM 
scenarios 

Limits at the Tevatron 

The DO and CDF Collaborations have recently published cross section limits from the Higgs search at 
the Tevatron in the channel where at least three bottom quarks are identified in the final state (660, (j) 
66) [88] and in the inclusive channel with T+r~ final states {pp cj) —>■ r+r") [89]. The CDF Col- 
laboration has also done analyses seaiching for a charged Higgs boson in top-quark decays [50]. While 
the cross section for a SM Higgs boson is significantly below the above limits, a large enhancement of 
these cross sections is possible in the MSSM. It is therefore of interest to interpret the cross section lim- 
its within the MSSM parameter space. One usually displays the limits in the M^i-tan (3 plane. As the 
whole structure of the MSSM enters via radiative corrections, the limits in the Myi-tan /3 plane depend 
on the other parameters of the model. One usually chooses certain benchmark scenarios to fix the other 
MSSM parameters [52, 53]. In order to understand the physical meaning of the exclusion bounds in the 
M4-tan /? plane it is important to investigate how sensitively they depend on the values of the other 
MSSM parameters, i.e. on the choice of the benchmark scenarios. 

Limits from the process bb(j), (j) ^ hh 

The DO Collaboration has presented the limits in the M4-tan j3 plane obtained from the hhcf), (/> — > 66 
channel for the m™^^ and no-mixing scenarios as defined in Ref. [52]. The m™^^ scenario according to 
the definition of Ref. [52] reads 

mt = 174.3 GeV, MsusY = 1000 GeV, = -200 GeV, Ma = 200 GeV, 

X^^ = 2MsusY (FD calculation), = VqMsvsy (RG calculation), 

A = At, m-g = 0.8 MsusY . (4.2.38) 

The no-mixing scenario defined in Ref. [52] differs from the mf^'^^ scenario only in 

Xt = Qi (FD/RG calculation) . (4.2.39) 

The condition = At implies that the different mixing in the stop sector gives rise to a difference 
between the two scenarios also in the sbottom sector. The definition of the ni™^^ and no-mixing scenarios 
given in Ref. [52] was later updated in Ref. [53], see the discussion below. 
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For their analysis, the DO Collaboration has used the following approximate formula [88], 



a{bb(l)) X BR{(p bb) = 2 a{bb(l))sM 



tan^ (3 9 



(4.2.40) 



which follows from Equation 4.1.35 and the discussion in Section 4.1. The cross section a{bb(l))sM has 
been evaluated with the code of Ref. [41], while Af, has been calculated using CPsuperH [1 13]. From the 
discussion in Section 4.1 it follows that the choice of negative values of fi leads to an enhancement of the 
bottom Yukawa coupling and therefore to an enhancement of the signal cross section in Equation 4.2.40. 
For tan f] = 50 the quantity Af, takes on the following values in the m™^^ and no-mixing scenarios as 
defined in Equations 4.2.38, (4.2.39), 



no-mixing seen., = -200 GeV, MsusY = 1000 GeV, tan (3 = 50 : Ab = -0.10 . (4.2.42) 

While the 0{as) contribution to A;,, see Equation 4.1.32, is practically the same in the two scenarios, 
the 0{at) contribution to Af, in the m™**^ scenario differs significantly from the one in the no-mixing 
scenario. In the m™^^ scenario the 0{at) contribution to A;, is about as large as the 0{as) contribution. 
In the no-mixing scenario, on the other hand, the 0{at) contribution to Af, is very small, because At is 
close to zero in this case. Reversing the sign of fi in Equations 4.2.41, (4.2.42) reverses the sign of A;,, 
leading therefore to a significant suppression of the signal cross section in Equation 4.2.40 for the same 
values of the other MSSM parameters. 

The predictions for btxp, (j) ^ bb evaluated with FeynHiggs have been compared with the exclusion 
bound for a x BR as given in Ref. [88]. As mentioned above, in our analysis we use the full Higgs cou- 
plings obtained with FeynHiggs rather than the approximate formula given in Equation 4.2.40. Similar 
results would be obtained with CPsuperH. 

The impact on the limits in the M^-tan (3 plane from varying /x while keeping all other parameters 
fixed can easily be read off from Equation 4.2.40. For a given value of the CP-odd mass and tan (3, the 
bound on a{bb(j)) x BR((/) —>■ bb) provides an upper bound on the bottom-quark Yukawa coupling. The 
main effect therefore is that as ji varies, the bound on tan (3 also changes in such a way that the value of 
the bottom Yukawa coupling at the boundary line in the A/^i-tan (3 plane remains the same. 

The dependence of the limits in the M^-tan [3 plane obtained from the process bb(j), (p ^ bb on 
the parameter fi is shown in Figure 4.2.11. The limits for /x = —200 GeV in the m™^^ and no-mixing 
scenarios, corresponding to the limits presented by the DO Collaboration in Ref. [88], are compared with 
the Umits arising for different fi values, fi = +200, ±500, ±1000 GeV. Figure 4.2.11 illustrates that the 
effect of changing the sign of /i on the limits in the M^-tan (3 plane obtained from the process bb<p, (p 
bb is quite dramatic. In the mf^^^ scenario the exclusion bound degrades from about tan /3 = 50 for 
Ma = 90 GeV in the case of jJL = -200 GeV to about tan /3 = 90 for Ma = 90 GeV in the case of 
/i = ±200 GeV. We extend our plots to values of tan [3 much larger than 50 mainly for illustration 
purposes; the region tan /3 » 50 in the MSSM is theoretically disfavoured, if one demands that the 
values of the bottom and r Yukawa couplings remain in the perturbative regime up to energies of the 
order of the unification scale. The situation for the bottom- Yukawa coupling can be ameliorated for large 
positive values of /i due to the A;, coiTcctions. The curves for n = ±500, ±1000 GeV do not appear- 
in the plot for the m™'^^ scenario, since for these /z values there is no tan (3 exclusion below tan (3 = 



^max scenario, = -200 GeV, tan [3 = 50 



0.21 ,(4.2.41) 
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Fig. 4.2.11: Change in the limits obtained from the bbcp, (j> ^ bb channel in the m™^" (left) and no-mixing (right) benchmark 
scenarios for different values of fj,. The value jj, = —200 GeV was chosen by the DO Collaboration in Ref. [88]. The other 
curves indicate the corresponding limits for ^ = +200, ±500, ±1000 GeV. The curves for /j, — ±500, ±1000 GeV (^ — 
±1000 GeV) do not appear in the left (right) plot for the m™^" (no-mixing) scenario, since for these fi values there is no tan (3 
exclusion below tan /3 = 130 for any value of Ma- 

130 for any value of Ma- On the other hand, the large negative values of fi shown in Figure 4.2.11, 
/i = —500, —1000 GeV, lead to an even stronger enhancement of the signal cross section than for 
fi = —200 GeV and, accordingly, to an improved reach in tan/3. It should be noted that for fi = 
—500, —1000 GeV the bottom Yukawa coupling becomes so large for tan/? » 50 that a perturbative 
treatment would no longer be reUable in this region. 

In Ref. [53] the definition of the m^^^ and no-mixing scenarios given in Ref. [52] has been up- 
dated. The sign of fi in the m™^^ and no-mixing scenarios has been reversed to /i = ±200 GeV in 
Ref. [53]. This leads typically to a better agreement with the constraints from {g — 2)^. Furthermore, 
the value of AfsusY in the no-mixing scenario was increased from 1000 GeV [52] to 2000 GeV in or- 
der to ensure that most of the parameter space of this scenario is in accordance with the LEP exclusion 
bounds [6, 47]. 

Another scenario defined in Ref. [53] is the "constrained-m™'^^" scenario. It differs from the 



mi 



scenario as specified in Ref. [53] by the reversed sign of Xt, 



OS 



-2MsusY (FDcalc), X, 



MS 



Ve MsusY (RG calc), /i = +200 GeV . (4.2.43) 



For small Ma and minimal flavor violation this results in better agreement with the constraints from 
BR(6 — > sj). For large tan/3 one has At ~ Xt, thus At and rrig have opposite signs. This can lead 
to cancellations in the two contributions entering Af,, see Equation 4.1.32. In contrast to the m™'^^ 
scenario, where the two contributions entering A;, add up, see Equation 4.2.41, the constrained-m™'^^ 
scenario typically yields relatively small values of A^ and therefore a coiTcspondingly smaller effect on 
the relation between the bottom-quark mass and the bottom Yukawa coupling, e.g. 

constrained-m)^'"^ scenario, /x = +200 GeV, tan (3 = 50 : Ab = -0.001 . (4.2.44) 

For large values of |;u| the compensations between the two terms entering A;, are less efficient, since the 
function / in the second term of Equation 4.2.41 scales Uke 1///^ for large |//|. 
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The impact of the benchmark definitions of Ref. [53] on the limits in the M^-tan /3 plane aiising 
from the bb(j), (j) ^ bb channel has been analyzed in Ref. [54]. The effect of changing AfsusY = 
1000 GeV to MsusY = 2000 GeV in the no-mixing scenario for /i = ±200 GeV results in substantially 
weaker (stronger) limits for ;U = +(— )200 GeV Also the constrained-jn^^^ scenario has been analyzed 
in Ref. [54]. As expected the variation of the exclusion bounds with a variation of /U is much weaker than 
in the other scenarios. 

Limits from the process pp ^ <j) ^ t^t~ 

The limits obtained from the — > — > r+r~ channel by the CDF Collaboration were presented in the 
Myi-tan [3 plane for the m™^^ and no-mixing scenarios as defined in Ref. [53] and employing two values 
of the /i parameter, ^ = ±200 GeV. According to the discussion in Section 4. 1, the limits obtained from 
the pp ^ (p ^ T~^T~ channel are expected to show a weaker dependence on the sign and absolute 
value of fi than the limits arising from the bb(p, (p ^ bb channel. On the other hand, for large values of 
tan /3 and negative values of fi, the large corrections to the bottom Yukawa coupling discussed above can 
invalidate a perturbative treatment for this channel. 

The MSSM prediction for a{pp (/)) x BK{(j) t^t^) as a function of tan/3 has been eval- 
uated by the CDF collaboration using the HIGLU program [84] for the gluon fusion channel. The pre- 
diction for the bb ^ (j) + X channel was obtained from the NNLO result in the SM from Ref. [36], and 
[a X BR]j^ggj^ / W >^ ^I^lsM calculated with the FeynHiggs program [55, 112]. While the full 
connection to the bottom Yukawa coiTcction was taken into account in the bb ^ 4> + X production chan- 
nel and the (f) t~^t~ branching ratios, the public version of the HIGLU program [84] does not include 
the Ab connection for the bottom Yukawa coupling entering the bottom loop contribution to the gg ^ (j) 
production process. In order to treat the two contributing production processes in a uniform way, the A;, 
conection should be included (taking into account the 0{as) and the 0{at) parts, see Equation 4.1.32) 
in the gg ^ (p production process calculation. For the large value of MgusY chosen in the m™^^ and 
no-mixing benchmark scenarios other higher-order contributions involving sbottoms and stops can be 
neglected (these effects are small provided MsusY ^ 500 GeV). 

In Ref. [54] a comparison of the "partial A^" and the "full Aj," results has been performed. It was 
shown that the inclusion of the A;, corrections everywhere can lead to a variation of A tan/? ~ 10 in 
the m™'^^ scenario, but has a much smaller effect in the no-mixing scenario. Following our analysis, the 
CDF Collaboration has adopted the prescription outlined above for incorporating the A^ conection into 
the gg ^ (j) production process. The limits given in Ref. [89] are based on the MSSM prediction where 
the Afe conection is included everywhere in the production and decay processes (see e.g. Ref. [117] for 
a previous analysis). 

We next turn to the discussion of the sensitivity of the limits obtained from the — > — > T~^r~ 
channel (including the Af, conection in all production and decay processes) on the sign and absolute 
value of ^. As discussed above, similar variations in the exclusion limits will occur if the absolute values 
of ^, rrig, At and MsusY are varied, while keeping the ratios appearing in A;, constant. The results 
are given in Figure 4.2.12 for the m™^^ scenario (left) and the no-mixing scenario (right). In the m™'^^ 
scenario we find a sizable dependence of the tan/3 bounds on the sign and absolute value of fi.^^ The 

'"For jj, — —300 GeV the curve stops at around tan/3 — 75 because the bottom Yukawa coupUng becomes very large, 
leading to instabilities in the calculation of the Higgs properties. For the same reason, even more negative values of /i are not 
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effect grows with Ma and, for the range of parameters explored in Figure 4.2. 12, leads to a variation of 
the tan /3 bound larger than A tan j3 ~ 30. In the no-mixing scenario the effect is again smaller, but it 
can still lead to a variation of the tan (3 bounds by as much as A tan [5 ~ 10. 
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Fig. 4.2.12: Variation of the limits obtained from the pp — > — > r r channel at the Tevatron in the m™™ (left) and 
no-mixing (right) benchmark scenarios for different values of ^. 

The results obtained in the constrained-m™'^^ scenario are again very robust with respect to varying 
/i. All values of /i result practically in the same tan (3 exclusion bounds [54]. 

Prospects for Higgs sensitivities at the LHC 

The most sensitive channels for detecting heavy MSSM Higgs bosons at the LHC are the channel pp 
H/A + X, H/A T~^T~ (making use of different decay modes of the two r leptons) and the channel 
tH^,H^ — > tUt (for mfj± > mt) [118, 119]. We consider here the pai^ameter region Ma ^ Mz, 
for which the heavy states H, A are widely separated in mass from the light CP-even Higgs boson h. 
Here and in the following we do not discuss search channels where the heavy Higgs bosons decay into 
supersymmetric particles, which depend very sensitively on the model parameters [120, 121, 119], but 
we will comment below on how these decays can affect the searches with bottom-quarks and r-leptons 
in the final state. 

Discovery region for the process pp H/A + X, H/A t^t~ 

To be specific, we concentrate in this section on the analysis carried out by the CMS Collaboration [122, 
119]. Similar results for this channel have also been obtained by the ATLAS Collaboration [118, 123, 
124]. In order to rescale the SM cross sections and branching ratios, the CMS Collaboration has used 
for the branching ratios the HDECAY program [68] and for the production cross sections the HIGLU 
program [84] {gg H/A) and the HQQ program [125] {gg hhH). In the HDECAY program the A;, 
corrections are partially included for the decays of the neutral Higgs bosons (only the 0{as) contribution 
to At is included, see Equation 4.1.32). The HIGLU program (see also the discussion in Section 4.2) and 

considered here. 
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HQQ, on the other hand, do not take into account the coiTcctions to the bottom Yukawa coupling. ' ^ The 
prospective 5a discovery contours for CMS (corresponding to the upper bound of the LHC "wedge" 
region, where only the light CP-even Higgs boson may be observed at the LHC) have been presented 
in Refs. [122, 119] in the My^-tan/? plane, for an integrated luminosity of 30 fb~^ and 60 fb~^. The 
results were presented in the m™*^^ scenario and for different fi values, /i = —200, +300, +500 GeV. 
It should be noted that decays of heavy Higgs bosons into charginos and neutralinos open up for small 
enough values of the soft supersymmetry-breaking parameters M2 and fi. Indeed, the results presented in 
Refs. [122, 1 19] show a degradation of the discovery reach in the M^-tan (3 plane for smaller absolute 
values of fi, which is due to an enhanced branching ratio of H, A into supersymmetric particles, and 
accordingly a reduced branching ratio into r pairs. 

We shall now study the impact of including the A;, corrections into the production cross sections 
and branching ratios for different values of fx. The inclusion of the A;, corrections leads to a modification 
of the dependence of the production cross section on tan (3, as well as of the branching ratios of the Higgs 
boson decays into t~^t~. For a fixed value of Ma, the results obtained by the CMS Collaboration for 
the discovery region in tan (3 can be interpreted in terms of a cross section limit using the approximation 
of rescaling the SM rate for the pp —>■ H + X, H —>■ t^t~ process by the factor 

tan^ 0OUS X ^n"^-^r-)cu, + mA^r + r-)cu. , 

BR{H^T + T-)sM 

In the above, tan /3cMS refers to the value of tan f3 on the discovery contour (for a given value of Ma) 
that was obtained in the analysis of the CMS Collaboration with 30 fb~^ [119]. These tan/3 values 
as a function of Ma correspond to the edge of the area in the M^-tan /? plane in which the signal 
pp — > H/A + X, H/A T'^T" is visible (i.e. the upper bound of the LHC wedge region). The 
branching ratios BR{H — > r + r~)cMS and BR(A t + t~)cms in the CMS analysis have been 
evaluated with HDECAY, incorporating therefore only the gluino-sbottom contribution to A^. 

After including all A;, corrections, we evaluate the pp H/A + X, H/A t^t^ process by 
rescaling the SM rate with the new factor, 

X ^ f- , (4.2.46) 



(1 + A,)2 BR(i7^r + r-)sM 

where A;, depends on tan /?. The quantities have been evaluted with FeynHiggs, allowing also decays 
into supersymmetric particles. The resulting shift in the discovery reach for the pp —>■ H/A+X, H/A —>■ 
T'^T" channel can be obtained by demanding that Equation 4.2.45 and Equation 4.2.46 should give the 
same numerical result for a given value of Ma- 

This procedure has been carried out in two benchmark scenarios for various values of ji. The 
results are shown in Figure 4.2A3 for the m™'*^ scenario (left) and for the no-mixing scenario (right). 
The comparison of these results with the ones obtained by the CMS Collaboration [122, 119] shows 
that for positive values of fi the inclusion of the supersymmetric radiative corrections leads to a slight 
shift of the discovery region towards higher values of tan (3, i.e. to a small increase of the LHC wedge 
region. For = —200 GeV the result remains approximately the same as the one obtained by the CMS 



"since HQQ is a leading-order program, non-negligible changes can also be expected from SM-QCD type higher-order 
corrections. 
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Collaboration. Due to the smaller considered tan /? values compared to the analysis of the Tevatron 
limits in Section 4.2, the connections to the bottom Yukawa coupling from A;, are smaller, leading to a 
better perturbative behavior. As a consequence, also the curves for fi = —500, —1000 GeV are shown in 
Figure 4.2.13. 
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Fig. 4.2.13: Variation of the 5a discovery potential for the pp — + H/A + X, H/A t^t^ process at the LHC in the m™"'^ 
(left) and no-mixing (right) benchmark scenarios for different values of fi. 

The change in the upper Umit of the LHC wedge region due to the variation of /i does not exceed 
A tan /3 ~ 8. As explained above, this is a consequence of cancellations of the leading A;, effects in 
the Higgs production and the Higgs decay. Besides the residual Af, corrections, a further variation of the 
bounds is caused by the decays of the heavy Higgs bosons into supersymmetric particles. For a given 
value of f-i, the rates of these decay modes ai^e strongly dependent on the particular values of the weak 
gaugino mass parameters M2 and Mi. In our analysis, we have taken M2 = 200 GeV, as established 
by the benchmai^k scenarios defined in Ref. [53], while Mi ~ 100 GeV. In general, the effects of the 
decays H/A — > Xi'Xj; X^xf only pl^Y ^ role for Ma ^ \n\ + Mi. Outside this range the cancellations 
of the Afo effects result in a very weak dependence of the rates on /i. The combination of the effects 
from supersymmetric radiative corrections and decay modes into supersymmetric particles gives rise to 
a rather complicated dependence of the discovery contour on fi, see Ref. [54] for more details. 

Discovery region for the process 

For this process we also refer to the analysis carried out by the CMS Collaboration [119, 126]. The 
corresponding analyses of the ATLAS Collaboration can be found in Refs. [118, 127, 128]. The results 
of the CMS Collaboration were given for an integrated luminosity of 30 fb~^ in the My^-tan (3 plane 
using the m™^^ scenario with p = —200 GeV. No A;, corrections were included in the gh tH^ 
production process [129] and the — > tu-j- decay [68]. 

In Figure 4.2.14 we investigate the impact of including the Af, corrections into the production 
and decay processes and of varying fi. In order to rescale the original result for the discovery reach 
in tan /3 we have first evaluated the tan /3 dependence of the production and decay processes. If no 
supersymmetric radiative corrections are included, for a fixed Ma value, the discovery potential can be 
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Fig. 4.2.14: Variation of the 5a discovery contours obtained from the tH^,H^ tUt channel in the m™"'' (left) and 
no-mixing (right) benchmark scenarios for different values of jj,. 

inferred by using a rate approximately proportional to 

tan2 /3cMS X BRiH^ ^ rz^^cMS • (4.2.47) 

Here tan PcMS is given by the edge of the area in the M^i-tan /3 plane in which the signal tUt is 

visible, as obtained in the CMS analysis. The BR(if^ — > rz^T-)cMS has been evaluated with HDECAY. 

The rescaled result for the discovery contour, including all relevant Af, corrections, is obtained by 
demanding that the contribution 

j^^^ X BR(F± ^ , (4.2.48) 

where Af, depends on tan /?, is numerically equal to the one of Equation 4.2.47. The quantities in 
Equation 4.2.48 have been evaluated with FeynHiggs. 

This procedure has been carried out in two benchmark scenarios for various values of fi. The 
results are shown in Figure 4.2. 14 for the m™*^^ scenario (left) and for the no-mixing scenario (right). 
As a consequence of the cancellations of the leading A^ effects in the Higgs production and the Higgs 
decay the change in the discovery contour due to the vaiiation of f^i does not exceed A tan /3 ~ 10(6) 
in the m™'^^ (no-mixing) scenario. Also in this case there is a variation of the contour caused by decays 
into supersymmetric particles that, as in the neutral Higgs boson case, are only relevant for small values 

of |;U|. 

4.3 Benchmark Scenarios 

The benchmark scenarios defined in Ref. [53], which were mainly designed for the search for the light 
CP-even Higgs boson h in the CP-conserving case, are also useful in the search for the heavy MSSM 
Higgs bosons H, A and H^. In order to take into account the dependence on fi, which as explained 
above is particularly pronounced for the bbcj), 4> —>■ bb channel, we suggest to extend the definition of the 
j^max ^j^j no-mixing scenarios given in Ref. [53] by several discrete values of fi. The scenarios defined 
in Ref. [53] read 



41 



no-mixmg: 



mt = 174.3 GeV, MsusY = 1000 GeV, /x = 200 GeV, Ma = 200 GeV, 

X^os ^ 2MsusY (FD calculation), Xf^ = \/6MsusY (Recalculation) 

Ab = At,mg = 0.8 MsusY • (4.3.49) 

mt = 174.3 GeV, MsusY = 2000 GeV, ^ = 200 GeV, Ma = 200 GeV, 

Xt = (FD/RG calculation) Ab = At, m~g = 0.8 MsusY • (4.3.50) 



constrained m^^^: 

mt = 174.3 GeV, MsusY = 1000 GeV, = 200 GeV, Ma = 200 GeV, 

X^^ = -2MsusY (FD calculation), = -VgMsusy (Recalculation), 

Ab = At, m-g = 0.8 MsusY • (4.3.51) 

The constrained-m™'^^ scenario differs from Equation 4.3.49 only by the reversed sign of Xt. While the 
positive sign of the product {n Ma) results in general in better agreement with the [g — 2)^ experimental 
results, the negative sign of the product {fi At) yields in general (assuming minimal flavor violation) 
better agreement with the BR(6 57) measurements. 

Motivated by the analysis in Section 4.2 we suggest to investigate the following values of 

^JL = ±200, ±500, ±1000 GeV , (4.3.52) 

allowing both an enhancement and a suppression of the bottom Yukawa coupling and taking into ac- 
count the limits from direct searches for charginos at LEP [130]. As discussed above, the results in the 
constrained-m™^^ scenario are expected to yield more robust bounds against the variation of n than in 
the other scenarios. It should be noted that the values /i = —500, —1000 GeV can lead to such a large 
enhancement of the bottom Yukawa coupling that a perturbative treatment is no longer possible in the 
region of very large values of tan j3. Some care is therefore necessary to assess up to which values of 
reliable results can be obtained, see e.g. the discussion of Figure 4.2. 12. 

The value of the top-quark mass in Ref. [53] was chosen according to the experimental central 
value at that time. We propose to substitute this value with the most up-to-date experimental central 
value for mt. 



4.4 Conclusions 

In this paper we have analyzed the impact of supersymmetric radiative coiTcctions on the current MSSM 
Higgs boson exclusion limits at the Tevatron and the prospective discovery reach at the LHC. In partic- 
ular, we have studied the variation of the exclusion and discovery contours obtained in different MSSM 
benchmark scenarios under changes of the higgsino mass parameter /i and the supersymmetry breaking 
parameters associated with the third generation squarks. These parameters determine the most important 
supersymmetric radiative corrections in the large tan [3 region that are associated with a change of the ef- 
fective Yukawa couplings of the bottom quarks to the Higgs fields (since the squarks are relatively heavy 
in the considered benchmark scenarios, other squark-loop effects are sub-dominant). These corrections 
had been ignored or only partially considered in some of the previous analyses of Higgs searches at 
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hadron colliders. We have shown that their inclusion leads to a significant modification of the discovery 
and exclusion regions. 

We have investigated the exclusion bounds obtained from the Tevatron searches for non SM-lrke 
Higgs bosons in different channels. For the bb(p, (p ^ bb channel {(j) = h, A) we find that the effects 
of the supersymmetric radiative corrections on the exclusion bounds in the M^i-tan /5 plane are quite 
dramatic. While in the m™^^ scenario the current data allow to rule out values of tan /3 > 50 (35) for 
Ma ~ 100 GeV if the higgsino mass parameter is chosen as /i = —200 GeV (—1000 GeV), hardly any 
bound on tan P can be set if positive values of fi are chosen. The shifts are smaller, but still important, 
for the no-mixing benchmark scenario. We have shown that the constrained-m™'^^ scenario yields results 
that are much more stable against variations of /i than the other benchmark scenarios. 

For the inclusive channel with t~^t~ final states, pp ^ (f) ^ t~^t~, compensations occur between 
large corrections to Higgs production and decay, so that the limits in the M^-tan /? plane obtained from 
this channel turn out to be less affected by varying /i than the ones from the associated production with 
bottom quarks. Nevertheless we have found that the exclusion limit is shifted by up to A tan /3 = 30 
as a consequence of choosing different input values for ^u. We have investigated the impact of including 
the dominant supersymmetric radiative corrections to the gluon fusion production process, which had 
previously been omitted. The inclusion of these corrections leads to a shift of up to A tan /? = 10 in the 
exclusion limit. Following our analysis, the CDF Collaboration has adopted the prescription outlined in 
this paper for incorporating the correction into the gg ^ (p production process. The Tevatron experiments 
are expected to collect further data at higher luminosities, up to 4-8 fb~^, in the next few years. This 
will extend the Tevatron MSSM Higgs boson discovery and exclusion reach in the M^-tan /3 plane to 
lower values of tan /?, decreasing the sensitivity of the obtained bounds to variations of the low energy 
supersymmetry mass parameters. 

For the LHC we have analyzed the channels pp — > H/A + X, H/A — > t~^t~ - and tH^,H^ 
TVr, which yield the best sensitivities in the search for heavy MSSM Higgs bosons. Accordingly, the 
discovery contours for these channels in the M^-tan (3 plane determine the boundary of the region where 
only the (SM-like) light CP-even Higgs boson can be detected at the LHC. Since the discovery contours 
for the LHC ai^e at smaller values of tan /? compared to those accessible via the cuiTcnt exclusion bounds 
at the Tevatron, the impact of the tan /3-enhanced supersymmetric corrections is less pronounced in this 
case. We have studied the effect of including the dominant supersymmetric corrections, which had been 
omitted in the analyses of the production processes at the LHC, and their variation with the relevant 
parameters. Possible decays of the heavy MSSM Higgs bosons into charginos and neutralinos have been 
taken into account. We have found that the prospective discovery contours at the LHC are shifted by up 
to A tan /? < 10. 

Based on our analysis of the sensitivities of the searches for MSSM Higgs bosons at the Tevatron 
and the LHC we have defined benchmark scenarios for the analysis of MSSM Higgs-boson searches at 
hadron colliders. They are based on a generalization of similar benchmark scenarios proposed for the 
searches for SM-like MSSM Higgs bosons at the Tevatron and the LHC. 
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5 Sensitivity of CDF's Higgs Boson Searches 
Contributed by: T. Junk for the CDF Higgs Group 

The search for the Standard Model (SM) Higgs boson is one of the central pieces of the cuiTcnt 
High Energy Physics program. The SU{2) x U{1) gauge model of electi'oweak interactions makes a 
number of predictions which have been experimentally verified to high precision, but its validity depends 
on the breaking of this symmetry to the C/(1)em symmetry group at low energies. Many differing 
proposals of the details of this symmetry breaking have been advanced, most of which predict one or 
more observable scalar bosons. If the minimal SM Higgs mechanism describes nature, then precision 
electroweak data [5] provide evidence that the scalar Higgs boson should be lighter than about 200 GeV, 
with a preferred value at around 115 GeV. Direct searches at LEP have excluded [131] a SM Higgs 
boson with a mass below 1 14.4 GeV. If there is a SM Higgs boson with a mass between 115 GeV and 
~ 200 GeV it is produced in pp collisions at the Tevatron, and, with enough data, it should be possible 
to exclude or discover such a particle. 

Data are being accumulated by the Tevatron experiments CDF and D0, whose runs are expected 
extend until 2009. CuiTcntly, more than 1 fb~^ of data have been recorded by each experiment, al- 
though the Higgs boson searches reported here are based on approximately 300 pb^^ of data. The exact 
luminosities used in the channels is listed in Table 5.0.4. 

With 300 pb~^ of data and the expected signal-to-background ratios in the channels, the SM Higgs 
boson hypothesis cannot be tested for any value of mn- Nonetheless, with additional data, and improve- 
ments to the analyses, sensitivity at the 95% CL level may be obtained for mn up to 180 GeV, assuming 
the design integrated luminosity of 8 fb~^ is collected with good quality by both detectors, according to 
a 1999 study [4]. An updated study [132] was conducted in 2003 to check the earlier projections with 
more realistic simulations and preliminary data samples which could be used to caUbrate some back- 
grounds. The later study did not consider searches for Higgs bosons with uih greater than 130 GeV, 
and also did not include the effects of systematic uncertainties on the amount of luminosity required to 
test for Higgs bosons. Each report includes calculations of the estimated amounts of luminosity required 
for a combination of all of CDF's channels and D0's channels to exclude at the 95% CL, assuming a 
Higgs boson is not present, as well as the luminosity requirements for a combined 3a evidence and 5a 
discovery. The luminosity thresholds ai^e shown in Figure 5.0.15 for the 1999 study and in Figure 5.0.16 
for the 2003 study. 

The CDF channels as they stand as of the Summer 2005 conferences are not as powerful as those 
assumed in the two sensitivity studies. The following sections provide a snapshot of the sensitivity of 
the CDF channels separately and combined, as of the October 2005 TeV4LHC workshop, with plans for 
improvement. 

5.1 Sensitivity by Channel 

The expected signal and background rates and shape distributions were collected from each of the chan- 
nel analysis teams and combined using the CLg technique [137, 138] to find the expected limits on the 
cross-section multiplied by the branching fractions. Candidate information was not included in the com- 
bination, so the observed limit of the combination is not computed. All of the observed limits in the 
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Fig. 5.0.15: SUSY/Higgs Working Group estimations of thie luminosity required for 95% exclusion, 3cr evidence, and 5cr 
discovery for the combined CDF+D0 search channels. (2000). 




Fig. 5.0.16: Higgs Sensitivity Working Group estimations of the luminosity required for 95% exclusion, 3cr evidence, and 5cr 
discovery for the combined CDF+D0 search channels, compared against the earlier SUSY/Higgs Working Group's calculation. 
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Table 5.0.4: Integrated luminosities by channel. 



Channel 


f Cdt (pb-^) 


Reference 




319 


[133] 


ZH uvbh 


289 


[134] 


gg^ H ^ W+W- 


360 


[135] 


W^H W^W+W- 


194 


[136] 



Table 5.1.5: Relative systematic uncertainties by channel. Errors from the same source are considered correlated, across 
channels, and between signal and background. The "uncorrelated" errors are uncorrelated across channels and between signal 
and background. 





Channel 




W^H - 


-* e^ubb 


ZH - 


-> ui/bb 


ZH - 


-+ e+e-bb 


gg-* H ^ W+W- 


VK±// - 


-* w^w+w- 


Source 


s[%] 


b [% ] 


s[%] 


b[%] 


s [% ] 


b [% ] 


s [% ] b [% ] 


s[%] 


b [% ] 


lumi 


6 


6 


6 


6 


6 


6 


6 6 


6 


6 


b-tag s.f. 


5 




64 


1.9 


15 


15 




0.37 




lepton ID 


5 








7 


7 








lepton trig 


0.6 








1 


1 




24 




PDF 


1 














1.5 




ISR 


3 














3.0 




FSR 


7 








7 






3.2 




JES 


3 




7.8 


3.5 












Jet model 


1.4 


















ui/bb trig 






3 


1.5 












uubb veto 






2 


2 












uncorrelated 




15 


2 


22.1 




9 


6 7 


3.7 


66 



channels are close to expectations, the observed limit of the combination is expected to be close to the 
expected combined limit. 

The W^H l^ubb Channel 

The results of the W^H — > l^ubb search are described in [133]. The reconstructed mass distribution 
in the single-tagged analysis is used in computing the expected limits, with each bin counted as an 
independent counting experiment. Systematic errors are taken on the background and signal rates, but 
the shapes are not varied. Each bin is assumed to have fully correlated systematic uncertainties with all 
other bins of the mass distribution. The systematic uncertainties are detailed in Table 5.1.5. Acceptances 
and signal distributions are linearly interpolated [139] between the supplied test points at which Monte 
Carlo samples are available. The observed and expected cross-section times branching ratio limits are 
shown at the 95% CL in Figure 5. 1. 17 as a function of mn- 

The ZH — > uuhb Channel 

The results of the ZH uvbb search are described in [134]. The reconstructed mass distribution was 
not provided for combination, but the numbers of events for the expected signal and background after 
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Fig. 5.1.17: The observed and expected 95% CL limits on the production cross-section times the Higgs decay branching ratio 
as a function of niH for the W^H I'^vbb channel. The limits are compared with the SM prediction. 
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Fig. 5.1.18: The observed and expected 95% CL limits on the production cross-section times the Higgs decay branching ratio 
as a function of mn for the ZH vvbh channel. The liimts are compared with the SM prediction. 



a mass window cut which moves with the Higgs boson mass under test are used. They are linearly 
inteipolated between the model points listed in [134]. The systematic uncertainties on the signal and 
background are detailed in Table 5.1.5. The observed and expected cross-section times branching ratio 
limit is shown at the 95% CL in Figure 5. 1 . 1 8 as a function of itih, and compared to the SM expectation. 



The ZH e+£-bb Channel 

The ZH i^£~bb channel is still in development and the analysis is still in its "blind" stage. Hence, 
data candidate information is not yet available. The current status is described in [140]. The selection 
starts with a very clean sample of Z — > i'^i" decays, identifying isolated leptons with mu close to 
mz, and two or three jets, at least one of which must be b-tagged. The systematic uncertainties on the 
signal and background are detailed in Table 5.1.5. The neural net has seventeen input variables described 
in [140]. The most powerful ones are the invariant mass of the two leading jets taken together, the event 
Ht (which is the scalar sum of all the Pt's of the observed particles), and the Et of the leading jet. The 
median expected limit on the cross-section times the branching ratio for this process is approximately 
2.2 pb for 300 pb~^ of data. This expected limit is lower than that for other channels mainly due to the 
very small background prediction. It must be compared, however, against a much smaller SM signal 
expectation. 



48 



NN Output 




ZH 
Zbb 



o.oa 




Fig. 5.1.19: The distribution of the neural net discriminant function for the ZH — ^ bb channel, shown separately for the 
signal and for the major backgrounds, Zbb and Z + 2 partons. The data in this channel are still blind. 

The gg ^ H ^ W+W' Channel 

The results of the gg ^ H ^ W^W~ search are described in [135]. The histograms of A(/)^£ are used 
as the discriminant variable input to the limit calculation - each bin is a separate counting experiment. 
The shapes are interpolated [139] between rriH points, as are the signal rates and background rates. 
The analysis uses mj;/-dependent cuts, and so the background rates depend on the itih under test. The 
systematic uncertainties on the signal and background are detailed in Table 5.1.5. 

The median expected 95% CL cross-section times branching ratio limit is shown in Figure 5.1.20 
as a function of niH compared to the SM expectation and to the computation of [135]. 

The W^H W^W+W- Channel 

The results of the W^H W^W+W~ search are described in [136]. It is a single counting 
experiment - there are no discriminant variables whose histograms have different s/h ratios to use. The 
acceptance is interpolated between the mu points listed in [136]. The systematic uncertainties on the 
signal and background are detailed in Table 5.1.5. For this calculation, the data statistical uncertainty 
on the residual conversion background is treated as independent of the other errors on the background 
and the errors add in quadrature instead of linearly as they do in [136]. Furthermore, the FSR systematic 
uncertainty is almost certainly truly uncorrelated with other channels' FSR uncertainty, but it has been 
treated as correlated. As is seen below, the entire systematic error treatment in this channel matters little 
to the sensitivity. 

The observed cross-section times branching ratio limit is shown at the 95% CL in Figure 5.1.21 as 
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Fig. 5.1.20: The observed and expected 95% CL limits on the production cross-section times the Higgs decay branching ratio 
as a function of mn for the gg ^ H ^ W^W~ channel. The limits are compared with the SM prediction, and also the 
prediction of a model with a heavy fourth generation of SM-like fermions. 



a function of itlh compared to the SM expectation and to the computation of [136]. 
5.2 Sensitivity of the SM Channels when Combined 

The observed 95% CL limits in all of CDF's SM Higgs channels are shown, compared with SM pre- 
dictions, and also compared with observed limits from D0's channels, in Figure 5.2.22. The different 
searches search for different processes which have different rates, and thus contribute differently to the 
combined sensitivity. It is somewhat easier to compare the channels' sensitivity to a SM Higgs when 
the ratio of the limit in each channel to the SM prediction is formed. This ratio is shown for the same 
collection of CDF and D0 channels in Figure 5.2.23. 

The CLs method is used on the collection of CDF's five SM Higgs boson search channels to 
compute the multiplicative scale factor S95 on the total signal which can just barely be expected to be 
excluded in a median experimental outcome. This procedure doesn't make much physical sense for scale 
factors exceeding unity, as there isn't a well-motivated physical model which scales all of the production 
mechanisms for SM Higgs bosons in the same way, but it provides a technical benchmark of how far we 
are from the SM in our sensitivity. The results of this combination are shown in Figure 5.2.24. It must 
be shown as a multiplicative factor of the SM prediction because of the different SM predictions used for 
each search channel. 



5.3 Necessary SM Channel Improvements 

The current channels as we have them are insufficient to test for the presence or absence of the Standard 
Model Higgs boson, even if the projected 8 fb~^ of data are collected. Improvements must be made to 
increase the acceptance, reduce the background, and to separate the selected events into disjoint subsets 
with different s/b ratios, and to combine them together. Furthermore, the results must be combined with 
D0. 
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Fig. 5.1.21: The observed 95% CL limit on the production cross-section times the Higgs decay branching ratio as a function 
of mn for the W^H W^W^W channel. The limits are compared with the SM prediction, and also the prediction of a 
"bosophilic" (also known as "fermiophobic") model. 
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Fig. 5.2.22: The observed 95% CL limits on tiie production cross section times the Higgs decay branching ratio for each of the 
five search channels, compared with D0's limits, and also compared with SM expectations. 
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Fig. 5.2.23: CDF and D0's observed 95% CL limits on the production cross section times the Higgs decay branching ratio, 
divided by the corresponding SM predictions, for each of the five search channels. 
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Fig. 5.2.24: The expected 95% CL limit on the multiplicative scale factor of SM Higgs boson production for CDF's five SM 
Higgs boson search channels combined, as a function of mn, assuming the absence of a Higgs boson. The yellow and green 
bands show the ±1(t and ±2a expectations, which fluctuate depending on the possible data which may be observed. 
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Table 5.3.6: Luminosity factors expected from analysis improvements, separated by channel. 



lllllji VJ V tllltllL 


vv ±1 ^ -c ly uu 


7,H ^ uDhh 


ZH ^ e~^P~hh 


Tn zLT RpsoliitioTl 


1.7 


1.7 


1.7 






1 s 


1 s 


Forward B-tags 


1.1 


1.1 


1.1 


Forward Leptons 


1.3 


1.0 


1.6 


Neural Nets 


1.75 


1.75 


1.0 


Track-Only Leptons 


1.4 


1.0 


1.6 


WH signal in ZH 


1.0 


2.7 


1.0 


Product of above 


8.9 


13.3 


7.2 


CDF+D0 Combination 


2.0 


2.0 


2.0 


All Combined 


17.8 


26.6 


14.4 



The Higgs Sensitivity Working Group report [132] lists changes which can be made to the analyses 
which can get us to the desired level of sensitivity. Much of this work has already been done to improve 
our resolutions, to increase our lepton acceptance to the forward region, and to develop neural nets. 
But the work has been done by a variety of different people separated in space, time, and institution. 
The work of many groups must be collected together in the analysis channels in order to achieve the 
sensitivity reported in [4, 132]. 

The factors on the expected amount of luminosity needed to get exclusion at the 95% CL, 3(7 
evidence and ha discovery can be computed for most of the improvements rather easily. For acceptance 
increases, the background ought to increase as the signal acceptance increases. In fact, it should increase 
faster, because as we expand our acceptance to forward regions of the detector or to include leptons of 
lower quality, a larger fraction of background is expected to creep in. For this estimation, the estimations 
are taken from the HSWG report's Sections 2.3 and 4.2 (for the Neural Net factor). A listing of improve- 
ments and their factors in luminosity is given in Table 5.3.6. It is assumed in the luminosity projections 
that the systematic uncertainties will scale inversely with the square root of the integrate luminosity. 
Furthermore, accounting of the shape uncertainties may make the systematic errors larger. 

The neural net factor of 1.75 is not uniformly applicable to all channels, as the ZH — > £~^£~bb 
channel estimations already take advantage of a neural net. The forward lepton acceptance improvement 
cannot strictly be multiplied by the track-only lepton factor since the forward tracking is not sufficient. 
Nonetheless, a naive product of the factors from the analysis improvements is approximatley 20. The 
analysis improvements will not be made all at once - work is ongoing to develop and chai^acterize the 
techniques. 

5.4 SM Sensitivity Projections 

Assuming that the acceptances of the channels are increased and neural nets or other advanced techniques 
are used to reduce the backgrounds, the projected reach of of the Tevatron SM Higgs search program is 
estimated. It is assumed that the systematic uncertainties scale inversely with the square root of the 
integrated luminosity, and that D0 contributes channels with the same sensitivity as CDF's and that they 
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Integrated Luminosity (fb'^) 

Fig. 5.4.25: The evolution of the expected significance of an excess in the data if a Standard Model Higgs boson is present with 
a mass of 115 GeV. The yellow (light) interior band shows the ±lcr distribution of the expected significance, and the green 
(darker) exterior band shows the ±2(t range around the expectation. CDF and D0 are combined, and the foreseen sensitivity 
improvements have been assumed. The integrated luminosity is per experiment. 

are combined together. Figure 5.4.25 shows how the significance of an excess of events is expected to 
develop, as a function of the integrated luminosity collected per experiment, assuming a SM Higgs boson 
is present with a mass mn = 115 GeV. The actual evolution of such an excess, if a signal is actually 
present, will be more of a random walk as data are collected, so the figure also includes the width of 
the expected distribution. Figure 5.4.26 shows the evolution of the probability of seeing a 2a, a 3a, or a 
5a excess in the combined data when searching for a SM Higgs boson of mass mn = 115 GeV, if it is 
truly present, as a function of the luminosity collected by each experiment. After collecting 8 fb~^ per 
experiment, it is 10% likely that a 5cr excess will be observed if mn is truly 115 GeV. 

5.5 The MSSM H/h/A t+t' Sensitivity 

CDF has published its search for H/h/A t~^t~ search, using 310 pb^i of Run 2 collision data [89]. 
Tau pairs are selected in which one tau decays leptonically, and the other decays semi-hadronically. Kine- 
matic selection requirements were designed to separate tau pairs from M^+jets and QCD backgrounds, 
in which jets are misidentified as taus. The dominant remaining background is Z —>^ t'^t" production. 
In order to separate H/h/A —>■ t^t" from this and other backgrounds, the invariant mass of the visible 
tau decay products is formed, shown in Figure 5.5.27. The reconstructed invariant mass of Higgs boson 
signal events peaks near the signal mass, with a width which grows rapidly with increasing Higgs boson 
signal mass. This is offset by the fact that the background is very small for large reconstructed masses. 
The observed and expected limits on the production cross section times the decay branching ratio to tau 
pairs is shown in Figure 5.5.28. 

This cross-section Umit can be interpreted in the MSSM; we choose to represent it as an exclusion 
in the (m^, tan P) plane in the inh — max and 7io — mixing MSSM benchmark scenarios [141]. This 
interpretation benefits from the fact that for larte tan /?, two Higgs bosons (either h and A, or H and A), 
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Fig. 5.4.26: The fraction of experiments expected to make an observation of a 115 GeV SM Higgs boson if it is truly there, 
as a function of the integrated luminosity. CDF and D0 are combined, and the foreseen sensitivity improvements have been 
assumed. Separate curves are shown for the fraction of experiments observing a > 2(t excess in the data, a > 3(t excess, or a 
> 5(T excess. 
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Fig. 5.5.27; The invariant mass of the reconstructed tau decay products in the MSSM H t t~ search. The data (points) are 
compared to a sum of background predictions. A Higgs boson signal of mass niA ~ 140 GeV, with a production cross section 
just at the exclusion threshold, is shown. 
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Fig. 5.5.28: The 95% CL limit on the production cross section times the decay branching ratio for Higgs bosons decaying to 
tau pairs, using 310 pb~^ of CDF data, as a function of the Higgs boson mass. 



are nearly degenerate in mass and contribute rougly equally to the expected signal. CDF's observed 95% 
CL limits are shown in Figure 5.5.29, along with projected CDF+D0 combined sensitivity contours for 2, 
4, and 8 fb~^ of data collected by both CDF and D0. The large improvement in sensitivity at larger Higgs 
boson masses comes from the fact that the background rate is very low for large invariant-mass tau pairs. 
For a search with a large background rate, the expected signal limit is roughly inversely proportional to 
the square root of the integrated luminosity, while for searches with very small backgrounds, the expected 
limit is roughly inversely proportional to the integrated luminosity. 



58 




m^(GeV/c2) 

Fig. 5.5.29: The observed 95% CL limits in the tau channel in the (mA,tan/3) plane, for the m£f-max MSSM benchmark 
scenerio and also the no-mixing benchmark scenario, using 310 pb~^ of CDF data. Projections are shown for the expected 
combined CDF+D0 exclusion reach for 2, 4, and 8 fb^^ per experiment. 
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6 Two-Loop EW Corrections to Higgs Production 
Contributed by: U. Aglietti, R. Bonciani, G. Degrassi, A. Vicini 

We study the impact of the two-loop electroweak coiTections on the production of a Higgs boson 
via gluon-fusion in proton-proton collisions at LHC energies. We discuss the prescritpion to include 
the corrections to the hard scattering matrix element in the calculation of the hadronic cross-section 
a{p + p H + X). Under the hypothesis of factorization of the electroweak corrections with respect to 
the dominant soft and coUinear QCD radiation, we observe an increase of the total cross-section from 4 
to 8%, for nih < 160 GeV. This increase is comparable with the present QCD uncertainties originating 
from hard scattering matrix elements. 

6.1 Introduction 

We study the impact of the two-loop electroweak corrections on the production of a Higgs boson via 
gluon-fusion in proton-proton collisions at LHC energies. We discuss the prescritpion to include the 
corrections to the hard scattering matrix element in the calculation of the hadronic cross-section a-{p + 
p ^ H + X). Under the hypothesis of factorization of the electroweak corrections with respect to the 
dominant soft and coUinear QCD radiation, we observe an increase of the total cross-section from 4 to 
8%, for rrih < 160 GeV. This increase is comparable with the present QCD uncertainties originating 
from hard scattering matrix elements. 

The Higgs boson is one of the missing ingredients of the Standard Model and its discovery rep- 
resents one of the most important physics goals of the LHC. This goal will be achieved only if we can 
predict with high accuracy all the production cross sections of this particle and if we understand in detail 
the different decay channels and the relative backgrounds. 

At the LHC, the gluon-fusion is the dominant production mode over the entire range of interest- 
ing values of the mass of the Higgs particle (100 GeV < uih ^ 1 TeV). In particular, in the range 
100 GeV < rriH ^ 2mt this production mode is larger by almost one order of magnitude with respect to 
the next important mechanism, the vector boson fusion. It is, therefore, very important to have a precise 
prediction of its cross section and a reliable estimate on the remaining theoretical accuracy. 

The total cross section for the Higgs boson production by gluon fusion in the LO approximation 
was calculated in the late seventies [142]. It is an ©(a^a) calculation, since the Higgs couples to 
the gluons only via an heavy-quark loop (the most important contribution is the one due to the loop 
of top). For what concerns the higher orders, the calculation of the NLO QCD corrections have been 
done in the infinite rrit (mass of the top) limit in [143], and, with the full quark mass dependence, in 
[20]. Besides of the fact that the infinite mt approximation should be valid in the Higgs mass range 
rriH ^ 300 GeV, it has been noticed [144] that this approximation works also for values of mn beyond 
the top threshold, and up to masses of 0(1 TeV). The total effect of the NLO QCD corrections is the 
increase of the LO cross section by a factor 1.5-1.7, giving a residual renormalization/factorization scale 
dependence of about 30%. The unexpected size of the NLO QCD radiative corrections made in such a 
way that the electroweak corrections, evaluated in the infinite mt limit in [145, 146] and turning out to 
amount of less that 1%, were totally neglected. The attention was driven by the evaluation of the NNLO 
QCD corrections, carried out in the infinite mt limit by several groups [147, 148, 149, 150, 151, 152]. 
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Fig. 6.2.30: Lowest order (a) and generic NLO-EW (b), (c), (d) Feynman diagrams. The solid lines are fermions. The wavy 
lines are gauge bosons [V — W,Z). 

The calculation shows a good convergence of the perturbative series: while the NNLO corrections are 
sizable, they aie, nevertheless, smaller that the NLO ones. Moreover, the NNLO corrections improve 
the stability agaist renormalization/factorization scale variations. The effect due to the resummation of 
soft-gluon radiation has been included in [153], and the remaining theoretical uncertainty, due to higher- 
order QCD corrections, has been estimated to be smaller than 10%. Finally, several efforts were also 
devoted to the calculation of QCD radiative corrections to less inclusive quantities, such as the rapidity 
distribution, recently evaluated at the NNLO [154, 155], or the transverse momentum (qt) distribution 
[156, 157, 158, 159, 160, 161, 162, 163, 164, 165, 166, 167, 168, 169, 170, 171, 172, 173, 174, 175, 
176, 177, 178], which, in [179, 180], is evaluated using the fixed-order perturbative results up to NLO in 
QCD and the resummation up to the NNLL. 

Motivated by this accurate scenario, the NLO electroweak corrections to the gluon fusion were 
again taken into account recently. In [181, 182] the contribution to the partonic cross section due to the 
light fermions were calculated. It turned out that they are sizeable. In particular, in the intermediate 
Higgs mass range, from 1 14 GeV up the the 2niw threshold, these corrections increase the LO partonic 
cross section by an amount of 4-9%. For larger values of the mass of the Higgs, uih > 2mw, they 
change sign and reduce the LO cross section; however, in this region the light-fermion corrections are 
quite small, reaching at most a -2%. In [183], also the remaining electroweak corrections due to the top 
quark were calculated as a Taylor expansion in m|^/(4m^). They are valid for inn ^ 2mw, range in 
which they have opposite sign with respect to the light-fermion connections. However, the corrections 
due to the top quark are smaller in size, reaching at most a 15% of the light-quark ones. 

The impact of the NLO electroweak corrections on the hadronic cross section has not been dis- 
cussed yet. We present here the effect of their inclusion in the calculation at the hadronic level. 

6.2 Inclusion of the Two-Loop Electroweak Corrections 

The partonic gluon fusion process occurs, in lowest order, via one-loop diagrams, as the one depicted 
in Fig. 6.2.30 (a); in the loop run only the top and the bottom quarks, because of the Yukawa suppres- 
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sion of the lighter quarks. The NLO-EW connections ai^e schematically represented by the diagrams in 
Figs. 6.2.30 (b), (c) and (d). In particular, in Figs. 6.2.30 (b) and (c) the WWH/ZZH couplings avoid 
the Yukawa suppression, and, therefore, in these diagrams the fermionic line represents all the possible 
flavours: light flavours, evaluated in [181, 182], and top quark, evaluated in [183]. In Fig. 6.2.30 (d), 
instead, the fermionic line can represent only the top quark [183]. 

At the hadronic level, we consider the Higgs boson production at the LHC, and therefore in proton- 
proton collisions. The hadronic cross section can be written as: 



where the partonic processes initiated by partons (a, b) are convoluted with the corresponding parton 
densities fi^p{x, M^), (i = a, h), evaluated at a scale M. The effect of the higher order QCD and EW 
corrections is described by the two functions xQCD-oniy 6ew> obtained by factoring the lowest 
order cross section uq. 

In the partonic cross section, QCD and EW corrections have been factorized. This ansatz is valid 
up to subleading higher order corrections which start at the 3-loop level (i.e. 0{aas) with respect to the 
lowest order). The factorization of the QCD initial state collinear divergences holds for the hard process 
described by the electroweak NLO corrections, following from general arguments of the factorization 
theorems and from the universal nature of the initial state collinear radiation. In fact, the whole set 
of EW corrections is characterized by a scale M^, much harder than the one typical of the leading 
collinear emission. In addition, in the limit of light Higgs, the EW corrections can be expanded as a 
Taylor series in powers of Jn^/M^ and the EW connections vertex becomes effectively pointlike. In this 
regime the factorization of the QCD collinear divergences becomes rigorous. For heavier Higgs masses, 
the factorization should still be valid only in leading order, due to the modifications induced by the EW 
form factor. 

6.3 Numerical Results 

The hadronic proton-proton cross section has been calculated at LHC energy, in NNLO-QCD accuracy, 
i.e. setting 6ew = 0, using the MRST2002 NNLO parton distribution functions [184]. The theoretical 
uncertainty due to the choice of the renormalization scale /i and of the factorization scale M has been 
canonically estimated by setting M = fi equal to m/j/2, mh,2mh respectively. The predictions, shown in 
Fig. 6.3.31 (dotted lines), vary by approximately ±8% with respect to the central value. This uncertainty 
is further reduced when including the effect of the resummation of all the initial state soft gluon radiation 



The two-loop electroweak conections have been added according to Eq. (6.2.53) and setting 
M = fj, = nih. The light fermion conections can be evaluated for any choice of nih, whereas the 
top quark contribution has been computed by means of a Taylor expansion and is limited to the region 
mh < 160GeV. The hadronic cross section increases from 4 to 8%, for nih < 160 GeV. As we can 





(6.2.53) 



[153]. 
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Fig. 6.3.31: The cross section ap+p^n+x, in pb, is plotted as a function of the mass of the Higgs boson, between 114 
GeV and 300 GeV. The dotted lines describe the band of NNLO-QCD uncertainty, for three values of the QCD factoriza- 
tion/renormalization scale jj, = mh/2,mH,2mh- The solid line is the NNLO-QCD (jj. — m,h) with the two-loop EW 
corrections, according to Eq. (6.2.53). The two-loop EW corrections include also the top-quark effect, for mn < 155GeV, but 
only the light quarks contribution for larger values of ruh ■ 
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ruh (GeV) 




nih (GeV) 




rrih (GeV) 


5ew 


mh (GeV) 


5ew 


114 


1.048 


136 


1.062 


158 


1.077 


180 


1.020 


116 


1.049 


138 


1.063 


160 


1.069 


182 


1.010 


118 


1.050 


140 


1.065 


162 


1.063 


184 


1.010 


120 


1.051 


142 


1.066 


164 


1.049 


186 


1.002 


122 


1.053 


144 


1.068 


166 


1.041 


188 


0.997 


124 


1.054 


146 


1.069 


168 


1.035 


190 


0.994 


126 


1.055 


148 


1.071 


170 


1.031 


192 


0.991 


128 


1.056 


150 


1.073 


172 


1.028 


194 


0.989 


130 


1.058 


152 


1.074 


174 


1.026 


196 


0.987 


132 


1.059 


154 


1.076 


176 


1.024 


198 


0.986 


134 


1.060 


156 


1.077 


178 


1.022 


200 


0.985 



Table 6.3.7: Rescaling factor Sew as a function of the Higgs boson mass. 



observe in Fig. 6.3.31, the effect of the electroweak corrections is an increase of the cross section by an 
amount which is of the same order of magnitude of the NNLO-QCD theoretical uncertainty, and possibly 
larger than the uncertainty estimated after the resummation of soft gluon radiation. The main source of 
uncertainty on the hadronic cross section remains in the accurate determination of the parton distribution 
functions of the proton. 

The effect of the NLO-EW corrections is of great interest, because it enhances the most important 
Higgs production mechanism and, in turn, affects the absolute number of events of all the Higgs decay 
modes. 

Following Eq. (6.2.53), the NLO-EW corrections can be implemented as a simple rescaling of the 
QCD hadronic cross section. This multiplicative factor is presented in Table 6.3 as a function of and 
can be fitted, in the range 114 GeV < niH ^155 GeV, by the following simple formula: 

^Ewimn) = 1.00961 + 6.9904 • 10"^ nih + 2.31508 • 10"^ m| . (6.3.54) 

The computation of the NLO-EW corrections to the gluon fusion process has been described in detail in 
[181, 182, 183]. The analytical expression of the probability amplitude has been expressed in terms of 
Generalized Harmonic PolyLogarithms (GHPL) [185] and has been implemented in a FORTRAN rou- 
tine^2_ The GHPL can be evaluated numerically in several different ways: by direct numerical integration 
of the basic functions, by power expansions or by solving the associated differential equations. We have 
checked that these fully independent approaches agree. 

6.4 Conclusions 

In conclusion, the calculation of the QCD corrections to the production of a Higgs boson via gluon-fusion 
has reached a very high level of accuracy; it is now mandatory the inclusion of the two-loop electroweak 
corrections, whose typical size for rrih < 160 GeV is larger than 5%, comparable or larger than the 

'^available upon request from the authors 
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QCD uncertainty. The main source of uncertainty on the hadronic cross section remains in the accurate 
determination of the parton distribution functions of the proton. 
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7 Higgs Resummation 
Contributed by: B.J. Field 



Resummation when combined with fixed-order calculations yield the most accurate theoretical 
definition for differential cross-sections for use by the experimental community to coiTcctly determine 
and unambiguously reconstruct the mass of any Higgs bosons that may exist in nature. There has been 
excellent progress in recent years at obtaining higher-order fixed-order results for all the relevant Higgs 
production processes. It is therefore necessary to follow these calculations with resummed calculations 
to complete our theoretical understanding ahead of future collider data. 

Resummation of processes involving both scalar and pseudoscalar Higgs bosons have been studied 
extensively[186, 187, 188, 189, 190, 144, 191, 192, 193, 194, 195, 196, 197, 23, 198, 199, 200, 201, 202]. 
In order to understand the small-p^- behavior of a differential cross-section, we need to employ the 
resummation formalism[203, 204, 205]. Resummation will smooth out any numeric instabilities and 
divergencies that occur in fixed order calculations. 

There are several ways to approach resummation of both total cross-sections and differential quan- 
tities. Here we will discuss only the differential quantities as they are more relevant for experimental 
concerns. The formalism of resummation is different depending on what kind of observable one is inter- 
ested in studying. Although much of the work of resummation has focused on inclusive production, it is 
possible to study exclusive processes[206]. 

For inclusive processes where typically only one particle is produced at lowest order, the standard 
formalism allows a resummed differential cross-section to be written as an inverse Bessel transformation 
in terms of an expansion in of both universal and process-dependent terms. For instance. 
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= ^ / dxi j dx2 / db-Jo{bpt) 
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and the Higgs mass is M| = Q"^, dcj) is any unintegrated phase space of the system under consideration, 
and (T^g '^^ is the lowest order cross-section with a cc initial state which is therefore defined at pt = 0. 
Since this is an inclusive process, the transverse momentum and rapidity in the differential cross-section 
{pt , y) are that of the Higgs boson produced. The integration variable b is the impact parameter, S is the 
center-of-momentum energy of the hadronic system, and the /j/^ are the parton distribution functions 
for a parton i in hadron A. The constant bo is written in terms of the Euler-Mascheroni constant = 
0.57721 . . . as 6o = 2e~'^^ to simplify the coming coefficient expressions. The coefficients Cab are 
process dependent and can be written as power series to be described below. jQ{bpt) is the first order 
Bessel function. The Sudakov form factor Sc, which makes the integration over the Bessel function 
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convergent, can be written as, 

Sc{Q,b) = expj- 



III? ' 

The coefficient functions A^., B^, and Cab can be written as power series in as as 



(7.0.55) 
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(7.0.56) 



(7.0.57) 



The A^c \ A^c \ and i?c^^ coefficients have been shown to be universal. 

Although this particular method of handling resummed differential cross-sections is useful, it ig- 
nores several of the more interesting channels where a Higgs boson would be produced in association 
with another particle, and therefore has nothing to say about the other particles in the process. There has 
been interest in a Higgs bosons produced in association with weak vector bosons, light and heavy quarks, 
and combinations of these particles as can easily be seen in the cross-sections that have been compiled 
for this workshop. In particular, there is continuing theoretical [4 1 , 207, 38, 37, 208, 39, 209, 40] and 
experimental[88] interest in a Higgs produced in association with bottom quarks, particularly in super- 
symmetric models where bottom quarks can play a role equal to or greater than top quarks in Higgs 
processes. 

To understand resummed processes for a more general particle configuration can be calculated 
with a different mechanism[206] which can be used for inclusive, exclusive, and the resummation of 
pair-invariant-mass quantities. This formalism is setup as a unified approach that allows one to perform 
the resummation based on the color flow of the chosen process as well as several previously calculated 
quantities. 

If we define a generic plus-distribution related object in terms of the variable that becomes soft 
(54) or zero (1 — z) on threshold as 



\n\si/M 

Si 



or. 



1 



(7.0.58) 



then at next-to-leading-log (NLL) accuracy we can define a total or differential cross-section, where we 
generically call the threshold variable Xth, as 



vr 



da = da^ — l C32?i(xth) + C2Po(2;th) + ci6{xi<ty] 



(7.0.59) 



where da can be either a total or differential quantity of interest and da^ is the total or differential 
Bom cross-section. There exists a similar equation of the NNLL corrections that will not be reproduced 
here but can be found in the literature [206]. The coefficients Cj can be calculated in terms of the color 
Casimir invariants of the partons involved in the process, the soft anomalous dimension matrix, and the 
kinematic invariants of the process. At higher orders, the expressions become more complicated, but 
straightforward to calculate in a unified manner for several different processes. 
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The results of the Higgs calculations in the literature tells us several qualitative facts about re- 
summed Higgs processes. First, the resummed quantities are finite at small-p^ and removed the fixed- 
order divergencies. We also find smaller scale uncertainties at higher-orders as expected. 

The two primary inclusive Higgs processes that have been studied are gg ^ ^ and bh — > <1>, where 
$ is generically any Higgs boson of interest. Some sample resummation calculation for inclusive Higgs 
from Refs [199, 200] are shown in Figure 7.0.32. 

Aside from the usual observations for these differential cross-sections, we can see that the total 
uncertainty in the magnitude of the cross-section at the peak, from parton distribution functions and a 
scale variation by a factor of two, is approximately 35%. This level of theoretical uncertainty is on par 
with similar fixed order calculations. 
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Fig. 7.0.32: Figure 7.0.32a shows tlie transverse momentum spectrum for a scalar Higgs boson produced via gluon-fusion at 
the LHC for \y\ < 2.5. The resummed curve is the NLL resummation, and the perturbative curve is the NLO fixed order 
calculation. The NLO fixed order calculation diverges in the negative direction at small pt. This piece of the differential 
cross-section is not shown for clarity. These two curves cross at approximately pt = 100 GeV/c and stay very close thereafter. 
Figure 7.0.32b shows the errors associated with the CTEQ6.1M PDF set for the inclusive resummation bb H. The variation 
is approximately 8 — 12%. 

Beyond the inclusive processes, there has been an excellent experimental use of tagged bottom jets 
to constrain the (M^, tan (3) plane in supersymmetric models[88], where Ma is the pseudoscalar mass 
and tan (3 is the ratio of the vacuum expectation values of the up- and down-quark sectors. However, the 
fixed-order results for the Higgs plus bottom quark process show a numeric instability in the region of 
the cut on the bottom quark[40] when calculated in the five flavor number scheme (5FNS). 

Using the unified approach for the resummation of exclusive quantities, the instabilities that oc- 
curred in the fixed-order calculation can be smoothed out in the region around the cut in the bottom 
quark transverse momentum in a fixed-order calculation become well behaved. However, it is important 
to understand that the exclusive resummation formalism is built on 2 ^ 2 kinematics and there is no 
way to add a cut in the transverse momentum of one of the outgoing particles without setting the differ- 
ential cross-section below that point to zero. However, one could use the smooth resummed calculation 
to further consti'ain the {Ma, tan 15) limits as the transverse momentum cut on the bottom jet is further 
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pushed down into the region where the fixed order calculation becomes increasingly unreliable. This 
investigation is currently underway. 
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8 Hadronic Higgs Production with Heavy Quarks at the Tevatron and the LHC 
Contributed by: S. Dawson, C.B. Jackson, L. Reina, D. Wackeroth 

A light Higgs boson is preferred by precision fits of the Standard Model (SM) and also theoret- 
ically required by the Minimal Supersymmetric extension of the Standai^d Model (MSSM). The pro- 
duction of a Higgs boson in association with a heavy quark and antiquark pair, both tt and bb, at the 
Tevatron and the Large Hadron Collider (LHC) will be sensitive to the Higgs-fermion couplings and can 
help discriminate between models. 

The associated production of a Higgs boson with a pair of tt quarks has a distinctive signature and 
can give a direct measurement of the top quark Yukawa coupling. This process is probably not observable 
at the Tevatron, but will be a discovery channel at the LHC for Mh < 130 GeV. The associated production 
of a Higgs boson with a pair of bb quarks has a small cross section in the Standard Model, and can be 
used to test the hypothesis of enhanced bottom quark Yukawa couplings in the MSSM with large values 
of tan (3. Both the Tevatron and the LHC will be able to search for enhanced bbh production, looking 
for a final state containing no bottom quarks (inclusive production), one bottom quark (semi-inclusive 
production) or two bottom quarks (exclusive production). 

The rates for tih production at the Tevatron and the LHC have been calculated at NLO QCD 
several years ago[43, 42, 210, 211, 44, 212]. The theoretical predictions for bbh production at hadron 
colliders involve several subtle issues, and depend on the number of bottom quarks identified in the final 
state. In the case of no or only one tagged bottom quark there are two approaches available for calculating 
the cross sections for bbh production, called the four flavor number schemes (4FNS)[37, 39] and five 
flavor number scheme (5FNS)[207]. The main difference between these two approaches is that the 4FNS 
is a fixed-order calculation of QCD corrections to the gg and gg-induced bbh production processes, while 
in the 5FNS the leading processes arise from bg (bg) and bb initial states and large coUinear logarithms 
are resummed using a pertubatively defined bottom quark Parton Distribution Function (PDF). Very good 
agreement is found for the NLO QCD corrected cross sections for bb Higgs associated production when 
the two schemes are compared[40, 208]. 

In the following sections, we present numerical results at NLO QCD for tih and bbh production at 
the Tevatron and the LHC. If not stated otherwise, numerical results have been obtained in the 4FNS. We 
emphasize theoretical uncertainties from scale and PDF uncertainties and also present differential cross 
sections at NLO for bbh production in the case when two 5 quarks are tagged. 

8.1 Results for tth Production 

The observation of a tth final state will allow for the measurement of the tth Yukawa coupling. If 
Mh < 130 GeV, pp — > tih is an important discovery channel for a SM-like Higgs boson at the LHC {^/s = 
14TeV) [213,214]. Given the statistics expected at the LHC, tt/i, with /i 66, r+r~, T^+VF", 77 
will be instrumental for the determination of the couplings of the Higgs boson. Precisions of the order 
of 10-15% on the measurement of the top quark Yukawa coupling can be obtained with integrated lumi- 
nosities of 100 fb~^ per detector[215, 216, 217, 218]. 

The impact of NLO QCD corrections on the total cross section for pp, pp tih production in the 
Standard Model is illustrated in Fig. 8.1.33[212, 44, 43, 210] and Fig. 8.1.34[212, 210]. The dependence 
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Fig. 8.1.34: Dependence of glo.n Lo(pp,pp tth) on the renormalization/factorization scale ^, at \fs = 2 TeV (l.h.s.)and 
^/s = 14 TeV (r.h.s.), for Mn = 120 GeV. 



of the total cross sections on the renormaUzation and factorization scales is strongly reduced at NLO as 
shown in Fig. 8.1.34. The numerical results at NLO are obtained using CTEQ4M (Fig. 8.1.34 (l.h.s.)), 
CTEQ5M (Fig. 8.1.34 (r.h.s.)), and CTEQ6M (Fig. 8.1.33) parton distribution functions. The NLO 
cross section is evaluated using the 2-loop evolution of as{iJi) with a^^^{Mz) = 0.116 (Fig. 8.1.34 
(l.h.s.)) and a^^°{Mz) = 0.118 (Fig. 8.1.34 (r.h.s.)) and Fig. 8.1.33), and nit = 174 GeV. The 
renormalization/factorization scale dependence, uncertainty on the PDFs, and the eiTor on the top quark 
pole mass, mt, are estimated to give a 15-20% uncertainty. 

8.2 Results for hhh Production 

The hbh production processes are only relevant discovery modes in the MSSM with large tan /?. To 
a good approximation, the predictions for the MSSM rates can easily be derived from the Standard 
Model results by rescaling the Yukawa couplings[40]. The dominant MSSM radiative correction to hhh 
production can be taken into account by including the MSSM connections to the hbh vertex only, i.e. 
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by replacing the tree level Yukawa couplings by the radiative coiTccted ones. We follow the treatment 
of the program FeynHiggs [219, 112] and take into account the leading, tan/3 enhanced, radiative 
corrections that are generated by gluino-sbottom and chargino-stop loops. For large tan (3, the bottom 
quark Yukawa coupling is enhanced and the top quark Yukawa coupling coupling is strongly suppressed, 
resulting in a MSSM bbh cross section that is about three orders of magnitude larger than the Standard 
Model cross section. For the Tevatron, we calculate the production rates for the lightest MSSM Higgs 
boson, h^, while for the LHC we consider the rate for the heavier neutral Higgs boson, H^}^ 

In the numerical evaluation of cross sections for the exclusive and semi-inclusive channels (bbh 
and bh + bh production), it is required that the final state bottom quarks have px > 20 GeV and pseudora- 
pidity I r/ 1 < 2.0 for the Tevatron and \ r]\< 2.5 for the LHC. In the NLO real gluon emission contributions, 
the final state gluon and bottom quarks are considered as separate particles only if their separation in the 
pseudorapidity-azimuthal angle plane, AR = y/ (Ar/)^ + (A(^)2, is larger than 0.4. For smaller values of 
AR, the four momentum vectors of the two particles aie combined into an effective bottom/anti-bottom 
quark momentum four-vector. 

If not stated otherwise, the numerical results at NLO are obtained using CTEQ6M PDFs, the 2- 
loop evolution of as{fi) with a^^'^ (Mz) = 0.118, and the MS renormalization scheme for the bottom 
quark mass and Yukawa coupling with 2-loop renormalization group improved MS masses. The bottom 
quark pole mass is chosen to be mi, = 4.62 GeV. 

8.3 Total Cross Sections for bbh Production 

We present total cross section results at NLO in the 4FNS in Fig. 8.3.35 for associated bb Higgs produc- 
tion in the MSSM with tan /3 = 40. The bands represent the theoretical uncertainty due to the residual 
scale dependence. They have been obtained by varying the renormalization (/i^) and factorization (fif) 
scales independently from ^o/4 to fio, where ^o = "ife + Mh/2. 

If the outgoing bottom quarks cannot be observed then the dominant MSSM Higgs production 
process at large tan/3 is gg — > {bb)h (the curve labelled '0 b'). The inclusive cross section is experi- 
mentally relevant only if the Higgs boson can be detected above the background without tagging bottom 
quarks. At the LHC, this process can be identified at large tan /3 by the decays to fi'^fi" and r+r^ for 
the heavy Higgs bosons, and A^, of the MSSM. At the Tevatron this process, with — > r+r", has 
been used to search for the neutral MSSM Higgs boson. If a single bottom quark is tagged then the final 
state is bh or bh (the curve labelled '1 b'). Although requiring a b quark in the final state significantly 
reduces the rate, it also reduces the background. A recent Tevatron study [88] used the search for neutral 
MSSM Higgs bosons in events with three bottom quaifc in the final state (bh^ + bh^ production with 

bb) to impose limits on the tan /3 and M^o parameter space. 

Finally, we show the fully exclusive cross sections for bbh production, where both the outgoing 
b and b quarks are identified (the curve labelled '2 b'). The exclusive measurement corresponds to the 
smallest cross section, but it also has a significantly reduced background. Moreover, both the exclusive 
and semi-inclusive bbh production modes are the only ones that can unambiguously measure the bottom 

"We assume Msusy = 1 TeV, M§ = 1 TeV, At = At = 2 TeV (h°). At = At = 25 GeV (H"), fi = Ah = 200 GeV 
{h°), and = M2 = 1 TeV (H'^). For M,,o = 120 GeV, the bbh° coupling is enhanced by a factor of 33 relative to the 
SM coupling, while for Mjjo between 200 and 800 GeV, the bbH'^ coupling is enhanced by a factor of 27 relative to the SM 
coupling. 
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Fig. 8.3.35: Total cross sections for pp,pp bbh in the MSSM in the 4FNS at NLO for the Tevatron and the LHC in the 
MSSM with tan P — 40 and with 0,1 or 2 6 quarks identified. The Tevatron (LHC) plot is for the lightest (heaviest) neutral 
Higgs boson, h" (H°). The error bands have been obtained by varying the renormalization and factorization scales as described 
in the text. 
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Fig. 8.4.36: The relative corrections daNLo /doLo — 1 for the pr distribution of the bottom or anti-bottom quark with the largest 
Pt (pT,max) (left) and of the SM Higgs boson (p^) (right) to bbh production in the SM at the Tevatron (with ^/s = 2 TeV and 
/X = 2/io) and the LHC (with = 14 TeV and fi = 4/io). 



quark Yukawa coupling. 

8.4 Differential Cross Sections for bbh Production 

In assessing the impact of the NLO corrections it is particularly interesting to study the kinematic distri- 
butions. In Figs. 8.4.36 and 8.4.37 we illustrate the impact of NLO QCD corrections on the transverse 
momentum and pseudorapidity distribution of the SM Higgs boson and the bottom quark by showing the 
relative correction, daNio/daio — 1 (in percent) for the exclusive case {bbh where both b quarks are 
observed). For the renormalization/factorization scale we choose fi = 2^q at the Tevatron and fi = 4/io 
at the LHC, with fiQ = rrih + Mh/2, and use the CTEQ5 set of PDFs. As can be seen, the NLO QCD 
corrections can considerably affect the shape of kinematic distributions, and their effect cannot be ob- 
tained from simply rescaling the LO distributions with a K-factor of cjjvio/o'Lo = 1-38 it 0.02 (Tevatron, 
fi = 2fio) and (T^rio/c^Lo = l-H ± 0.03 (LHC, /i = 4/io).''^ 

'^The kinematic distributions have been calculated within the Standard Model and using the on-shell scheme for the definition 
of the b quark mass, but we see a similar behavior when using the A/5' bottom quark Yukawa coupling. 
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Fig. 8.4.37: The relative corrections (Ignlo /doLo — 1 for the r] distribution of the bottom quark r)b (left) and of the SM Higgs 
boson {r]h) (right) to hhh production in the SM at the Tevatron (with y/s = 2 TeV and ^ — 2/io) and the LHC (with \/s = 14 TeV 
and /i = 4/io)- 
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Fig. 8.5.38; Normalized cross sections for Higgs production with one 6 jet at the Tevatron (l.h.s) and the LHC (r.h.s) showing 
the uncertainty from PDFs for both the gg (4FNS) and bg (5FNS) initial states. 



8.5 PDF and Renormalization/Factorization Scale Uncertainties 

A major source of theoretical uncertainty for cross section predictions comes from the PDFs. We study 
the uncertainties of semi-inclusive bh production rates from the uncertainties in the PDFs using the 
CTEQ PDF sets[32]. First, the central value cross section do is calculated using the global minimum 
PDF (i.e. CTEQ6M). The calculation of the cross section is then performed with the additional 40 sets 
of PDFs to produce 40 different predictions, cjj. For each of these, the deviation from the central value 
is calculated to be Aaf = |crj — cjo| when cJi^cJo. Finally, to obtain the uncertainties due to the PDFs 

the deviations are summed quadratically as Acj^ = \Jj2i ^^f"^ ^''^'^ the cross section including the 
theoretical uncertainties arising from the PDFs is quoted as \ ^^^t • 

In Fig. 8.5.38, we plot the normalized total SM NLO cross sections for semi-inclusive hh produc- 
tion, calculated in the 5FNS {hg hh) as implemented in MCFM [220] and in the 4FNS {gg b{b)h), 
and compare their respective uncertainties due to the PDFs. We see that, at both the Tevatron and the 
LHC, the PDF uncertainties are almost identical for both the gg and bg initial states. 

In Figs. 8.5.39 and 8.5.40 we compare the uncertainties from residual scale dependence and the 
PDFs on the example of bg — > bh (5FNS) at the Tevatron and LHC respectively [40]. Here, we perform 
the comparison for both the total cross section (left) and the total cross section normalized to the central 
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Fig. 8.5.39: Comparison between theoretical uncertainties due to scale dependence and uncertainties arising from the PDFs at 
the Tevatron for semi-inclusive bh production in the Standard Model. In the right-hand plot, both uncertainty bands have been 
normalized to the central value of the total cross section ao. 



value calculated with CTEQ6M (right). Similar results are obtained in the 4FNS. 

From Fig. 8.5.40 one can see that, at the LHC, the theoretical uncertainty is dominated by the 
residual scale dependence. Due to the large center of mass (c.o.m.) energy of the LHC, the gluons 
and bottom quarks in the initial state have small momentum fraction (x) values and, hence, small PDF 
uncertainties typically in the 5-10% range. 

In contrast, due to the smaller c.o.m. energy, the PDF uncertainties at the Tevatron (Fig. 8.5.39) are 
comparable and even larger than the uncertainties due to residual scale dependence over the full Higgs 
mass range. The smaller c.o.m. energy results in higher-x gluons and bottom quarks in the initial state 
which corresponds to large PDF uncertainties in the 10-30% range. 

8.6 Conclusion 

The NLO cross sections for tih and bbh have been presented for the Tevatron and the LHC with emphasis 
on the renormalization/factorization scale and PDF dependences. 
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Fig. 8.5.40; Comparison between theoretical uncertainties due to scale dependence and uncertainties arising from the PDFs 
at the LHC for semi-inclusive bh production in the Standard Model. In the bottom plot, both uncertainty bands have been 
normalized to the central value of the total cross section ao- 
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9 Heavy-flavor effects in supersymmetric Higgs boson production at Iiadron colliders 
Contributed by: A. Belyaev, S. Berge, P.M. Nadolsky, F.I. Olness, CP. Yuan 



The Higgs sector may be represented by one complex scalar- doublet, as it is economically realized 
in the Standard Model (SM), or by two or more doublets, as it takes place in the Minimal Supersymmetric 
Standard Model (MSSM) and its extensions. An important feature of MSSM is that, for large values of 
tan /?, the Yukawa couplings of the ft-quai^ks to the neutral Higgs bosons mu (where mu = h, H, or 
A) are strongly enhanced compared to the SM bbHsM Yukawa coupling. Consequently, production of 
supersymmetric Higgs bosons in bb fusion can have a large cross section in supersymmetric extensions 
of the Standard Model [221, 222, 35, 223, 224]. 

The partonic processes contributing to the inclusive Higgs boson production with enhanced bbniH 
coupUng are represented by (a) bb — > mu; (b) gb — > mub; and (c) gg — > bbrriH scattering. The three 
processes (a,b,c) all give rise to the same hadronic final states, with two S-mesons appearing in different, 
but overlapping, regions of phase space. The distinction between the three processes depends very much 
on the factorization scheme adopted for the QCD calculation, as has been recently reviewed in Ref. [225]. 

As shown in Refs. [62, 76], the correct model for the transverse momentum distribution of the 
Higgs boson is crucial for unambiguous reconstruction of the Higgs boson mass in the rriH — > rr decay 
channel. It is also important for discriminating the signal events from the backgrounds by examining the 
qx distribution of the Higgs boson in mnbb associated production, followed by mu bb decay [4]. The 
transverse momentum {qx) distributions of Higgs bosons may be sensitive to the mass m;, of the bottom 
quark when qx is comparable to nih- In Refs. [226, 201] , we study the effect of the initial-state multiple 
parton radiation and heavy-quark masses on the transverse momentum distribution in the bb — > niH 
process. Here we summarize the results of those two papers. 

9.1 Transverse Momentum Resummation for Massive Quarks 

The resummed differential cross section for inclusive production of Higgs bosons in scattering of initial- 
state hadrons A and B takes the form [205] 

da f°° hdh 



dQ'^dydq, 



2 



f°° hdh 

/ —— Jo[qTh)W{h,Q,XA-,XB,rni,) + Y{qT,Q,y,mb), (9.1.60) 
Jo 



where y is the rapidity of the Higgs boson, xa,b = Qe^^ / are the Born-level partonic momentum 
fractions, S is the square of the center-of-mass energy of the collider, and Jo(gTb) is the Bessel function. 
The resummed form factor W is given in impact parameter (b) space and factorizes as 

W{h,Q,XA,XB,mb) = I ^ag)e-'5(b.Q.-^) p ./^(x^, b, m^) ^ysCxB, b, m,), (9.1.61) 

where the summation is performed over the relevant parton flavors j and k. Here, ajj] is a product of the 
Born-level prefactors, e"'^^'^''^'"'*) is an exponential of the Sudakov integral 

A(as(ll),mb)lni ^ 



S{h,Q,mb) = / -TTT- A{as{fJ.),mb)ln( ^ ] + B{as{iJ.),mb) 



(9.1.62) 
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Fig. 9.1.41: The bottom-quark distributions Vb/piXjhjmt) in the proton vs. the impact parameter b. The solid and dashed 
curves correspond to the S-ACOT and massless ("ZM-VFN") factorization schemes, respectively. 

with bo = 2e~'^^ 1.123, and Vj/A{x, b, mi) are the b-dependent parton distributions for finding a par- 
ton of type j in the hadron A. In the perturbative region (b^ ^ Ag^^), the distributions Vj/A{x, b, mh) 
factorize as 

Vj/A{x,h,mb)\^^^j^-2 = X] / -T'^j/<^^^/^'^''^f>,fJ'F) fa/Ai^fJ-p) (9.1.63) 

a=g,u,d,... ^ 

into a convolutions of the Wilson coefficient functions Cj jaix, b, mi,, ^f) and /c^^-integrated parton dis- 
tributions fa/A{^j fJ-p)- The Sudakov exponential and b-dependent parton densities resum contributions 
from soft and collinear multi-parton radiation, respectively. Y = PERT — AS Y is the difference between 
the finite-order cross section (PERT) and its asymptotic expansion in the small-gT- limit (ASY). 

The Higgs cross sections depend on the mass mf, of the bottom quark. The distributions Vj /a{x, b, mh) 
for the heavy quarks (j = c, b) cannot be reliably evaluated at all impact parameters if a conventional 
factorization scheme, such as the zero-mass variable-flavor number (ZM-VFN, or massless) scheme, is 
used. The reason is that mf, acts as an additional large momentum scale, which, depending on the value of 
b, introduces large logarithms ln"(mf,b) or non-negligible terms oc (mfcb)". The situation encountered 
here is reminiscent of the heavy-quai^k contributions to the DIS structure functions Q^), which are 
not adequately described by the conventional factorization schemes at either small or large momentum 
transfers (see, for instance, [227]). To work around this complication, Ref. [228] proposed to for- 
mulate the CSS formalism in a general-mass variable flavor number (GM-VFN) scheme [229], which 
correctly evaluates the heavy-quark mass effects at all momentum scales. Among all GM-VFN factor- 
ization schemes, the S-ACOT scheme [229, 230] was found to be well-suited for the efficient calculation 
of the CSS resummed cross sections. In particular, in this heavy-quark CSS (CSS-HQ) formalism [228] 
the dependence on m^ is dropped in all 0{as) terms in Eq. (9.1.60) except for Vb/A{x, b, mi,). 

The dependence of the bottom-quark paiton density Vi,/p{x,h,mi,) on the impact parameter is 
shown in Fig. 9.1.41. The ZM-VFN parton density Vi,/p{x,h,mi,) is not properly defined below the 
threshold fip = fni, (or above b = bo /mi,). It was continued to large b in the previous calculations 
using an effective "ZM-VFN" approximation described in Ref. [226]. The S-ACOT parton density 
Vh/p^x, b, mf,) is well-defined at all b. It reduces to the ZM-VFN result at b ^ bo/mi, and is strongly 
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Fig. 9.1.42: Transverse momentum distribution of on-shell Higgs bosons in the bb 
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ruH channel at (a) the Tevatron and 



(b) LHC. The solid (red) lines show the qr distribution in the massive (S-ACOT) scheme. The dashed (black) lines show the 
distribution in the massless ("ZM-VFN") scheme. The numerical calculation was performed using the programs Legacy and 
ResBos [231, 232] with the CTEQ5HQ1 parton distribution functions [233]. The bottom quark mass is taken to be rrib — 
4.5 GeV. 



suppressed at b 3> bo /mi,. The suppression is caused by the decoupling of the heavy quarks in the 
parton densities at /Ui? much smaller than (b much larger than bo/mf,). Consequently the impact of 
the non-perturbative contributions from b > 1 GeV^^ is reduced in the heavy-quark channels compared 
to the light-quark channels. 

The massless ("ZM-VFN") calculation therefore underestimates the true behavior at b > 0.1 GeV~ 
and small qx- This effect can be seen in Fig. 9.1.42, which displays da/dqT for bb mu boson produc- 
tion at (a) the Tevatron and (b) LHC.'^ At the Tevatron, the qx maximum shifts in the "ZM-VFN" approx- 
imation to larger qx by about 2 GeV out of 11.7 GeV (about 17%). For a Higgs mass Mh = 200 GeV, 
the maximum of da/dqT shifts by about 1.9 GeV out of 12.7 GeV. At the LHC, the difference between 
the "ZM-VFN" and S-ACOT calculations is smaller compared to the Tevatron, because the influence of 
the b > 0.1 GeV^^ region is reduced at smaller momentum fractions x probed at the LHC [234]. The 
maximum of the qx distribution shifts in the "ZM-VFN" approximation by about 1.3 GeV (9% out of 
14.1 GeV) to larger q^. The results for other Higgs masses can be found in Ref. [226]. 



9.2 Numerical Comparison with PYTHIA 

The full qx dependence of the bb rriH process is affected by constraints on phase space available for 
QCD radiation (less relevant at small qx). We illustrate the interplay of various effects by comparing 
the CSS-HQ resummation to the PYTHIA Monte Carlo program [235]. We focus on production of the 
CP-odd Higgs particle A for tan /? = 50 (predictions for the other Higgs bosons can be obtained by 
rescaling the bbA coupling). 

As compared to the CSS-HQ formalism, the PYTHIA calculation does not include contributions 

'^Fig. 9.1.42 does not specify the overall normalization of qT distributions. It is valid for both Standard Model and super- 
symmetric Higgs bosons, since at leading order the supersymmetric result can be obtained by rescaling the Standard Model 
bbil/s A/ coupling: Qi/l^^jj A} ~ {— sin a, cos q, sin/375}(;fi,^/cos /3. The net effect of mj, on gr distributions will be the 
same for the SM and MSSM Higgs bosons, up to an overall normalization constant. 
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Fig. 9.2.43: qr distributions for production of 100 GeV CP-odd Higgs bosons A via bb fusion in the Tevatron Run-2. The solid 
and dashed curves correspond to the lowest-order Vl^-term W{1, 1, 0) (with functions A{as{jl)) and B{as{fi)) evaluated at 
0{as)) and PYTHIA. 

generated from the C-functions and y-term, and it evaluates the soft parton contributions at 0{as). 
Therefore, we start by comparing the PYTHIA qt distribution to the resummed VF-term 1, 0) in 
Eq. (9.1.60), with the functions A, B, and C in Eqs. (9.1.62), (9.1.63) being evaluated at orders a^, as, 
and a^, respectively. The orders of as in A, B, and C are shown as the arguments of W{1, 1, 0). 

It is evident from Fig. 9.2.43 that the shapes of W{\, 1, 0) and PYTHIA qt distribution are very 
different, though the integrated rates {i.e., the areas under the two curves) are about the same. The qt 
distribution from PYTHIA is narrower and peaks at lower qx than 1, 0). The large discrepancy 
between the two curves is in contrast to the case of W and Z production via light-quark scattering, 
where the above two calculations predict similar, though not identical, qx distributions [231]. 

A closer examination reveals that additional features must be implemented in the resummed cross 
section in order to reliably describe the distributions of Higgs bosons produced via hh fusion. 

• The kinematical effects account for a large part of the disparity between 1, 0) and PYTHIA. 
The bottom-quark PDF is a rapidly decreasing function of x in the probed range of x. Conse- 
quently, approximations for the true partonic kinematics (especially those made for the light-cone 
momentum fractions x) may have a strong impact on the rate of hh scattering. This feature should 
be contrasted to the behavior of the light-quark PDF's in W and Z production, which include a 
substantial valence component and vary slower with x. As a result, the kinematical approxima- 
tions ai^e less consequential in the W and Z case. 

When PYTHIA generates QCD radiation, the kinematical distributions of the final-state parti- 
cles, including the quarks and gluons from the QCD showering, are modified to satisfy energy- 
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Fig. 9.2.44: Comparison of qr distributions predicted by TOT(l), PERT(l) and PYTHIA, for Higgs boson produced via bb 
fusion at (a) the Tevatron Run-2 and (b) LHC, for Ma =100 and 300 GeV respectively. 

momentum conservation at each stage of the showering. In the resummation calculation, informa- 
tion about the exact parton kinematics is included in the finite-order term (PERT). The resummed 
cross section is therefore expected to be closer to PYTHIA once the 0{as) finite term, PERT(l)- 
ASY(l), is included. In the W{\, 1, 0) calculation, the emitted gluons ai^e assumed not to carry any 
momentum at all in the soft limit. To compensate for small, but nonzero energy of the soft gluon 
emissions, we introduce a "kinematical correction" (KC) in the W and ASY terms. This correction 
modifies the minimal values of partonic momentum fractions xa and xb in order to account for 
reduction of phase space available for coUinear QCD radiation at large qt- 
• The lowest-order cross section W{\, 1,0) does not evaluate effects of the bottom-quark mass, 
which is first included in the C-function of order a^. Also, additional, though not complete, 0{a1) 
contributions arise in the Sudakov form factors inside PYTHIA when the next-to-leading order 
PDF's ai^e used. To account for both features, we evaluate the W term at one higher order (2,2,1) 
and include the m;, dependence using the CSS-HQ scheme. 

Thus, our full prediction TOT(l) is obtained by adding W^^^""^(2,2,l) (evaluated in the CSS-HQ 
formalism with the kinematical correction) and PERT(l), and subtracting ASYkc(I)- It is shown for 
Ma = 100 GeV at the Tevatron in Fig. 9.2.44(a) and Ma = 300 GeV at the LHC in Fig. 9.2.44(b). 
TOT(l) (solid line) is compared to the fixed-order prediction PERT(l) (dashed) and the PYTHIA pre- 
diction (dot-dashed). As one can see, the results for Tevatron and LHC are qualitatively similar. TOT(l) 
is closer to the PYTHIA prediction than W{\, 1, 0), though the two distributions are not identical. The 
PYTHIA qx distribution peaks at lower qx than TOT(l). In the large qT region, the TOT(l) rate is lai^ger 
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Fig. 9.2.45: Comparison of tlie integrated rates, deduced form Fig. 9.2.44, as a function of tfie minimal qr value taken in the 
integration over gr at the Tevatron Run-2 (left) and LHC (right) for Ma =100 and 300 GeV, respectively. 

than the PYTHIA rate. 

Finally, Fig. 9.2.45 shows the integrated cross section as a function of the minimal qt in the 
calculation for the Tevatron (left) and LHC (right). This is another way to illustrate the differences in the 
shapes of qt distributions obtained in the resummation, fixed-order, and PYTHIA calculations. 

9.3 Conclusion 

Multiple parton radiation in 6-quark scattering is conspicuously sensitive to effects of large bottom-quark 
mass TTife and phase-space constraints on coUinear emissions. Both jn^ dependence and phase-space 
dependence tangibly modify the shape of Higgs qx distributions in the bb — > inn processes. The two 
types of effects were consistently implemented within the CSS resummation formalism for heavy-quark 
scattering [228, 226, 201], realized in a massive (GM-VFN) factorization scheme. These corrections act 
on different qx regions. When the dependence on mt, is taken into account, the position of the peak in the 
da I dqr distribution shifts to a lower qr value, leaving the rate at large qx essentially unchanged. The 
kinematical correction is effective in the high-gy region, where it largely reduces the Higgs production 
rate. As a result, we obtain an improved prediction for the full qx spectrum of Higgs bosons, an important 
piece of information needed for the future Higgs searches. 
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10 Higgs Signal for h ^ aa at the Tevatron 

Contributed by: M. Carena, T. Han, G.-Y. Huang, C.E.M Wagner 

The elucidation of the mechanism leading to the origin of mass of all observed elementary particles 
is one of the main goals in high energy physics. The simple Standai^d Model (SM) picture, based on the 
spontaneous breakdown of the electroweak symmetry by the vacuum expectation value of an elementary 
Higgs field, seems to lead to a picture that is consistent with all experimental observables, provided the 
Higss boson mass is smaller than about 250 GeV. Moreover, the best fit to the precision electroweak 
observables measured at the LEP, SLC and Tevatron experiments lead to values of the Higgs mass of the 
order of or smaller than the present bound coming from direct searches at LEP, mn^^ ^114 GeV. 

In spite of the extraordinary good agreement of the experimental observations with the standard 
model predictions, there are many theoretical motivations to go beyond the Standard Model description. 
Several extensions of the Standard Model exist in the literature, and in most of them the Higgs sector 
is extended to a more complicated structure, often including at least two Higgs doublets. The require- 
ment of preserving the good agreement with experimental data can be easily fulfilled in extensions, like 
supersymmetry, in which the effect of the additional particles on the precision electroweak observables 
rapidly vanish with increasing values of the new particle masses. Independently of the particular exten- 
sion, the direct and indirect limits on the Higgs mass must be revised. In particular, the direct search for 
Higgs bosons may be affected by additional decay modes that are beyond the ones analysed by the LEP 
collaborations. 

As an example, let us consider the minimal supersymmetric extension of the Standard Model 
(MSSM). In the MSSM, there is an additional Higgs doublet, leading, in the absence of CP-violation 
in the Higgs sector, to two CP-even and one CP-odd Higgs boson states. At large values of tan /3, the 
ratio of the v.e.v. of the two Higgs doublets, one of the CP-even Higgs bosons acquires Standard Model 
properties, while the second Higgs boson may only be produced in association with the CP-odd Higgs 
boson state. In addition, the masses of the non-standard CP-even Higgs and the CP-odd Higgs are close 
to each other. Under these conditions, the mass bound on the SM-lrke CP-even Higgs is similar to the 
SM one, while the CP-odd and the second CP-even Higgs boson mass bound reads nih > 90 GeV [236]. 

In this note, we will depart from these simple assumptions, by breaking the mass relations that ap- 
pear in the simplest supersymmetric models, and studying the consequences of such modifications of the 
parameters of the theory. Indeed, while it has been a common belief that the Higgs boson will be eventu- 
ally discovered at the upcoming LHC experiments, one would like to fully utilize the potential to search 
for the Higgs bosons at the Tevatron in these non-conventional scenarios as well. Non-standard mass re- 
lations are already present in extensions of the MSSM including an additional singlet (NMSSM) [237], 
or when explicit CP- violations exist in the Higgs sector [238]. In these cases, the SM-like Higgs (h) 
may dominantly decay into a pair of lighter Higgs (a), often the CP odd state. Therefore it is possible 
that the Higgs escaped detection at the LEP experiments by avoiding the usual decay modes such as 
h 2b, 2r, WW* and ZZ*, and the lower limit on Higgs mass should be re-evaluated [239]. We 
are interested in analysing the sensitivity of the Tevatron experiments in the search for a light, SM-lrke 
Higgs boson with such an exotic decay mode. In particular, we will consider the case when the Higgs 
boson decays into a pair of scalars h — > aa, which in turn cascade into a heavy fermion pair a ^ bb and 
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a TT, respectively. 

The dominant production of the Higgs boson at hadron colliders comes from the gluon fusion, 
but the above channel would encounter huge SM backgrounds. We therefore consider the Higgs signal 
produced in association with a.W or Z boson, in the hope that the leptonic decays of the weak bosons 
will provide a clean trigger, and will significantly reduce the background as well. The events being 
searched are 

Wh^ luuaa (10.0.64) 
y liyi,bb,TT 

Zh^ in-,aa ^J^^^j^^'^^ (10.0.65) 

with I = e,iJL. 



10.1 Parameter Choices 

We would like to perform a relatively model-independent search for the typical signal. The direct search 
for a Higgs boson with SM-like couplings to the gauge bosons, in a model and decay mode-independent 
way, leads to a lower bound on rrih of about 82 GeV [240]. Moreover, the proposed search is expected 
to become inefficient for nih > 130 GeV, since the standard decays into the WW* and ZZ* channels 
are still expected to be dominant. Therefore, the optimal setting to detect the Higgs decaying into an aa 
pair is to have the mass rrih within the range of 90 — 130 GeV. The choice for ttIq can be more flexible. 
As long as nia > 2m{, to kinematically allow the decay a — > bb, our analyses are rather insensitive to the 
mass choices (see below for a more detailed analysis of this question). 

In a generic model, the Wh/Zh production rate differs from that in the SM. The change can be 
characterized by a prefactor K^^y^^ (^hzz^^ where Khvv is the coupling strength of Higgs to vector 
boson V relative to that in the SM. The production cross section can thus be written in terms of the SM 
result with an overall factor to account for the modification of the coupling 

a{Vh) = Klyya^^\Vh). (10.1.66) 

We are interested in the range of k ~ 0.5—1.0, so that this Higgs contributes to the electroweak symmetry 
breaking and consequently the associated productions are still sizeable. 

In order for the h ^ aa decay to be dominant and thus escape the LEP bound, BR{h — > aa) is 
required to be close to unity. For instance, in the NMSSM, BR{h — > aa) > 0.9 turns out to be very 
general in terms of the naturalness of c in the trilinear coupling term {cv/2)haa [241]. Moreover, if the 
down quark and lepton coupling to the Higgs is proportional to their masses, then BR{a bb) and 
BR{a Tf) are set to be 0.92 and 0.08, respectively. In general, however, the relations between the 
coupling and the masses may be modified by radiative corrections, which can lead to a large increase of 
the BR{h tt) [77]. The representative values and the ranges of the parameters are summarized in 
Table 10.1.8. 

Including the decay branching fractions, for instance for ai bb, 02 tt, we obtain the cross 
section as 

o-(aa) = ulyy a^^\Vh) BR{V) 2BR{h aa)BR{a bb)BR{a Tf). (10.1.67) 
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representative 


considered 




parameters 


value 


range 


masses 


mh 
ma 


120 
30 


90-130 
20-40 


coupling 


l^hVV 


0.7 


0.5-1.0 


branching 
fractions 


BR{h aa) 
BR{a bb) 
BR{a Tf) 


0.85 
0.92 
0.08 


0.8-1.0 
0.95-0.70 
0.05-0.30 


252r 




0.061 


0.019-0.42 



Table 10.1.8: Parameter choices for h —* aa decays. 

where BR{V) = 0.213 (0.067) is the leptonic branching ratio of the decay of W (Z). into I = e, ji. 
The overall factor modifying the SM result in Eq. (10.1.67), i.e. 

= 2nlyyBR{h aa)BR{a bb)BR{a tt), (10.1.68) 

corresponds to the process-dependent factor defined in the DELPHI search [239]. Our parameter 
choice (range) is equivalent to a C|^2r 0.061 (0.019—0.42), consistent with the bounds for a large 
range of our mh,rna choices set forth in Ref. [239] A value of 0.061 for is assumed for all 
numerics from here on, unless explicitly noted otherwise. 

10.2 Signal Event Rate 

The associated production of pp — > Wh usually features a lai^ger cross section than that of Zh, and the 
leptonic branching fraction of W is about 3 times larger than Z's. For illustration purposes, we choose 
to concentrate on the Wh channel henceforth. 

The Standard Model rate of a Higgs produced in association with a leptonically decaying W is 

a^^^{Wh) BR{W Ivi) ~ 85 (24) fb (10.2.69) 

at ^ = 1.96 TeV for nih = 90 (130) GeV. 

Including the branching fractions and couplings, the cross section of the signal in Eq. (10.1.67) is 

cr{aa) ~ 5.3 (1.5) fb for = 0.061 (10.2.70) 

as illustrated in Fig. 10.2.46. The solid curve on top represents the total cross section for Vh production, 
with V decaying leptonically, but without any cuts. The dashed curve represents the cross section after 
adjusting for the couplings and branching fractions, as in Eq. (10.1.67). Cross sections for Zh are also 
plotted for completeness. 

The events have yet to pass the acceptance cuts, or to have the taus and 6's tagged. Both bring 
significant reductions to the event rate. Our challenges are to retain as many events as possible, and to 
control the backgrounds from various sources. 

'^Conversion between C|f,2T- and C%, involves a factor BR{a bb)/2BR{a tt). 



85 



Wh Production/Detection at Tevatron Zh Production/Detection at Tevatron 
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Fig. 10.2.46: Cross sections of Higgs signal at the Tevatron in the 262r channel produced by Higgs-strahlung with a leptonically 
decaying W (left) or Z (right). rUa = 30 GeV is assumed except for the two curves where ma ~ (mh — 10 GeV)/2. 
= 0.061 is understood. 



10.3 Background and Acceptance Cuts 

We look for events with 5 particles plus missing energy in the final states: hhrflvi. We wish to trigger 
the events by the isolated lepton I, tag the &'s and r's, and demand significant missing transverse energy 
{]^t) in the events. With neutrinos in the decay products, tau momenta cannot be fully reconstructed. 
Therefore we cannot reconstruct the invariant masses or rrih ~ rnhh-rr- Instead, the signal should 
appear as a peak in the mf,{, plot, around the value of nia- 



Acceptance Cuts 

The following cuts are employed to mimic the detector acceptance: 



PT 


> 


10 GeV 


for b, 




> 


10 GeV 


for Wh only 


H 


< 


3.0 


for 6, 


AR 


> 


0.4 


for hh, bl^ 




> 


20 GeV 


for bb, 
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and 



PT > 10, 8, 5 GeV 
\7]\ < 1.5 
AR > 0.4 



for Th,Te,T^ 

for T 

for TT, rb, tI^ 



(10.3.72) 



> 10 GeV 



for TT. 



where re,r^ and Th stand for the decays of r — > eUgUr, r ^J-f^fr, and r — > hadrons + i/,-, respectively. 
Lower cuts on hb and rr invariant masses are to eliminate the large number of background events from 
soft photons and gluons. 

The momentum of the tau-lepton cannot be fully reconstructed since all tau decays involve at least 
one neutrino, therefore the cuts on tau are applied to the visible decay products, and are decay-mode 
dependent. After these acceptance cuts, 10 — 25% of the signal events survive, and the cross section 
becomes 0.6 (0.4) fb for mh = 90 (130) GeV with the given set of input parameters, or ^ 0.06. 
The cross sections passing acceptance ai^e plotted in Fig. 10.2.46 versus the Higgs mass, represented 
by the circled curve. There would be a few events to several tens of events with a few fb^^ integrated 
luminosity, for ~ 0.019 - 0.42. 

Irreducible Background 

The dominant source of the irreducible background the bb pair from a virtual gluon splitting, the rf pair 
from an intermediate Z* /^* and the charged lepton plus missing energy from a W boson. Our Monte 
Carlo simulations with MadEvent [242] show that the Z* is almost on-shell. It can be readily removed 
with a cut on rr invariant mass. However as we shall see below, due to the small size of the background, 
we can afford not to do so. A rr pair from a virtual photon can be more easily confused with the signal, 
but such a background is further suppressed by the electromagnetic couplings. 

After applying the acceptance cuts, the irreducible background is estimated to be around 0.01 fb, 
which is very small compared to the signal size. It is entirely absent given the luminosity at the Tevatron. 

Reducible Background and Tagging 

Taus and &'s need to be identified. During the identification (tagging), signal events ai^e lost due to the 
tagging efficiency and to additional cuts. For example, the efficiency of tagging a single bottom is ai^ound 
50% in the region pT > 15 GeV, and falls off rapidly as we approach lower px- Tagging (hadronic) taus 
faces the same problem. Therefore tightening the kinematical cuts on 6's and taus are necessary to assure 
favorable tagging efficiencies. We decide to tag one b and one tau. 

The b- and hadronic r-tagging efficiencies are taken to be [243, 244] 



Outside these kinematical regions, the tagging efficiencies drop sharply [243, 244]. 

Reducible background arise from jets mis-identified as 6's, or hadronically decaying taus. The 
mistag rate per jet is taken to be around 0.5 — 1.0% (0.5%) for tau (b) [243, 244]. In addition, the 
experiments cannot distinguish directly produced electrons (muons) from leptonically decaying taus. 



b — tagging : 50% for E-;^*' > 15 GeV and \r]jet\ < 1-0 
T- tagging: 40% for E^s > 20 GeV and \r]\ < 1 .5 . 



(10.3.73) 
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• The background due to misidentified bottom comes from the process 2r2j7 + which has a 
cross section of 5 fb. Considering the mistag rate and the additional cuts, it contributes 0.01 fb to 
the background events. 

• The background due to misidentified r differs from the decay modes of r's: 

- For TiTh2bl^T C^l'^bTh^r), it comes from 2r26j with from the leptonic decays of both 
taus. The contribution is estimated at 0.003 fb. 

- For T}iTh2bl^T, the background comes from 2j2bl^T estimated at 50 fb after events of the 
bb and jj resonances around the Z mass ai^e rejected. It's further reduced by a factor of 
0.01 — 0.02 from the tau-mistag rate, and a factor of 0.8 due to b-tagging. 

This results in a background rate of 0.4 — 0.8 fb. In a continuum distribution of mbh, it is at 
or below the level of the signal. We notice that the 5 jets are harder in this background than in 
the signal (see Fig. 10.3.47). Imposing a upper px cut of 50 GeV will reduce the background 
by a factor of about 4, while the signal is minimally affected. 

• The backgrounds from both a mistagged tau and a mistagged b mostly come from the ^jl^r 
events, which has a cross section of about 16 pb. After the cuts and folding in the mistag rates, 
this contributes 0.3 — 0.6 fb of background events. It can be further reduced by imposing upper 
Pt cuts, similarly to the 2j2bl^T background. 
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Background P Distribution 



10 



b or b 
X or X 



o 



■B 



10- 



10-^ 




100 




Fig. 10.3.47: Transverse momenta distribution of the b, r or jets in the Wh signal (left) and background (right) events. 

After carefully tightening the cuts, the reducible background can be a factor of a few to ten smaller 
than the signal, but unfortunately, the cuts and the tagging efficiencies together reduce the signal greatly 
to about 0.08 fb for Wh and 0.03 - 0.05 fb for Zh, with ~ 0.06. With an optimistic ~ 0.42, 
the cross section is 0.55 fb, we would expect to see about a couple of signal events with an integrated 
luminosity of a few fb~^. 

To illustrate a most optimistic situation, we explore the mass relations of and rrih. The signal 
loss is mainly due to the softness of the b and r's, therefore most events are rejected from the lower pr 
threshold. Increasing would stretch the pt distributions to the higher p^ end. To achieve this without 
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significantly affecting the decay phase space of h, we set 

ma = {mh-WGeY)/2, (10.3.74) 

which resulted in almost doubling the signal rate, as seen from the curve with triangles in Fig. 10.2.46. 
The cross section (~ 0.2 fb for = 0.06, and ~ 1.4 fb for = 0.42) is still challenging for 
observation with the Tevatron's projected luminosity. 

10.4 Summary 

The search for a Higgs boson with couplings to the gauge bosons of the order of the SM-one, and 
decaying into two lighter CP-odd Higgs bosons states may be performed at the Tevatron collider in the 
Wh lvaa{2b2T) channel. For 90 < nih < 130 GeV, we found a sizable number of events, with 
negligible in^educible background. However, further cuts and tagging h and r, necessary to remove 
the much larger reducible background worsen the signal event rate to the order 0.1 fb for a value of 
~ 0.06, a factor determined by the product of the relevant branching ratios times the ratio of the 
Higgs production cross section to the SM-one. Therefore, the signal observation becomes statistically 
limited. For an optimal choice of the value of the CP-odd Higgs mass m^, the signal rate may be twice 
as large. With favorable couplings and branching fractions, the factor can be as large as 0.42, which 
would enhance the signal rate by a factor of 7 to around 1 .4 fb. It can be further improved by another 
40 — 60% by combining Zh events with the Wh events, leading to a possible observation of a few events 
for a Tevatron luminosity of the order of a few fb~^. 



89 



11 The Z ^ bb decay as a b-jet energy calibration tool 
Contributed by: T. Dorigo, J. Donini 

We use a sample of Z ^ 66 decays free collected by the CDF experiment to assist in both a precise 
measurement of the energy scale of 6-quark jets and a determination of the 6-jet energy resolution. The 
more precise determination of 6-jet energy scale helps all precision measurements of the top quark mass 
and a determination of the 6-jet energy resolution is important for the search of a low-mass Higgs boson. 
This technique could also prove useful at the LHC. 

11.1 Introduction 

Since their discovery in 1983[245, 246, 247, 248], W and Z bosons have been studied at hadronic 
colliders only using their leptonic decays. As a matter of fact the hadronic decays of these particles are 
generally so difficult to separate from the huge QCD background that, after the extraction of a nice mass 
bump in the dijet mass distribution by the UA2 collaboration in 1987[249, 250], little more has emerged. 

In fact, at the Tevatron things are more complicated for the direct observation of hadronic decays of 
vector bosons. With respect to the SppS, the Tevatron's higher center-of-mass energy is a disadvantage 
for once, because in the face of a four-fold increase in signal cross section, the iiTcducible background 
from QCD processes yielding jet pairs increases by over an order of magnitude, due to the steeply falling 
gluon PDF g{x). 

Using Run I data, hadronic W decays were successfully used by the CDF and DO experiments in 
the discovery and measurement of the top quark both in the single lepton and fully hadronic final state; 
the W — > qq' decay was used both in the event selection and as a constraint in kinematical fits to extract 
the top quark mass. A handful of dijet masses peaking at 80 GeV were also directly observed in a subset 
of high-purity events[251]. In Run II, with increased sample sizes, it has become possible to exploit 
the hadronic decay of W bosons in top events even more, by using them for a direct calibration of the 
energy measurement of light-quark jets in the reconstruction of the tt decay[252]. That technique has 
allowed a significant reduction of the systematic uncertainty arising from the knowledge of the jet energy 
scale, which is by now the largest contribution to the top mass measurement error. 

For the Z boson, which is not produced in top decays and whose inclusive cross section in pp 
collisions is three times smaller than that of the W, the extraction of hadronic decays is even more 
complicated; only the decay to 6-quark pairs reaches the level of observability, thanks to the significant 
reduction of QCD processes provided by the distinct signature of 6-quark jets. Indeed, a small signal of 
Z ^ bb decays was extracted by CDF in Run I data exploiting the semileptonic decay of 6 quarks with an 
inclusive muon trigger of low Pr[253]. The signal was too small to allow any study of 6-jet energy and 
resolution, but its demonstrated observability in the Tevatron environment gave hope to the searches for 
the analogous signature of a low-mass Higgs boson decay, and spurred the development of a dedicated 
trigger for Run II, capable of collecting a large Z signal without the need to rely on the semileptonic 
decay of 6 quarks. 

A large-sized signal of Z — > 66 decays free from selection biases allows both a precise measure- 
ment of the energy scale of 6-quark jets and a determination of the 6-jet energy resolution. The reduction 
of the uncertainty in the 6-jet energy scale helps all precision measurements of the top quark mass, while 
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a determination of the 6-jet energy resolution is important for the search of a low-mass Higgs boson. The 
signal, most notably, opens the doors to a direct test of algorithms that attempt to increase the resolution 
of the 6-jet energy measurement. These algorithms are a critical ingredient for the observability of the 
Higgs boson at the Tevatron if Mh < 135 GeV. 

11.2 Triggering on Z ^ bb decays 

In Run II CDF benefits from a hai^dwai^e tracker using silicon detector hits at the second trigger level, the 
Sihcon Vertex Tracker (SVT)[254]. 

The SVT works by comparing the pattern of hits in the five layers of silicon sensors of the Sil- 
icon Vertex Detector (SVX)[255] to those expected by charged tracks of given transverse momentum, 
azimuth, and impact parameter, which are stored in 256 associative memory chips. Use of the 12-fold 
azimuthal symmetry of the SVX reduces the number of needed patterns and allows a parallelization of 
the task of finding track candidates and performing linearized fits. On average as little as 15 fis are 
needed to process an event and determine the impact parameter of tracks with a resolution of 35 /im. 
The efficiency to reconstruct fiducial tracks with Pt > 2 GeV is close to 90%. 

Using SVT information as well as calorimetric input, the Z bb trigger selects events containing 
two back-to-back Et > 10 GeV jets and two Pt > 2 GeV ti'acks whose impact pai^ameter with respect 
to the beam line is larger than 160 fim; a veto on forward jets with Et > 3 GeV is also applied to 
reduce QCD backgrounds. These requirements have an efficiency of about 4% on Z — > 66 decays, and 
they result in an effective cross section lower than 100 nb, which corresponds to a manageable rate for 
machine luminosity up to L = 10'^^ cm~^s~^. 

As the luminosity grows, so does the average number of multiple interactions occurring in the 
same bunch crossing. Since the Z signal can only be isolated in clean events with two back-to-back jets 
and little extra jet activity, it is reasonable to foresee a dynamic prescaling of the trigger, which should 
anyway allow the collection of at least 2 fb^^ of data with the base data collection plan of Tevatron's 
Run II.A sample of 80,000 signal events is thus achievable. 

11.3 Preliminary Run II results 

A signal of Z decays to 6-quai-k pairs has been observed in 333 pb~^ of CDF data collected by the trigger 
described above. After a reconstruction of jets with a.R = 0.7 cone algorithm[256], events were selected 
by requiring two jets of raw transverse energy exceeding 20 GeV in the rapidity interval \ri\ < 1.5, both 
of them containing a secondary vertex (6-tag) reconstructed by the SecVtX algorithm[257]. 

After those requirements the Z signal is still buried in a very large background consisting predom- 
inantly of QCD direct bb production, which needs to be reduced further. 

Most direct bb pairs are produced at the Tevatron by gluon fusion, whose high color charge in the 
initial state and color flow topology are distinctive characteristics. To exploit the smaller probability of 
QCD radiation from the initial state quarks in Z boson production, the two leading jets are required to 
be back-to-back in azimuth within A<l>jj > 3, and events containing a third jet with raw Eij, > 10 GeV 
are discarded. 

The resulting sample of 86,000 events contains roughly 3400 Z boson decays. Their reconstructed 
dijet invariant mass can be fit using as a background template the mass distribution of dijet events which 
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L7-corrected Dijet Mass (GeV) Data/MC JES scale factor 



Fig. 11.3.48: The Z bb signal extracted by CDF with 333 pb^^ ofRun II data. Left: the dijet mass of events with two 6-tags 
is fit as the sum of a background template (in grey) and a signal template (in green). Top right: fit as a function of the fe-JES. 
Bottom right: number of signal events from the fit as a function of the 6- JES. 

do not contain secondary vertices, by accounting for the bias due to the non-flat 6-tag probability versus 
dijet mass using a correction function; the latter is obtained from events failing the kinematical require- 
ments. For the Z signal, 40 Monte Carlo-derived templates with a varying 6-jet energy scale factor (JES) 
from 0.8 to 1.2 in steps of 0.01 are used in turn. One thus obtains the dependence of fit on the 6- JES, 
from which a measurement of the latter can be extracted. Fig. 11.3.48 shows the results for the best fit, 
which corresponds to a 6- JES of unity, with a statistical uncertainty smaller than 2%. 

11.4 Prospects for the B-jet energy scale extraction 

The largest contribution to the total uncertainty in the top quark mass determination at the Tevatron 
originates from the knowledge of the jet energy scale, a factor which measures the discrepancy between 
the effect of detector response and energy corrections in real and simulated hadronic jets. The JES 
can be determined from studies of photon-jet balancing[258], but modeling and selection biases limit 
the accuracy of the method; a determination which is mostly statistics-limited comes instead from the 
measurement of ^ qq' decays in top events. With these methods, the Tevatron experiments can reach 
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a precision close to 1% on the top quark mass in Run II by reducing the JES uncertainty to a similar 
level. 

When dealing with 6-jets, however, one has to cope with several peculiarities of their fragmentation 
and decay properties, and with the different color connection of 6-jets and light quark jets in top quark 
decay. All these effects have to be accurately modeled if one is to use a generic JES factor extracted from 
jets not containing heavy flavors to the two 6-jets always present in a tt decay. A recent study predicts 
that the uncertainties in the modeling of fragmentation, decay characteristics and color flow may affect 
the knowledge of the 6- JES by as little as 0.6%[259], but a direct determination of that quantity is of 
course preferable. 

Due to the small cross section of production processes yielding events with a high-energy pho- 
ton recoiling against a 6-quark jet, a measurement of the 6- JES with balancing techniques is difficult, 
although both DO and CDF have recently started exploring that option. 

The preliminary Z —>■ bb signal extracted by CDF appears to demonstrate that the data have a 
sufficient statistical power to allow the determination of a precise 6- JES factor: one expects that 2 fb~^ 
of data will reduce the statistical uncertainty of template fits well below 1%. However, systematics are 
a concern: given the smallness of the signal to noise ratio of the selected sample, a meaningful determi- 
nation of the 6- JES from bump fitting requires that the background shape be modeled with the utmost 
accuracy, especially if its most probable value occurs close to that of the Z signal, as in Fig. 1 1.3.48. 

Reducing the 20 GeV threshold on raw transverse energy of the jets, which directly affects the 
peak position of the background distribution, is however problematic, since 1ow-£'t jets suffer from 
subtle trigger effects which are hard to model correctly; moreover, at very low Et it has been shown 
that cc production is a sizable contribution to the S VT-triggered dataset, and its presence complicates the 
understanding of the 6-tagging bias in the data. 

Recent studies have shown that using large Monte Cai^lo samples of direct bb production processes 
and careful parametrizations of trigger biases it is possible to check the background shape extracted from 
the data as described above, and reduce the associated systematic uncertainty in the determination of the 
6- JES. CDF plans to use that additional information to finally achieve a significant measurement of the 
6-jet energy scale with Z — > 66 decays, thus justifying several years of investigation of the extraction of 
a well-known signal. 

11.5 B-jet energy resolution studies 

The mass resolution of pairs of 6-jets has been duly stressed as one of the critical factors in the search for 
a light Higgs boson decay at the Tevatron. While the 1999 study of the Tevatron Higgs Working Group[4] 
could only make the educated guess that a aM/^hh ~ 10% relative mass resolution was attainable with 
a dedicated effort, the Higgs Sensitivity Working Group[132] went as far as producing some evidence 
that such precision was indeed reachable, by a careful use of several corrections in series, followed by 
the exploitation of the coiTclations between kinematic variables measured in WH Ivbb events and 
the induced biases in the dijet mass measurement (see Fig. 11.5.49). 

Indeed, when compared to any selection applied on the data to increase the signal to noise ratio, an 
algorithm that reduces the width of a resonance sitting on top of a large background has the obvious ad- 
vantage of keeping intact the size of the signal. If signal significance is on the yardstick, a 20% decrease 
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Fig. 11. 5. 49: Dijet mass distribution for pairs of &-jets after different levels of jet energy corrections: raw jet energies (top left), 
energies corrected with the HI algorithm (top right), then after the subsequent application of &-specific corrections (bottom 
left), and finally after the use of the hyperball algorithm, a method that corrects the dijet mass accounting for the correlation 
with event observables. The red points describe the behavior of WH Ivbb events; the stacked histograms are Monte Carlo 
simulations of Vy+ jets (purple), production (cyan), single top production (green and yellow), and W Z production (blue). 

of aM/Mhb can be shown to have the same effect of a 20% increase in collected luminosity[132]. 

The resolution in the transverse energy of generic jets can be measured with 7 — jet events. 
Those events have in fact constituted the basis of CDF studies of an algorithm exploting both the 
tracker, the shower max detector, and the calorimeter to increase the precision of the transverse en- 
ergy measurement[260]. As Fig. 1 1.5.50 shows, a 30% improvement in the resolution of generic jets can 
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Fig. 1 1.5.50: Jet energy resolution measured in photon+jet data as a function of photon Pt- The blue points show the resolution 
of the standard jet energy corrections, the red points show the results of a dedicated algorithm exploiting information from many 
subdetectors. 



be achieved by a combined use of the information from different subdetectors. Unfortunately, the lack of 
sizable samples of data containing a 6-quark recoiling against an energetic photon prevents a data-driven 
study of the 6-jet energy resolution, and a check of the effect of generic coiTcctions applied to &-jets. 

The development of a 6-specific algorithm, aimed at increasing the Et resolution on jets contain- 
ing 6-tags, highly profits from the availability of a statistically significant sample of Z bb decays: one 
can then both check the effectiveness of any recipe and measure the resulting mass resolution. 

As was noted above, 6-jets are different from generic jets originated from light quarks or gluons 
in several aspects. 23% of 6-quarks decay semileptonically, and more soft leptons are yielded by the 
following charm quark decay; the large mass of the 6-quark produces tracks with significant transverse 
momentum with respect to the jet axis; and finally, ^-quarks have a hard fragmentation function, which 
may translate in an average detector response different from that of generic jets. The total effect of 
these peculiarities is a worse Et resolution for 6-jets and a significant negative bias, mainly due to the 
neutrinos from semileptonic decay. 

S-jets which are tagged by a vertex-finding algorithm are also different from an experimental 
point of view, since the detection of a displaced vertex allows the measurement of several ancillary 
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characteristics: the distance between primary and secondary vertex, the total charge of tracks forming 
the secondary vertex, the total transverse momentum and combined mass of the charged decay tracks. 

All these observable quantities can be exploited by algorithms detecting the correlation between 
their values and the average bias on the jet Et measurement. For instance, the presence of a muon in a 
jet is strongly correlated with the resulting calorimeter response, such that the muon Pt can be used with 
success to increase the Et resolution. The best results are obtained when all correlations are exploited 
together, by finding the most probable bias in the Et measurement as a function of the value of all 
observed jet variables. 



2-component fit to double tagged ZBB data 
• Double tagged events 

Background: 0.9601 81 ± 0.006582 

Z signal: 3394 ± 515 evts 

Fit result: x^/dof= 1 .772352 
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Fig. 11.5.51: Relative Et resolution on 6-tagged jets from a QCD bb Monte Carlo simulation, after generic jet corrections (top 
left, red points) and after fe-specific jet corrections using all jet observables (bottom left, blue points). Right: a comparison of 
the Et resolution obtained with the two corrections. 

Preliminary results by the CDF collaboration have determined that the Et resolution for 6-jets can 
be improved by as much as 30% from the baseline resolution yielded by the application of standard, non- 
6-specific energy coiTcctions (see Fig. 11.5.51). Two algorithms ai^e being developed for that purpose, 
and the study on the Z ^ bb signal will prove their effectiveness in the neai" future. 
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12 Selected Topics in Standard Model Higgs searches using H W^W and H t^t 
Decays at the LHC 

Contributed by: B. Mellado, W. Quayle, S. L. Wu 

We study control samples for an in-situ determination of the major backgrounds to i7 — > W^W~ - 
l^l~vu, where a full jet veto is applied. We find that the theoretical uncertainty on the extrapolation of 
the QCD W~^W~ background from the control sample to the signal-like region is 5 % and that the im- 
pact of the singly-resonant top background on the extrapolation from a b-tagged control sample is about 
10 %. We therefore conclude that it is possible to perform a data-driven estimate of the background 
in the signal-like region with sufficient accuracy to achieve a. 5 a discovery of a 160 GeV Higgs with 
less than 2 fb^^ of integrated luminosity. We evaluated the Next-to-Leading Order corrections to Higgs 
production in the analysis H — > r+r~ in association with one High Pt jet. The ratio of NLO to LO 
cross-sections after the application of analysis cuts is in the range 1.5 1.6 for Higgs production via 
gluon-gluon fusion. 

12.1 Introduction 

We study control samples for an in-situ determination of the major backgrounds to — > W~^W~ — > 
l~^l~i'v, where a full jet veto is applied. We find that the theoretical uncertainty on the extrapolation 
of the QCD W~^W~ background from the control sample to the signal-like region is 5 % and that the 
impact of the singly-resonant top background on the extrapolation from a b-tagged control sample is 
about 10 %. We therefore conclude that it is possible to perform a data-driven estimate of the background 
in the signal-like region with sufficient accuracy to achieve a. 5 a discovery of a 160 GeV Higgs with 
less than 2 fb~^ of integrated luminosity. We evaluated the Next-to-Leading Order connections to Higgs 
production in the analysis H — > t^t~ in association with one High Pt jet. The ratio of NLO to LO 
cross-sections after the application of analysis cuts is in the range 1.5 1.6 for Higgs production via 
gluon-gluon fusion. 

The search for the Higgs boson called for by the Standard Model is arguably one of the most 
important topics in high-energy particle physics today. For the early observation of a Higgs boson of 
mass 115 < Mh < 135 GeV the most relevant final states involve H ^ j^f and H t'^t" [261]. 
For the range of masses 135 < Mh < 190 GeV the most promising decay mode of the Standard Model 
Higgs boson is H ^ W+W [262]. 

In this work, we discuss selected topics related to the search for the Standard Model Higgs boson 
at LHC using H — > W~^W~ and H — > t~^t~. In Section 12.2 we describe our Monte Carlo samples, 
event selection, methods for in-situ background determination for the channel H — > W~^W~ — > I'^l'vv 
with a full jet veto. Section 12.7 reports on studies of QCD higher order coiTcctions to Higgs signal 
production in the analysis of H —>^ t^t~ in association with one high Pt jet [263]. 

12.2 Selected Topics in H ^ W+W' l+l-uu Analysis 

This Section is subdivided into four sub-sections. Sub-section 12.3 describes the Monte Cai^lo samples 
and the analysis method used in the analysis for the search of the Higgs boson with H W^W~ 
l^l^uv when applying a veto on events with high Pt jets. We also discuss data-driven methods for 
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Table 12.3.9: Cut flows (in fb) for Mh = 160 GeV in the e/i channel. 



Sample 


gg^ H VBF 


tt 


EWWW 


gg WW 


qq WW 


Z ^ TT 


Primary 


1.86 0.03 


33.4 


0.08 


6.19 


121.0 


7.96 


b-tagged 


0.18 0.007 


17.02 


0.0001 


0.08 


1.51 


1.29 



Table 12.3.10: Cross-sections (in fb) in the two control samples discussed in Section 12.3 for Ain = 160 GeV, summed over 
lepton flavor. 

the extraction of the backgrounds. In sub-sections 12.4 and 12.5 we discuss the theoretical uncertainties 
in the background extraction procedures. In sub-section 12.6, we perform a brief comparison of three 
generators for the W'^W" background for validation purposes. 

12.3 Monte Carlo and Analysis Method 

We consider the following signal and background processes: 

• Higgs production. We model the gluon-initiated process with the generator provided in MC@NLO 
and normalize the cross-section for the signal to the values obtained used HIGLU [84]. The small 
contribution from Weak Boson Fusion (VBF) is modelled with Pythia [264, 265]. 

• QCD W^W~ production is modelled with the generator provided in MC@NLO version 3.1 [266, 
267]. A non-negligible number of W^W~ events come from gg —>■ W^W~ diagrams that are 
not included in MC@NLO; we model this contribution using the generator documented in [268]. 

• ti production. The (dominant) doubly-resonant contribution is modelled with MC@NLO. To es- 
timate the impact of the singly-resonant and non-resonant W~^W~bb contributions to the back- 
ground, we perform a comparison between leading-order calculations of pp — > W~^W~bb and 
pp^ti^ W+W~bb using MadEvent [269, 270]. 

• QCD Z/7 production, with Z — > ee/ ^ih/tt. We model this background with MC@NLO. 

Although we do not expect detector effects to be important in this calculation, it is convenient to simulate 
a detector using the last fortran-based release of ATLFAST, and we apply the jet energy corrections in 
ATLFAST-B [271].!'^ 

Our event selection consists of the f ollowing cuts: 

'^We also apply a small correction to the energy of jets for which HERWIG was used for the parton showering and hadroniza- 
tion; the correction is given by (1 — 5 x IQ^'^ P^* + 0.042) where the jet Pt is measured in GeV. 
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Fig. 12.3.52: The distribution of the azimuthal angle (in the transverse plane) between the leptons after cuts. 

• Trigger and Topology cuts. We require that the event has exactly two leptons with transverse 
momentum greater than 15 GeV in the region with |?7| < 2.5, and we apply a lepton identification 
efficiency of 90% for each lepton. The dilepton invariant mass is required to be less than 300 GeV. 

• Z rejection. The event is rejected if the leptons have an invariant mass between 82 and 98 GeV. 
We require a large missing transverse momentum Pt > 30GeV which is raised to 40 GeV if the 
two leptons have the same flavor. To reduce the nontrivial background from the decay Z — > rr ^ 
II + Pt, we calculate, using the collinear approximation, x^. and x^, the energy fractions carried 
by the visible decay products of the r leptons, and AI^-t, the invariant mass of the two r leptons. 
We reject the event if > 0, x^ > 0, and \Mrr - Mz\ < 25GeV. 

• Jet Veto. We reject the event if there are any jets with Pt > 30GeV anywhere in the detector, or if 
it contains any b-tagged jets with Pt > 20GeV and \r]\ < 2.5. We assume a b-tagging efficiency 
of 60% with rejections of 10 and 100 against jets from c quarks and light jets, respectively. 

• Transverse momentum of the Higgs candidate. We require that P^'^^" > 11.1 GeV. 

In the signal-like region, we apply three more cuts: we require that the dilepton mass have 6.3 < Mu < 
64.1 GeV, that the azimuthal opening angle between the leptons satisfy A(f)ii < 1.5 radians, and that 
the transverse mass obey 50 < Mt < Mh + 10 GeV. The cross-sections after successive cuts for a 
representative Higgs mass of 160 GeV in the e// channel are shown in Table 12.3.9. We also consider 
two control samples: 

• The primary control sample is defined the same way as the signal-like region, but with different 
cuts on the dilepton opening angle in the transverse plane and the dilepton invariant mass. We 
require A(j)u > 1-5 radians and 80 < Mu < 300 GeV; we remove the cut on the transverse mass. 

• The b-tagged control sample cuts are the same as in the primary control sample, except that instead 
of applying a b-jet veto, we require that there be a b-tagged jet with Pt between 20 GeV and 
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30 GeV; we also remove the lower bound on the dilepton invariant mass. 

Table 12.3.10 shows the cross-sections in these two control samples. In order to make meaningful esti- 
mates of systematic errors, it is helpful to define the following three quantities: 

• Ciww- The ratio of the QCD W'^W~ cross-section in the signal-like region over the QCD 
W^W~ cross-section in the primary control sample. 

• att '■ The ratio of the ti cross-section in the signal-like region over the ti cross-section in the b- 
tagged control sample. 

• : The ratio of the ti cross-section in the primary control sample over the ti cross-section in 
the b-tagged control sample. 

With these ratios taken from Monte Carlo, we estimate the number of ti events in the signal-like region 
as j\[^^9n.ai-iike _ auNi,_tagged ^nd the number of W^W~ background events in the signal-like region 
as 

-n.Tsiqnal—like T\Tprimary / i^Tprimary WW -\t n i. i ,i \ 

^WW =awwN^w -(^WwiNlotai - ottt Nb-tagged - small backgrounds) 

where ]S[^~*°'3a<^d number of events in the b-tagged control sample, Nf^^^"'^^ is the total number of 
events in the primary control sample, and the "small backgrounds" consist mostly of Drell-Yan events. 

12.4 Theoretical Uncertainties in the W^W~ Background 

We begin with the theoretical uncertainties in the extrapolation coefficient aww- Here, the theoretical 
eiTor is dominated by the uncertainty in the normalization of the gg W^W~ contribution; recent 
studies have shown that this contribution can be in excess of 30% for the cuts used in those studies [268, 
272]. 

We compute the the theoretical error as the sum in quadrature of the uncertainty due to the fit error 
in the parton density function parameterization and the uncertainty due to the choice of scale. To 
estimate the parton density function (PDF) uncertainty, we have used the CTEQ6 PDF set and its error 
sets; using equation (3) in [273], we find that the uncertainty in aww is 2.8%. To assess the uncertainty 
due to the choice of scale, we have varied the renormalization and factorization scales by factors of 
8.^^ We examine four choices of scale variations: Scale 1 has Qren 8Qren, Qfac Qfac/^\ Scale 

2 has Qren ^ <5ren/8, Qfac ^ ^Q fac, Scalc 3 haS Q^en ^ ^Qren, Qfac ^ ^Qfac, and Scalc 4 haS 

Qren Qren/8, Qfac ^ Qfac/8. Table 12.4.11 shows the cross-sections before and after cuts in the 
signal-like region and primary control sample for the gg — > W~^W~ and qq — > W^W~ contributions, 
with the central-value Q^ scales and the four modified scale choices. The largest variation in aww we 
observe is 4.1%, and we take this to be the theoretical en^or due to the choice of Q^ scale. The total 
theoretical uncertainty we calculate on the prediction of aww is therefore 5%. 

12.5 Theoretical Uncertainties in the Top Background 

We now turn our attention to the uncertainties in au and a^^ . Here, the most important question to ask 
is how to handle single top production. A procedure for generating both pp — > ti and pp — > Wt without 

'^This is an unusually large scale variation to choose; typically, a scale uncertainty will be quoted based on a scale variation 
of 2 or at most 4. Our motivation for this choice is the fact that we expect the K-factor for gg W^W^ to be large, since the 
K-factor for gg —> 77 has been calculated and it is slightly less than 2 [274]. 
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No cuts 


Sig. Reg. 


Cont. Samp. 




Scale Choice 


gg WW 


qq WW 


99 


qq 


99 m 


aww 


Central 


487.77 


11302.44 


6.45 


63.20 


6.38 


130.10 


0.5103 


scale 1 


239.93 


12862.82 


2.92 


69.25 


3.33 


143.83 


0.4904 


scale2 


1058.97 


9076.86 


14.5 


49.03 


13.46 


107.44 


0.5255 


scale3 


278.17 


11189.52 


3.81 


65.02 


3.54 


131.92 


0.5081 


scale4 


913.38 


11702.80 


11.1 


61.81 


12.66 


133.51 


0.4988 



Table 12.4.11: Cross-sections before and after cuts for the signal-like region and the Primary control sample, with the corre- 
sponding extrapolation coefficients, using the nominal assumptions and the 4 altered scale choices. 



Process 


Signal-like 


Cont. Samp. 


b-tagged 


att 


WW 


W'^W-hh 


13.34 


109.41 


47.13 


0.2829 


2.3211 


tt W+W-bb 


9.80 


80.77 


37.72 


0.2599 


2.1413 



Table 12.5.12: Cross-sections (in fb) and extrapolation coefficients for the tt background for various masses, using MadGraph 
to model the W^W~bb background. 

double-counting at leading order was presented in [275], and a calculation including off-shell effects and 
spin correlations in the W~^W~bb system at tree level was presented in [276]. Unfortunately, we know 
of no event generator available at the time of this writing which also takes into account the one-loop 
radiative corrections to W^W~bb production, so we will perform our uncertainty estimate at tree-level. 

In addition to the it Monte Carlo sample (from MC@NLO) that we have used in the other sections 
of this note, we have generated two separate W~^W~bb Monte Carlo samples using MadGraph. One 
includes only doubly-resonant top quark pair production, and the other includes the full W~^W~bb final 
state. For this generation, we have allowed the b-quarks to be generated with Pt as low as 1 GeV, 
and with pseudorapidity as high as 100. One would expect a disproportionately lai^ge contribution from 
the region where one b-quark is soft or forward, and we therefore feel it is likely that the single-top 
contribution is overestimated in our nonresonant W~^W~bb Monte Carlo. This is exactly what we want 
if we are to prove that our analysis is robust. We have applied the cuts for the signal-like region and both 
of the control regions to these two Monte Carlo samples to assess the importance of single-top production 
in this analysis. 

Table 12.5.12 shows the W~^W~bb background cross-sections in the signal-like region, the pri- 
mary control sample, and the b-tagged control sample obtained with the leading-order doubly-resonant 
ti and inclusive W^W~bb samples. We note that although the difference in the absolute cross-section 
given by the two samples is approximately 30%, the corresponding differences in the predictions of au 
and are only about 9%. It is worth noting that this figure is only a general guideline, since the exact 
values of au and are strongly dependent on the particulars of the b-tagging algorithm used. Our 

intent in this section is only to give a rough idea of what the theoretical uncertainty on the extrapolation 
from a b-tagged sample to a b-vetoed sample should be. In practice, this uncertainty should be addressed 
in detail using full detector simulation by any experimenter performing a. H ^ W'^W~ search like the 
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one outlined here. 

12.6 Comparison of MC @ NLO, Alpgen, and Sherpa 

In this subsection, we check the agreement among the predictions of the QCD W~^W~ background given 
by MC@NLO and by tree-level jet-parton matching algorithms like the ones in Sheipa and Alpgen. We 
begin with a few general comments about the generators under study. In this sub-section we ignore the 
process gg W'^W~. 

• All three generators ignore the contribution from gluon-initiated diagrams that contain a quark 
box. This contribution is not negligible; in practice, we treat the gluon-initiated contribution as a 
separate process modelled with a separate generator. 

• The matrix element calculations in MC@NLO and Alpgen were programmed by hand by their 
respective authors, while Sherpa uses an automated matrix element generator to write code to 
compute the (tree-level) matrix elements relevant to a particular process. There ai^e therefore some 
differences regarding which Feynman diagrams are included in the two calculations. In the case 
of this analysis, where we are concerned with the production of W pairs which decay leptonically, 
Sheipa includes the contribution from diagrams where two Z bosons are produced, with one Z 
decaying to leptons and the other to neutrinos. This leads to a spike in the dilepton invariant 
mass distribution in events with same-flavor leptons; this feature does not appear- to be present in 
MC@NLO and Alpgen. For this reason, we will consider only events with one electron and one 
muon in this section. 

• MC@NLO includes the contribution from loop diagrams in its calculation; Sherpa and Alpgen 
rely instead on jet-parton matching schemes like the one discussed in [277]. 

It is worthwhile to point out that although the treatment of soft hadronic physics in Alpgen, Sherpa, 
and MC@NLO are all quite different, the result is nevertheless similar for the three generators. Fig- 
ure 12.6.53 shows the distribution of the transverse momentum of the Higgs candidate (in the e — fi 
channel) given by Alpgen, Sherpa, and MC@NLO for the QCD W^W~ background. Although a 
detailed study of the eiTors on these distributions is beyond the scope of this work, we feel that the sim- 
ilarity among all three generators is encouraging. We note that the Alpgen and Sherpa samples predict 
a slightly lower cross-section for events with Higgs candidate Pt between roughly 10 and 25 GeV. This 
is no doubt an artifact of the jet-parton matching method, and we expect that the behavior of this region 
could be tuned by tuning the matching parameters in the respective generators (although such a tuning is 
not necessary for our analysis). 

It is also worthwhile to compare the predictions of variables related to spin correlations in the 
W^+W^" system, as these correlations are crucially important for the analysis. Figure 12.6.54 shows 
the distribution of the azimuthal angle between the leptons (in the e — fi channel) before the cuts in the 
next section; there is a slight difference in the shape of these inclusive distributions. The discrepancy is 
not serious at all; we believe it is a kinematic effect caused by the depletion in events with Higgs Pt 
between 10 and 25 GeV mentioned in the previous paragraph. Figure 12.6.55 shows the distribution of 
the dilepton opening angle in the transverse plane after all but the last three cuts of Section 12.3 have 
been applied; there is good agreement for this distribution among the various generators. Figure 12.6.56 
shows the dilepton invariant mass for events with one electron and one muon (before the cuts of the 
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Alpgen WW+0/1j 
Alpgen WW+0/1/2j 
Sherpa QCD WW+0/1j 
MCatNLO qq^ WW 



nl 1 1 1 I 1 I 1 I 1 1 1 I 1 I 1 I I 1 I I 1 . . I . rr.T-rr-T-w 

20 40 60 80 100 120 140 

Higgs Candidate (GeV) 

Fig. 12.6.53: The transverse momentum of the Higgs candidate in QCD W pair production as given by Alpgen, Sherpa, and 
MC@NLO. There is a small shift in the location of the peak, but the difference is not dramatic at all. 




Fig. 12.6.54: The distribution of the azimuthal angle between the leptons (in the e — jj, channel) in QCD W pair production as 
given by Alpgen, Sherpa and MC@NLO. 
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Fig. 12.6.55: The distribution of the azimuthal angle between the leptons (in the e — /i channel) in QCD W pair production 
as given by Alpgen, Sherpa and MC@NLO. This figure plots the distributions after all cuts except the cuts on Mu, A<j)ii, and 
Mt. 




Fig. 12.6.56: The distribution of the dilepton invariant mass (in the e — n channel) in QCD W pair production as given by 
Alpgen, Sherpa and MC@NLO. 
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analysis are applied); it is cleai^ from the figure that the distribution of this variable is also very similar in 
all three generators. 

12.7 QCD NLO Corrections for Higgs Production in H ^ t~^t~ in Association with One High 

In a recent publication it was demonstrated that the search for the Higgs boson using the H — > r+r" 
decay in association with one high Pt jet is a promising discovery channel at the LHC [263]. Next- to- 
Leading order corrections (NLO) were not evaluated for the final state considered in [263] . 

There are two main requirements in the analysis, which are instrumental in achieving a good 
signal-to-background ratio: the appUcation of a large cut on the Pt of the Higgs candidate {Pth > 
100 GeV) and the requirement that the invariant mass of the Higgs candidate and the leading jet in the 
event be very large (Mhj > 700 GeV). In addition to the two cuts just mentioned it was necessary to 
require that there be no additional hard jets (hadron level Pt > 20 GeV) in the central region of the 
detector (|?/| < 2). The latter is introduced to suppress the it background. 

It is meaningful to evaluate QCD higher order corrections to the signal process after the application 
of the cuts mentioned above. Apart from a chance in the overall normalization, the impact of extra jets 
in the final state on the analysis is not expected to be trivial. 

In order to evaluate QCD Next-to-Leading corrections, we use the MCFM program [278]. This 
package enables the user to apply cuts at the parton level. Next-to-Leading Order matrix elements to 
Higgs production in association with one jet are available. In this calculation the infinite top mass ap- 
proximation is used. In addition, NLO matrix elements for Higgs production via weak boson fusion are 
also available within MCFM. 

Figure 12.7.57 shows the Higgs Pt (plots on the left) and the invariant mass of the Higgs and 
the leading jet (plots on the right) for Higgs produced via gluon-gluon fusion. The upper plots in Fig- 
ure 12.7.57 show the distributions to Leading Order (LO, solid lines) and to NLO (dashed lines). The 
lower plots in Figure 12.7.57 show the ratio of the NLO to the LO cross-sections. 

Figure 12.7.57 illustrates that the QCD NLO corrections to the signal produced via gluon-gluon 
fusion in the region of the phase space where the Higgs boson will be searched for are large. The 
size of the NLO correction is larger than the correction before the application of cuts on the Higgs Pt 
and Mhj. The perturbative analysis shows that the NLO correction grows with Mhj- This can be 
understood qualitatively: for large values of Mhj a large Pt extra parton is likely to be present in the 
final state, providing extra transverse momentum to the system made by the Higgs and the leading jet 
and indirectly enhancing its invariant mass. 

As pointed out above in this Section, the application of a veto on extra hadronic activity is impor- 
tant for the suppression of the ti production. It is necessary to evaluate the QCD NLO corrections after 
the addition of this jet veto. Plots in Figure 12.7.58 show the same distributions as in Figure 12.7.57 
after the application of a veto on events with an extra parton with Pt > 30 GeV and \7]\ < 2. The 
plots on the left illustrate that after the application of a veto on extra high Pt partons changes size of the 
NLO corrections takes place. The ratio of NLO to LO cross-sections decreases with the Higgs Pt for 
Pt > 50 GeV. After the application of the veto the ratio of the NLO to LO cross-sections becomes flat 
as a function of Mhj instead of increasing, as illustrated in Figure 12.7.57. 
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Fig. 12.7.57: Higgs Pt (plots on the left) and the invariant mass of the Higgs and the leading jet (plots on the right) for Higgs 
produced via gluon-gluon fusion. The upper plots show the distributions to Leading Order (solid lines) and to Next-to-Leading 
Order (dashed lines). The lower plots show the ratio of the Next-to-Leading Order to the Leading Order cross-sections. The 
package MCFM was used to evaluate the cross-sections. Cross-sections are given in fb per bin. 
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Fig. 12.7.58: Same as Figure 12.7.57 after requiring a veto on an additional parton with Pt > 30 GeV and \-q\ < 2. 
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After the application of the cuts described above in this section the ratio of the NLO to the LO 
cross-sections remains large, in the range 1.5 1.6. 

A similar analysis was performed with the Higgs production via weak boson fusion. Before the 
application of a jet veto the impact of QCD NLO corrections on the Higgs Pt and Mhj distributions is 
mild and have little impact on the sensitivity of the channel. 

12.8 Summary 

We have proposed a method to estimate the normalization of the dominant backgrounds in the H — > 
W^W~ I'^l^uu channel using two control samples in the data, one b-tagged, and the other b- 
vetoed; in our approach, the systematic errors must be given in terms of the ratios aww^ cttt, and . 
We have computed the theoretical uncertainty on aww\ the result is 5%. We have shown that, for a 
b-tagging algorithm operating only on jets with Pt > 20 GeV and |r/| < 2.5, such that = 60% and the 
rejections against light quarks and c-quarks are 100 and 10 respectively, the effect of singly-resonant and 
non-resonant W~^W~bb diagrams is less than 10% on au and . A study using these uncertainties 
and this background extraction technique is in progress; the preliminary result is that a Higgs discovery 
at Mh = 160 GeV would require less than 2 fb^^ of integrated luminosity [279]. However, final 
calculations of the uncertainties on these last two extrapolation coefficients, as well as final results on 
the overall sensitivity of the search we have presented here, must be computed within the context of the 
LHC experiments. 

We have evaluated the QCD NLO corrections for signal in the Higgs search with H — > r+r" 
in association with one high Pt jet. The ratio of the NLO to LO cross-sections for Higgs production 
via gluon-gluon fusion is well above 2 for Higgs Pt > 50 GeV and increases with Mhj- The ratio 
drops substantially with the application of a veto on events with an extra parton with Pt > 30 GeV and 
\r]\ < 2. After the application of analysis cuts the ratio of the NLO to LO cross-sections for Higgs signal 
produced via gluon-gluon fusion is in the range 1.5 1.6. 
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13 An invisibly decaying Higgs at Tevatron and LHC 
Contributed by: H. Davoudiasl, T. Han, H.E. Logan 

The Higgs particle is the only missing part of the highly successful Standard Model (SM) of 
particle physics. The current experimental data from direct searches [131] and electroweak precision 
measurements [280, 281, 282] point to a Higgs mass in the range 114 GeV < rrih < 250 GeV. Thus, if 
the Higgs exists the Tevatron might detect it in the next several years and the LHC is expected to discover 
it. 

Most analyses assume that the Higgs will predominantly decay into detectable SM fields. How- 
ever, this may not be a good assumption if there are new weakly interacting particles with mass less 
than half the Higgs mass that couple to the Higgs with 0(1) strength. In this case, if rrih < 160 GeV 
~ 2 mw so that the Higgs partial width into SM particles is very small, the Higgs will decay predomi- 
nantly into the new weakly interacting particles. In particular, if these new weakly interacting particles 
are neutral and stable, the Higgs will decay invisibly. There ai^e many models in which this situation is 
realized, such as the Minimal Supersymmetric Standard Model (MSSM, with Higgs decays to lightest 
neutralinos), models with extra dimensions (with Higgs decays to Kaluza-Klein neutrinos [283]), and 
Majoron models [284, 285]. An invisible Higgs is also quite generic in minimal models of dark matter 
containing a stable singlet scalar- [286, 287, 288]. The combined LEP experimental bound on the mass 
of an invisibly-decaying Higgs boson is 114.4 GeV at 95% confidence level [289]. 

In this work, we study the discovery potential for the invisible Higgs hinv at the LHC and the 
Tevatron in the channels Z + hinv and hinv + 33 in Weak Boson Fusion (WBF). There have been a 
number of similar studies in the past [290, 291, 292, 293, 294, 295, 296, 297, 298]. We also examine the 
prospects for determining the mass of the invisible Higgs from production cross sections at the LHC. We 
show that the Z + hi^v channel gives a surprisingly good handle on the Higgs mass given 100 fb^^ of 
integrated luminosity. We also show how the Z + h\nv and WBF channels can be combined at the LHC 
to remove model assumptions from the Higgs mass extraction. A more detailed account of this study can 
be found in Ref. [299]. 

13.1 Production of /iinv via WBF at the Tevatron 

WBF production of the invisible Higgs was studied for the LHC in Ref. [295], which showed that WBF 
can provide significant signals for invisible Higgs discovery, even at low luminosity. Here, we will use 
their approach to show that WBF contributes significantly to the observation of hi^v at the Tevatron. 
Even though a 3(7 observation of a 120 GeV /linv in any single channel at the Tevatron is not possible 
with less than 12 fb^^ per detector, one can enhance the significance of the signal by combining data 
from various channels. At the Tevatron, an important production mode is Z + hi^v [294] and yields a 
somewhat larger significance than the WBF channel that we study. Combining these two channels and 
data from two Tevatron detectors, we show that a 3a observation of /linv with mh = 120 GeV can be 
obtained with 7 fb~^ of integrated luminosity per detector. 

At the LHC, the kinematic requirements for suppressing the backgrounds rely on the large energy 
and rapidity of the forward tagging jets characteristic of WBF at the LHC, together with the large rapidity 
coverage of the LHC detectors. Despite the more limited kinematic range and rapidity coverage at the 
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Tevatron, we show that the WBF production mode will indeed have a significant impact on the prospects 
for the observation of /ijnv at the Tevatron, before data from the LHC becomes available. 

The signal here is + 2j. A large background comes from vv) + 2j with the jets pro- 

duced via QCD. A smaller, but less reducible, background comes from Z{-^ vv) + 2j in which the Z is 
produced by WBF and the jets have kinematics similar to that of the signal. In addition, there are back- 
grounds from W{-^ Iv) + 2j, in which the lepton from the W decay is missed, and QCD backgrounds 
with fake Pt from missed jets in multi-jet events and jet energy mismeasurements in di-jet events. 

We generate the signal, h\^^ + 2j, the QCD and electroweak backgrounds with Z{-^ vv) + 2j, 
and the QCD background with W{-^ iu) + 2j for the Tevatron using Madgraph [269, 242]. We start 
with the following "minimal cuts": 

PTij)> 10 GeY, |r?(j)|<3.0, ARijj) > OA, Pt> 90 GeV. 

(13.1.75) 

The Pt > 90 GeV requirement provides a trigger. We take the calorimeter pseudo-rapidity coverage 
from, eg.,Ref. [300]. 

We impose "WBF cuts": we require that the two jets reconstruct to a large invariant mass, 

rrijj > 320, 340, 360, 400 GeV, (13.1.76) 
and are separated by a large rapidity gap, 

Ar]jj > 2.8. (13.1.77) 

These two cuts eliminate most of the QCD Z + 2j and W + 2j backgrounds, in which the jets tend to 
be softer and have a smaller rapidity gap, while preserving a significant fraction of the WBF signal. 

To reduce the W + 2j background further, we apply a lepton veto. We veto events that contain an 
isolated electron with [301] 

PtW > 8 GeV, \r]{i)\ < 3.0. (13.1.78) 

For simplicity, we apply the same veto to W decays to muons or taus. Loosening the veto requirements 
to pri^) > 10 GeV, \ri{i)\ < 2.0 increases the W + 2j background by about a factor of two. 

Background can also come from QCD multi-jet events with fake Pt due to mismeasurement of 
jets and jet activity escaping down the beampipe. We follow the techniques of a CDF study of Pt + 2j 
[302] to deal with this background. Please see Ref . [299] for more details. 

In Table 13.1.13 we show results for signal and background cross sections for the rrijj cuts given 
in Eq. (13.1.76). In Table 13.1.14 we show the resulting signal-to-background ratio and significance for 
10 fb-^ 

We find a signal significance of about l.6a with 10 fb^^ of luminosity at one Tevatron detector. 
This significance is not much less than that found in Ref. [294] for Z + /linv at the Tevatron, namely 
1.9cj with 10 fb~^ for nih = 120 GeV. Combining data from both Tevatron detectors, a 3a observation 
would require at least 12 fb~^ in the Z + /linv channel, or 18 fb~^ in the WBF channel. However, 
by combining these two channels, we find that a 3a observation of hi^v is possible with 7 fb^^ per 
detector, if the background can be determined to better than 10%. Thus, WBF provides an important 
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rrijj cut S(/ii„v + 2j) B(Z + 2i,QCD) B(Z + 2j,EW) B(H^ + 2j,QCD) 

320 GeV 4.1 fb 55fb 1.7 fb 7fb 

340 GeV 3.6 fb 43 fb 1.6 fb 5 fb 

360 GeV 3.2 fb 34 fb 1.4 fb 5 fb 

400 GeV 2.4 fb 21 fb 1.2 fb 2 fb 

Table 13.1.13: Signal and background cross sections for hinv + '2j at Tevatron Run 2, for = 120 GeV. The statistical 
uncertainty on B(Z + 2j,QCD) after cuts is roughly 10% due to our limited Monte Carlo sample. There is an additional 
background from QCD with fake Pt which is taken from Ref. [302] to be 5 fb; this represents a conservative overestimate of 
the fake Pt background. 



rrijj cut 


S (lOfb-i) 


S/B 


S/VB(10fb-i) 


320 GeV 


41 evts 


0.060 


1.6 


340 GeV 


36 evts 


0.066 


1.5 


360 GeV 


32 evts 


0.070 


1.5 


400 GeV 


24 evts 


0.082 


1.4 



Table 13.1.14: Number of signal events, signal-to-background ratio, and significance for hinv + 2j at Tevatron Run 2, for 
TUh = 120 GeV. We include the background from QCD with fake Pt of 5 fb [302] in S/B and S/\/B. 



second channel that brings an observation of /linv into the realm of possibihty at the Tevatron before the 
results of the LHC become available. Here, we note that there may be other production channels, such 
as gg — > hinvjj, that could contribute to the signal, even after the WBF cuts we have outlined. However, 
this could only enhance /i^^ production, making our results for the WBF channel a lower bound on the 
number of signal events. 

In Refs. [295, 303], it is claimed that vetoing additional soft jets in the central region improves 
the signal-to-background ratio by a factor of three at the LHC. If a similar background reduction could 
be achieved at the Tevatron, the prospects for /linv observation in the WBF channel would improve 
considerably: a 3a observation in the WBF channel alone would then be possible with 6 fb~^ per de- 
tector, with a signal-to-background ratio close to 1/5. Further discussion of background reduction is 
presented in Ref. [299]. We emphasize that we have not applied a central jet veto to obtain the results in 
Tables 13.1.13 and 13.1.14. 



13.2 Associated Z + hi^v Production at the LHC 

Discovery of the Higgs in the Z + h-mv channel was studied for the LHC in Refs. [293, 296]. This channel 
was also analyzed for the Tevatron in Ref. [294]. In Ref. [293], the Z+jet background at the LHC was 
found to diminish the significance of the signal considerably, and the electroweak backgrounds coming 
from WW and ZW final states were ignored. We update and refine the analysis of Ref. [293] by taking 
into account sources of background not included in that study and considering a wider acceptance range 
for the leptons. We show that, with the kinematic acceptance and the cuts we adopt, the prospects for the 
discovery of the invisible Higgs in Z + hmv at the LHC are brighter than presented in Ref. [293], even 
with the WW and ZW backgrounds included. Our results are consistent with those of Ref. [296]. 
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We consider the production process 



pp ^ i+r) + /linv ; i = e,fi, (13.2.79) 

at the LHC. We assume that the Higgs decays 100% of the time to invisible final states, and that the 
production cross section is the same as in the SM. Our results can be easily scaled for other invisible 
branching fractions or non-SM production cross sections. The signal rate is simply scaled by the pro- 
duction rate and invisible branching fraction: 

5 = 5o^5^, (13.2.80) 

where 5*0 is the signal rate from our studies, a/asM is the ratio of the nonstandard production cross 
section to that of the SM Higgs, and BRj„t, is the invisible branching fraction. Assuming that the SM is 
the only source of background, the luminosity required for a given signal significance then scales like 



where Cq is the luminosity required for a given significance found in our studies. 



(13.2.81) 



Signal for /ijnv 

As the signal is £~^i~PT, the most significant sources of background are 

i+£-)Z{^ i^u), W+{^ i+u)W-{^ ru), £+£-)W{^ iu), (13.2.82) 

(with the lepton from the W decay in ZW missed) and Z + jets final states with fake Pt [293, 294]. We 
simulate the signal and the first three backgrounds for the LHC using Madgraph [269, 242]. 

We start with the following "minimal cuts": 

Pt(^^) > 10 GeV, |r/(£^)| < 2.5, AR{£+e-) > 0.4, (13.2.83) 

where r] denotes pseudo-rapidity and AR is the separation between the two particles in the detector, 
AR = a/ (A?7)2 + (A(^)2; is the azimuthal angle. The electromagnetic calorimeter at both AT- 
LAS [304] and CMS [305] covers the range |ry| < 3; however, the electron trigger covers only |ry| < 2.5 
(2.6) at ATLAS (CMS). The pseudo-rapidity acceptance for dielectrons could be expanded by requir- 
ing only one electron within \r]\ < 2.5 and the other within \r]\ < 3. Meanwhile, the muon trigger 
covers \r]\ < 2.2 (2.1) at ATLAS (CMS), with muon identification and momentum measurement out to 
1 77 1 < 2.4. We require |r/(^^)| < 2.5 for both leptons, so that the larger acceptance for dielectron events 
compensates the smaller acceptance for dimuon events. 

Because we will cut on the invariant mass of the dilepton pair to keep only events in which the 
dileptons reconstruct to the Z mass, we imitate the effects of LHC detector resolution by smearing the 
electron momenta according to 

AE/E = , © 0.5%, (13.2.84) 

^E{ GeV) 
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Fig. 13.2.59: Missing pr distribution for Z(— » e'^e ) + hmv signal (solid lines, with mh = 120, 140 and 160 GeV top to 
bottom) and backgrounds from WW and ZZ (dotted lines) at the LHC, after applying the cuts in Eqs. (13.2.83), (13.2.85) and 
(13.2.86). 



with the two contributions added in quadrature. This smearing has a negUgible effect on our results. We 
have thus apphed the same smearing to the final state with muons. 

The WW background can be lai^gely eliminated by requiring that the invariant mass 7n£+^- 
is close to mz- 

\m^+^- - mz\ < 10 GeV. (13.2.85) 

Also, the and i~ from two different parent W bosons tend to be more back-to-back than the leptons 
in the signal. We therefore impose an azimuthal angle cut on the lepton pair, 

A(t)i+i- < 2.5 or 143°. (13.2.86) 

This cut also eliminates Drell-Yan backgrounds with fake Pt caused by mismeasurement of the lepton 
energies. 

Our third cut is on Pt- The number of £^i~PT signal events typically falls more slowly with Pt 
than those of the ZZ or WW backgrounds, as shown in Fig. 13.2.59. 
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The final state £~^i~)W{—>- iu), where the lepton from the W decay is missed, can be a 

potential background. However, the probability of missing the lepton from the W decay is small given 
the kinematic coverage at the LHC. To reduce this background, we veto events with a third isolated 
electron with 

PT > 10 GeV, \r]\ < 3.0. (13.2.87) 

For simplicity, we apply the same veto to W decays to muons or taus. This veto reduces the Z + W 
background to the level of 5-10 fb, so that it has Uttle effect on the significance of the signal. 

We also include the background from Z + jets with fake Pt- As shown in Ref. [293], events of 
the type Z+jets can constitute a significant background due to jet energy mismeasurements resulting in 
fake Pt, or when one or more jets are emitted outside the fiducial region of the detector and ai^e therefore 
missed. The majority of those events can be eliminated by applying a jet veto, but those in which the 
jet(s) ai^e soft and/or escape down the beampipe can fake Z + Pt events. A simulation of the latter 
requires simulating the detector effects, which is beyond the scope of our analysis. Instead, as explained 
in Ref. [299], we adopt the results for this background from Ref. [293]. 

At this point, we note that there are other potentially large sources of background that need to 
be addressed [294]. The background events from Z* t~^t~ — > £^l~ Pt are efficiently suppressed 
by our Z-mass cut on m^+i-, the Pt cut, and the cut on A(p£+£^ that requires that the leptons are not 
back-to-back. This can be seen from Table 2 in Ref. [296], where it is shown that, after cuts similar 
to those we use, the resulting background from a single Z is basically absent for the ZH production 
channel. The same conlcusion is reached for the W + jet background in the ZH channel, in Table 2 
of Ref. [296]. Hence, fake events from W{-^ £i/)+jet, where the jet is misidentified as a lepton of the 
appropriate charge and flavor, are also ignored in our analysis. 

Our results for the background and signal cross sections are tabulated in Table 13.2.15. The 
corresponding signal to background ratio, S/B, and significance, S/ \fB, are tabulated in Table 13.2.16. 
We see from Table 13.2.16 that a > 5cr discovery can be obtained for mh = 120 GeV with 10 fb^^ of 
integrated luminosity, even with our conservative estimate for the Z+jets background for Pt > 75 GeV. 
With 30 fb~^, discovery can be pushed out to rrih = 160 GeV. 















S{Z + /linv) 




Pt cut 


B{ZZ) 


BiWW) 


B(ZW) 


B(Z + j)* 


rrih = 120 


140 


160 GeV 


65 GeV 


48.0 fb 


10.6 fb 


10.2 fb 


22 fb 


14.8 fb 


10.8 fb 


7.9 fb 


75 GeV 


38.5 fb 


4.3 fb 


7.4 fb 


9fb 


12.8 fb 


9.4 fb 


7.0 fb 


85 GeV 


30.9 fb 


1.8 fb 


5.5 fb 




11.1 fb 


8.3 fb 


6.3 fb 


100 GeV 


22.1 fb 


0.6 fb 


3.6 fb 




8.7 fb 


6.8 fb 


5.3 fb 



Table 13.2.15: Background and signal cross sections for associated Z{~* ) + hinv production at the LHC, combining the 
ee and ii/i channels. 'Estimated from Ref. [293] (see text for details). 

The Z + /linv channel can thus be used at the LHC for rrih < 160 GeV to supplement the 
WBF channel [295], which has higher significance. WBF production of h-mv at the LHC was studied 
in Ref. [295], which concluded that with only 10 fb~^ of integrated luminosity, /iinv can be detected 
at the > 5cr level up to rrih ~ 480 GeV. They also showed that the invisible branching fraction of 
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nih = 120 GeV 




ruh = 140 GeV 


ruh = 160 GeV 


Pt cut 


S/B 


S/VB(10fb-i) SA 


/B(30fb-i) 


S/VB (30fb-i) 


S/VB (30fb"i) 


65 GeV 


0.22 (0.16) 


5.6 (4.9) 


9.8 (8.5) 


7.1 (6.2) 


5.2 (4.5) 


75 GeV 


0.25 (0.22) 


5.7 (5.3) 


9.9 (9.1) 


7.3 (6.7) 


5.4 (5.0) 


85 GeV 


0.29 


5.7 


9.8 


7.4 


5.6 


100 GeV 


0.33 


5.4 


9.3 


7.3 


5.7 



Table 13.2. 16: Signal significance for associated Z{—* £^£ ) + /iinv production at the LHC, combining the ee and fifj, channels. 
The numbers in the parentheses include the estimated Z+jets background discussed in the text. 



a 120 GeV Higgs can be constrained at the 95% confidence level to be less than 13% if no signal is 
seen in the WBF^ hi^v channel, again with 10 fb~^. However, we would like to emphasize that the 
Pt measurements in the process t^l~ Pt that we studied here are largely determined by pt{U), and 
the distribution will suffer much less from systematic uncertainties compai'ed to the WBF where Pt is 
determined mainly from the forward jets. 

Higgs boson mass 

The Z + /linv channel may also provide an interesting handle on the Higgs boson mass, as follows. The 
mass of an invisibly-decaying Higgs boson obviously cannot be reconstructed from the Higgs decay 
products. Unless the Higgs is also observed in a visible channel, our only chance of determining the 
Higgs mass comes from the nih dependence of the production process. Extracting rrih from the pro- 
duction cross section requires the assumption that the production couplings are the same as in the SM. 
(Non-observation of the Higgs in any visible final state implies that the invisible branching fraction is 
close to 100%.) 

The Higgs mass extraction from measurements of the production cross sections in Z + hi^v and 
WBF are shown in Tables 13.2.17 and 13.2.18, respectively. There are two sources of uncertainty in 
the signal: statistical and from background normalization. The statistical uncertainty is Aas/crs = 
y/S + B/S. We estimate the total background normalization uncertainty for Z + hi^v to be the same 
size as that of the dominant process involving Z vu: AB/B = AB{ZZ)/B{ZZ). We assume that 
this background can be measured via the corresponding channels in which Z — > £~^£~ and take the 
uncertainty to be the statistical uncertainty on the Z £+i~ rate: AB{ZZ)/B{ZZ) ~ 7.1% (2.2%), 
for an integrated luminosity of 10 (100) fb~^. In Tables 13.2.17 and 13.2.18 we quote the resulting 
uncertainty on the signal cross section, given by Aas/crs = (B/S) x AB/B. The total uncertainty 
[Aa s / cr s]tot, presented in Tables 13.2.17 and 13.2.18, is then the sum, in quadrature, of the statistical 
and background uncertainties, as well as the other uncertainties given in the table captions. We then have 
Arrih = {l/p)[Aas/crs]tot- 

The cross section for Z + /linv production falls quickly with increasing nih due to the s-channel 
propagator suppression. This is in contrast to the WBF production, which provides a > 5cr signal up to 
nih — 480 GeV with 10 fb~^ if the Higgs decays completely invisibly [295]. Thus, while the statistics 
are much better for the WBF measurement than for Z + /ijnv, the systematic uncertainties hurt WBF 
more because {da s / dnih) / s is much smaller for WBF than for Z + /linv- The Z + ^jnv cross section 
is therefore more sensitive to the Higgs mass than the WBF cross section. 
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ruh ( GeV) 


120 


140 


160 


p = {das/dmh)/as (1/ GeV) 


-0.013 


-0.015 


-0.017 


Statistical uncert. 


21% (6.6%) 


28% (8.8%) 


37% (12%) 


Background normalization uncert. 


33% (10%) 


45% (14%) 


60% (19%) 


Total uncert. 


40% (16%) 


53% (19%) 


71% (24%) 


Am/, ( GeV) 


30(12) 


35 (12) 


41 (14) 



Table 13.2.17: Higgs mass determination from Z + hi-n^ with 10 (100) fb~^, assuming Standard Model production cross section 
and 100% invisible decays. The signal and background cross sections were taken from Table 13.2.15 for Pt > 75 GeV. The 
total uncertainty includes a theoretical uncertainty on the signal cross section from QCD and PDF uncertainties of 7% [306] 
and an estimated lepton reconstruction efficiency uncertainty of 4% (2% per lepton) and luminosity normalization uncertainty 
of 5% [307]. 



ruh ( GeV) 


120 


130 


150 


200 


p = {das/dmh)/as (1/ GeV) 


-0.0026 


-0.0026 


-0.0028 


-0.0029 


Statistical uncert. 


5.3% (1.7%) 


5.4% (1.7%) 


5.7% (1.8%) 


6.4% (2.0%) 


Background norm, uncert. 


5.2% (2.1%) 


5.3% (2.1%) 


5.6% (2.2%) 


6.5% (2.6%) 


Total uncert. 


11% (8.6%) 


11% (8.6%) 


11% (8.6%) 


12% (8.8%) 


Am/, ( GeV) 


42 (32) 


42 (33) 


41(31) 


42 (30) 



Table 13.2.18: Higgs mass determination from WBF hinv with 10 (100) fb^^, assuming Standard Model production cross 
section and 100% invisible decays. The background and signal cross sections were taken from Tables II and III, respectively, 
of Ref. [295], and include a central jet veto. The total uncertainty includes a theoretical uncertainty from QCD and PDF 
uncertainties of 4% [30, 29], and an estimated uncertainty on the efficiency of the WBF jet tag and central jet veto of 5% and 
luminosity normalization uncertainty of 5% [307]. 
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ruh ( GeV) 


120 


140 


160 


r = cTs(^/i)/fTs(WBF) 


0.132 


0.102 


0.0807 


{dr/dmh)/r (1/GeV) 


-0.011 


-0.013 


-0.013 


Total uncert., Ar/r 


41% (16%) 


54% (20%) 


72% (25%) 


Am/, ( GeV) 


36 (14) 


43 (16) 


53 (18) 



Table 13.2.19: Higgs mass determination from the ratio method discussed in the text, with 10 (100) fb^^. The event rates for 
WBF were interpolated linearly for Higgs masses of 140 and 160 GeV, which were not given explicitly in Ref. [295]. Statistical 
uncertainties were obtained assuming SM signal rates. The total uncertainty includes theoretical uncertainties from QCD and 
PDF uncertainties of 7% for Z + /linv [306] and 4% for WBF [30, 29], and estimated uncertainties on the lepton reconstruction 
efficiency in Z + hinv of 4% (2% per lepton) and on the efficiency of the WBF jet tag and central jet veto of 5% [307]. The 
luminosity normalization uncertainty cancels out in the ratio of cross sections and is therefore not included. 

More importantly, however, taking the ratio of the Z + /linv and WBF cross sections allows for 
a more model-independent determination of the Higgs mass. This is due to the fact that the production 
couplings in Z + /linv (hZZ) and in WBF (contributions from hWW and hZZ) are related by custodial 
SU(2) symmetry in any model containing only Higgs doublets and/or singlets. The production couplings 
thus drop out of the ratio of rates in this wide class of models (which includes the IVISSM, multi-Higgs- 
doublet models, and models of singlet scalar dark matter), leaving dependence only on the Higgs mass. 
The resulting Higgs mass extraction is illustrated in Table 13.2.19. Assuming SM event rates for the 
statistical uncertainties, we find that the Higgs mass can be extracted with an uncertainty of 35-50 GeV 
(15-20 GeV) with 10 (100) fb~^ of integrated luminosity. The ratio method also allows a test of the SM 
cross section assumption by checking the consistency of the rrih determinations from the Z + /linv and 
WBF cross sections alone with the m/j value extracted from the ratio method. Furthermore, observation 
of the invisibly-decaying Higgs in WBF but not in Z + hi^v allows one to set a lower limit on rrih in this 
class of models. 

We note that the Pt distribution is also sensitive to mh'. larger nih results in a larger average 
Pt in Z + /linv events. At the LHC, the production cross section and Pt distribution may be the only 
experimental handles on the mass of a Higgs boson with no visible decays. 
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14 Studies of Spin Effects in Charged Higgs Boson Production with an Iterative Discriminant 
Analysis 

Contributed by: S. Hesselbach, S. Moretti, J. Rathsman, A. Sopczak 

We report on detailed Monte Carlo comparisons of selection variables to separate tbH^ signal 
events from the Standai^d Model tt background using an Iterative Discriminant Analysis (IDA) method. 
While kinematic differences exist between the two processes whenever mfj± / myy±, the exploration 
of the spin difference between the charged Higgs and the gauge boson becomes crucial in the 
particularly challenging case of near degeneracy of the charged Higgs boson mass with the mass. 
The TAUOLA package is used to decay the tau leptons emerging from the charged Higgs and boson 
decays taking the spin difference properly into account. We demonstrate that, even if the individual 
selection variables have limited discriminant power, the IDA method achieves a significant separation 
between the expected signal and background. For both Tevatron and LHC energies, the impact of the 
spin effects and mass on the separation of signal and background has been studied quantitatively. 
The effect of a hard transverse momentum cut to remove QCD background has been studied and it is 
found that the spin effects remain important. The separation is expressed in purity versus efficiency 
curves. The study is performed for charged Higgs boson masses between the mass and near the top 
mass. 

14.1 Introduction 

The importance of charged Higgs boson searches has in the recent years been emphasized [308, 309, 4, 
310] for LEP, a future International Linear Collider (ILC), the Tevatron and the Large Hadron Collider 
(LHC), as the detection of a charged Higgs boson would be a definite signal for the existence of New 
Physics beyond the Standai^d Model (SM). Chai^ged Higgs bosons naturally arise in non-minimal Higgs 
scenarios, such as Two-Higgs Doublet Models (2HDMs). A Supersymmetric version of the latter is the 
Minimal Supersymmetric Standard Model (MSSM). It is a Type II 2HDM with specific relations among 
neutral and charged Higgs boson masses and couplings, dictated by Supersymmetry (SUSY) [66]. 

The Tevatron collider at Fermilab is currently in its second stage of operation, so-called Run 2, 
with a center-of-mass (CM) energy of ^/s = 1.96 TeV. This machine will be the first one to directly 
probe charged Higgs boson masses in the mass range up to mfj± rrit. Starting from 2008, the LHC 
at CERN will be in a position to confirm or rule out the existence of such a particle over a very large 
portion of both the 2HDM and MSSM parameter space, m^± < 400 GeV, depending on tan j3, the ratio 
of the vacuum expectation values of the two Higgs doublets (see the reviews [31 1, 312, 313] and a recent 
study [314]). 

At present, alower bound on the charged Higgs boson massexists from LEP [315], ^ m^yi, 
independently of the charged Higgs boson decay Branching Ratios (BRs). This limit is valid within any 
Type II 2HDM whereas, in the low tan (3 region (below about 3), an indirect lower limit on mu± can 
be derived in the MSSM from the one on niA (the mass of the pseudoscalar Higgs state of the model): 
ml± ~ +m\> (130 GeV)^. 

If the charged Higgs boson mass mu± satisfies mu± < rrit — mi,, where rrit is the top quark mass 
and rrih the bottom quark mass, bosons could be produced in the decay of on-shell (i.e., Tt — > 0) top 
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(anti-)quarks t bH^, the latter being in turn produced in pairs via gg fusion and qq annihilation. This 
approximation is the one customarily used in event generators when < m^. Throughout this study 
we adopt the same notation as in Ref. [316]: charged Higgs production is denoted by qq, gg ^ tt ^ 
tbH^ if due to (anti-)top decays and by qq, gg tbH^ if further production diagrams ai^e included. In 
fact, owing to the large top decay width (Tt ~ 1.5 GeV) and due to the additional diagrams which do not 
proceed via direct tt production [317, 318, 319], charged Higgs bosons could also be produced at and 
beyond the kinematic top decay threshold. The importance of these effects in the so-called 'threshold' 
or 'transition' region {mfj± ~ nit) was emphasized in Les Houches proceedings [63, 35] as well as in 
Refs. [316, 320, 321, 322], so that the calculations of Refs. [317, 318] (based on the appropriate qq, gg 
tbH^ description) are now implemented in HERWIG [323, 324, 325, 326] and PYTHIA [327, 328]. A 
comparison between the two generators was carried out in Ref. [316]. For any realistic simulation of 
production with > the use of these implementations is important. In addition, in the mass 
region near the top quark mass, a matching of the calculations for the qq, gg — > tbH^ and gb — > tH^ 
processes might be required [328]. 

A chai^ged Higgs boson with m}j± < nit decays predominantly into a r lepton and a neutrino. 
For large values of tan /3 ( > 5) the corresponding BR is near 100%. For mfj± > nit, — > TU-j- IS 
overtaken by — > tb, but the latter is much harder to disentangle from background than the former. 
The associated top quark decays predominantly into a boson, or at times a second charged Higgs 
boson, and a b quark. The reaction 

qq,gg^tbH^ {t ^ bW) {H^ ^ T^Vr) (14.1.88) 

is then a promising channel to search for a charged Higgs boson at both the Tevatron (where the dominant 
production mode is qq) and the LHC (where gg is the leading subprocess). If the TVr decay 

channel is used to search for Higgs bosons, then a key ingredient in the signal selection process should 
be the exploitation of decay distributions that are sensitive to the spin nature of the particle yielding the r 
lepton {H^ in the signal or in the background), as advocated in Refs. [329, 330, 331, 332] (see also 
[333, 127]). The r spin information affects both the energy and the angular- distribution of the r decay 
products. 

In the search for a charged Higgs boson signal containing a r lepton, not only the magnitude of the 
production cross section is important, but also the efficiency of identifying the r lepton in the hadronic 
environment plays a crucial role. Since r leptons have a very short life-time {r^ 10"^ s), they decay 
within the detectors and can only be identified through their decay products. In about 35% of the cases 
they decay leptonically and about 65% of the times they do so hadronically. Both of these decay modes 
are usually addressed in charged Higgs boson searches by employing dedicated r lepton triggers. The 
identification of taus in hadronic pp collisions has recently been studied, e.g. Z — > t~^t~ events [334] 
and further details are given in [335]. 

It is the purpose of this note to outline the possible improvements that can be achieved at the Teva- 
tron and LHC in the search for charged Higgs bosons, with mass below the top mass and including the 
appropriate description of the spin effects in the rur decay. In order to quantify the spin effect an 

Iterative Discriminant Analysis (IDA) method has been applied, which is a powerful tool to separate sig- 
nal and background, even in cases such as the one presently under study when several selection variables 
with limited discriminant power ai"e present. 
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14.2 Tevatron energy 

We start by studying charged Higgs production qq, gg tbH^ with subsequent decays t bW, —>■ 
TVr at the FNAL Tevatron with y/s = 1.96 TeV. In the following we analyze hadronic decays of the 
boson and r lepton {W^ — > qq' , r — > hadrons+i''T-), which results in the signature 26+2j -|-Tjot+p™'^^ (2 
6 jets, 2 light jets, 1 r jet and missing transverse momentum). The most important irreducible background 
process is qq,gg — > tt with the subsequent decays t —>■ bW^ and t bW~, one boson decaying 
hadronically (W^ qq') and one leptonically (W^ tu^-), which results in the same final state 
particles as for the expected signal. 

14.2.1 Simulation and detector response 

The signal process qq, gg tbH^ is simulated with PYTHIA [327]. The subsequent decays t — > bW^ 
(or its charge conjugate), — > qq' and — > tu-j- are also carried out within PYTHIA, whereas the 
r leptons are decayed externally with the program TAUOLA [336, 337], which includes the complete 
spin structure of the r decay. The background process qq, gg tt is also simulated with PYTHIA with 
the built-in subroutines for tt production. The decays of the top quarks and bosons ai^e performed 
within PYTHIA and that of the r lepton within TAUOLA. 

The momenta of the final b and light quarks from the PYTHIA event record aie taken as the mo- 
menta of the con^esponding jet, whereas for the r jet the sum of all non-lep tonic final state particles as 
given by TAUOLA is used. The energy resolution of the detector and parton shower and hadronization 
effects are emulated through a Gaussian smearing (A(pt)/pt)^ = (0.80/^^)^ of the transverse mo- 
mentum pt for all jets in the final state, including the r jet [4]. As typical for fast simulation studies, 
no effects of underlying events, are simulated. Events are removed which contain jets with less than 20 
GeV transverse momentum^^, coiTcsponding to about \r]\ > 3. The transverse momentum of the lead- 
ing charged pion in the r jet is assumed to be measured in the tracker independently of the transverse 
momentum of the r jet. The identification and momentum measurement of the pion is important to fully 
exploit the r spin information. In order to take into account the tracker performance we apply Gaussian 
smearing on 1/pf with cr(l/pf )[TeV"^] = v^O.522 + 222/(pnGeV])2 sin6l^, ^^ere 9^, is the polar 
angle of the vr. The missing transverse momentum p^^^^ is constructed from the transverse momenta of 
all visible jets (including the visible r decay products) after taking the modelling of the detector into 
account. The generic detector description is a good approximation for both Tevatron experiments, CDF 
and DO. 

14.2.2 Expected rates 

For completeness we present a brief discussion of the expected cross section of the charged Higgs boson 
signature under investigation. The signal cross section has been calculated for tan /3 = 30 and mfj± = 
80, 100, 130 and 150 GeV with PYTHIA, version 6.325, using the implementation described in [328], in 
order to take the effects in the transition region into account. Furthermore, it has been shown in [316] that 
the signal cross section for tbH^ agrees with the one from the top-decay approximation tt tbH^ for 
charged Higgs boson masses up to about 160 GeV if the same factorization and renormalization scales 
are used. Thus, we have used everywhere in this study the factorization scale {nit + fnH± ) 1 4 and the 

"in order to be largely independent of the specific detector performance, no requirement on thie jet resolution is applied. 
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renormalization scale m^^i for both signal and background (i.e., those recommended in [328] as most 
appropriate for the tbH^ signal)^", since the primary purpose of our study is to single out variables that 
show a difference between our and data samples and that this can unambiguously be ascribed to 
the different nature of the two kinds of bosons (chiefly, their different mass and spin state). In addition, 
the running b quark mass entering in the Yukawa coupling of the signal has been evaluated at mfj± . This 
procedure eventually results in a dependence of our background calculations on tan /? and, especially, 
mjj± that is more marked than the one that would more naturally arise as only due to indirect effects 
through the top decay width. Hence, the cross sections have been rescaled with a common factor such 
that the total ti cross section is = 5.2 pb [338]. To be more specific, we have first calculated 

the total cross section fjP_!'°<i'PYTHiA ^ y^j^]^ (j^g built-in routine for ti production in PYTHIA for 
all mij± = 80, 100, 130 and 150 GeV and then calculated from this the respective rescaling factors 
c{mu±) = 5.2 pb/(T^j°'^'^^^^^ {mfj±) for each m^±. Then we have calculated the background cross 
section for 'mu± = 80 GeV into the final state with the signature 2b + 2j + rjot + pf"^^ by enforcing 
the respective decay channels in PYTHIA using the built-in routine for ti production and multipUed it 
with c(80 GeV). In the same manner we have calculated the signal cross sections with the PYTHIA 
routines for tbH^ production by enforcing the respective decay channels in PYTHIA and multiplying 
with the rescaling factors c{mH±) for = 80, 100, 130, 150 GeV. The resulting cross sections are 

given in Table 14.2.20 before (cr*'^) and after (a) applying the basic cuts p[^^^ > 20 GeV and the hard cut 
^miss > xoo GeV. For the four signal masses, the tbH^ and ti tbH^ cross section calculations agree 
numerically. 





qq, 99 tt 


qq, 99 - 


^tbH^ 


mH± (GeV) 


80 


80 


100 


130 


150 


cj*'^ (fb) 


350 


535 


415 


213 


85 


f7(fb)forj4''*' > 20 GeV 


125 


244 


202 


105 


32 


a (fb) for > (20, 100) GeV 


21 


30 


25 


18 


7 



Table 14.2.20: Tevatron cross sections of background qq,gg tt and signal qq,gg tbH for tan/3 = 30 and mj^± = 
80, 100, 130 and 150 GeV into the final state 26 + 2j + rjct + vT'"" before (cr"') and after (cr) the basic cuts (pt > 20 GeV for 
all jets) and the hard cut (pr'"" > 100 GeV). 



14.2.3 Event preselection and discussion of discriminant variables 

The expected cross sections of the 2b + 2j + rjct + p™'^^ signature are of the same order of magnitude 
for the signal and background reactions, as shown in Table 14.2.20. Thus, the same number of signal 
and background events is assumed for the analysis of different kinematic selection variables. For the 
signal 5 • 10^ events have been simulated with PYTHIA for each charged Higgs mass at the Tevatron 
energy of 1 .96 TeV using the built-in ti routine in the ti tbH^ approximation, while for the ti 
background also 5 • 10^ events have been simulated using the built-in ti routine. Then the basic cuts 
^^cts ^ 2Q GeV are applied. An additional hard cut on the missing transverse momentum p^^^^^ > 100 
GeV is used to suppress the QCD background, as for example demonstrated in Ref. [333]. After the 
additional anti-QCD cut about 28000 to 42000 signal events, depending on the simulated charged Higgs 

^"clearly, for a proper experimental study, factorization and renormalization scales for our background process qq, gg 
tt — > tbW^ ought to be chosen appropriately, i.e., unrelated to the charged Higgs boson mass. 
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bosons mass, and about 30000 ti background events remain. Other background reactions, for example 
W+jet production, ai^e expected to be negligible because they have either a much lower production cross 
section or ai^e strongly suppressed compai^ed to tt background, as quantified for example in Ref. [333]. 
In addition to the previous study (based on 5000 x BR(r — > hadrons) events each) [335], the present 
one applies an IDA method [339] to explore efficiencies and purities. As already mentioned, particular 
attention is devoted to the study of spin sensitive variables in the exploitation of polarization effects for 
the separation of signal and background events. 

Figures 14.2.60-14.2.67 show examples of the signal and background distributions of some of the 
kinematic variables used in the IDA method and the respective difference between signal and background 
distributions, namely: 

• the transverse momentum of the r jet, pi'"'' (Fig. 14.2.60), 

• the transverse momentum of the leading vr^ in the r jet, pf^ (Fig. 14.2.61) 

• the ratio pj'^/pp"* (Fig. 14.2.62), 

• the transverse momentum of the second (least energetic) b quark jet, (Fig. 14.2.63), 



• the transverse mass in the rjet + pf^^^^ system, mt = -J Ipl'^^pf^^^^ll — cos{A<f>)], where A(f) is the 



azimuthal angle between pp""* andp™'**'^ (Fig. 14.2.64)^\ 

• the invariant mass distribution of the two light quark jets and the second b quark jet, rrijjb^ 
(Fig. 14.2.65), 

• the spatial distance between the rjet and the second b quark jet, AR{t, 62) = \/ (Ac/))^ + (A?])^, 
where Acj) is the azimuthal angle between the r and 6 jet (Fig. 14.2.66) and 

• the sum of the (scalar) transverse momenta of all the quark jets, -f/jcts = Pt^ + Pt^ + Pt^ + Pt^ 
(Fig. 14.2.67). 

The distributions of signal and background events are normaUzed to the same number of 10^ events, in 
order to make small differences better visible. 

The signal and background distributions for the variables shown in Figs. 14.2.63-14.2.67 are as 
expected rather similar for m//± = m-\Y± and are hence mostly important to discriminate between signal 
and background in the IDA for 'mH± > m\Y± . Especially the transverse mass. Fig. 14.2.64, shows a large 
variation with the charged Higgs boson mass. However, the different spin of the charged Higgs boson and 
the W boson has a large effect on the rjet variables and pj (Figs. 14.2.60 and 14.2.61) resulting 
in significantly different distributions of signal and background even for mf^± = my/±. Moreover, the 
spin effects in the and distributions are con^elated which can be seen in Fig. 14.2.62 where 
the distributions of the ratio pj jp^"^ [329, 331, 332] show even larger differences. This highlights the 
importance of the additional variable pj (and hence pi /Pt"^), compared to a previous study [335]. 
The large separation power of this variable is indeed due to the different r polarizations in signal and 
background as can be inferred from the lower plots in Figs. 14.2.60-14.2.62. There the signal and 
background distributions for p^"^ , p^ and p^ jp^"^ are shown for reference samples where the r decay 
has been performed without the inclusion of spin effects with the built-in routines of PYTHIA and hence 
the differences between signal and background nearly vanish. 

^'strictly speaking this is not ttie transverse mass since there are two neutrinos in the decay chain of the charged Higgs boson 
we are considering, even so the characteristics of this mass are very similar to that of the true transverse mass. 
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14.2.4 Iterative discriminant analysis (IDA) 

The IDA method is a modified Fisher Discriminant Analysis [339] and is characterized by the use of a 
quadratic, instead of a linear, discriminant function and also involves iterations in order to enhance the 
separation between signal and background. 

In order to analyze our events with the IDA method, signal and background have been split in two 
samples of equal size. With the first set of samples the IDA training has been performed and then the 
second set of samples has been analyzed. We have used the following 20 variables in the IDA study: 
the transverse momenta pj , pf^^^, p^^, , pj}, pj-^, pi''; the transverse mass rrit', the invariant 
masses mjj, rrijji,^, nijjiy^, rrihb and s = mjjbbr; the spatial distances Ai?(r, bi), Ai?(r, b2), AR{T,ji), 
AR{T,j2); the total transverse momenta of all quai^k jets i/jets and of all jets i^aii = -f^jcts + In 
the analysis of real data, b-quai^k tagging probabilities and the reconstruction of t and W masses could 
be used to improve the jet pairing, and replace the allocation of least and most energetic 6-jet by a 
probabilistic analysis. 

The results of the IDA study are shown in Figs. 14.2.68 and 14.2.69 for the event samples with 
spin effect in the r decays for = 80, 100, 130, 150 GeV and for the reference samples without the 
spin effect for m^± = 80 GeV in order to illustrate the spin effect. In all plots of the IDA output variable 
the number of background events has been normalized to the number of signal events. Two IDA steps 
have been performed. Figure 14.2.68 shows the IDA output variable after the first step, where 90% of 
the signal is retained when a cut at zero is applied. The signal and background events after this cut are 
then passed to the second IDA step. Figure 14.2.69 shows the IDA output variable distributions after the 
second step. A cut on these distributions leads to the efficiency and purity (defined as ratio of the number 
of signal events divided by the sum of signal and background events) combinations as shown in the lower 
right plot in Fig. 14.2.69. These combinations define the working point (number of expected background 
events for a given signal efficiency) and the latter can be optimized to maximize the discovery potential. 
The difference between the dashed (no spin effects in r decay) and solid (with spin effects in r decay) 
lines for 'mfj± = 80 GeV in the lower right plot in Fig. 14.2.69 stresses again the importance of the spin 
effects to separate signal and background. 

In order to illustrate the effect of the hard cut on the missing transverse momentum (p™'^^ > 
100 GeV), which is imposed to suppress the QCD background, the final efficiency-purity plot of the 
IDA analysis is shown in Fig. 14.2.70 for = 80 GeV for two reference samples (red, long dashed: 
with spin effects in the r decay; red, dotted: without spin effects) without imposing the hard cut. The 
black lines (dashed and solid) are for the samples with the hard cut as also shown in the lower right plot 
in Fig. 14.2.69. As expected the achievable purity for a given efficiency decreases with the hard cut, 
therefore the spin effects become even more important to separate signal and background. In principle, 
by choosing the signal reduction rates in the previous IDA iterations, the signal and background rates in 
the final distributions can be varied appropriately. However, we have checked that a different number of 
IDA iterations and/or different efficiencies for the first IDA iteration have only a minor effect on the final 
result. 
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Fig. 14.2.60; pt distributions of the r jet for the tbH^ signal and the tt background for = 1.96 TeV (left) and the respective 
differences between signal and background (right). The lower plots show distributions without spin effects in the t decays. 
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Fig. 14.2.61 : pt distributions of the leading tt* from the r decay for the tbH^ signal and the tt background for a/s = 1.96 TeV 
(left) and the respective differences between signal and background (right). The lower plots show distributions without spin 
effects in the r decays. 
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Fig. 14.2.63: pt distributions of the second (least energetic) b quark jet for the tbH^ signal and the tt background for = 
1.96 TeV (left) and the respective differences between signal and background (right). 
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Fig. 14.2.64: Transverse mass mt = y^2pp°'p™''^'^[l — cos(A(^)] distributions of the Tjot +Pt""''' system (A(j)is the azimuthal 
angle between and pT^^^) for the tbH^ signal and the background for y/s — 1.96 TeV (left) and the respective 
differences between signal and background (right). 
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Fig. 14.2.65: Invariant mass distributions of the two light quark jets and the second (least energetic) b quark jet for the tbH^ 
signal and the ti background for ^/s = 1.96 TeV (left) and the respective differences between signal and background (right). 
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Fig. 14.2.69; Upper row, middle row and lower left figure: distributions of the IDA output variable in the second IDA step for 
90% efficiency in the first IDA step (corresponding to a cut at in Fig. 14.2.68) for the tbH^ signal (solid, red) and the tt 
background (dashed, blue). Lower right figure: efficiency as a function of the purity when not taking the spin effects in the 
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Fig. 14.2.70: Efficiency as a function of purity for mu± = 
80 GeV and ^/s = 1.96 TeV. The black lines are the re- 
sults after applying the hard cut pT"^^ > 100 GeV when not 
taking the spin effects in the r decay into account (dashed) 
and with spin effects in the t decay (solid), as also shown in 
Fig. 14.2.69. The red lines are the results without applying the 
hard cut on when not taking the spin effects in the r de- 
cay into account (dotted) and with spin effects in the r decay 
(long dashed). 

The simulation procedure and the emulation of the detector response are the same as those outlined in 
Sect. 2.1 for the Tevatron, as well as, for the preselection and IDA method, as described in Sects. 2.3 
and 2.4, respectively. Hence, only the expected LHC rates are discussed, followed by the description of 
changes in the distributions of the variables and the final IDA results. 

Unlike the case of the Tevatron, where only charged Higgs masses smaller than the top quark 
mass can be explored, and 2HDM/MSSM signatures practically rely on rur pairs only, at the LHC the 
phenomenology is more varied. Here, the search strategies depend strongly on the charged Higgs boson 
mass. If mf^± < mt — mi, (later referred to as a light Higgs boson), the charged Higgs boson can 
be produced in top (anti-)quark decay. The main source of top (anti-)quarks at the LHC is again tt 
pair production (cr^f = 850 pb at NLO) [340]. For the whole (tan P,mf^±) parameter space there is 
a competition between the bW^ and bH^ channels in top decay keeping the sum BR(t — > bW~^) + 
BR(t bH^) at almost unity. The top quark decay to bW^ is however the dominant mode for most 
of the pai^ameter space. Thus, the best way to search for a (light) chai^ged Higgs boson is by requiring 
that the top quark produced in the tbH^ process decays to a W^. While in the case of decays r's 
will be tagged via their hadronic decay producing low-multiplicity narrow jets in the detector, there are 
two different decays that can be explored. The leptonic signature bbH^W^ —>■ bbrvlv provides 
a clean selection of the signal via the identification of the lepton / = e, fi. In this case the charged 
Higgs transverse mass cannot be reconstructed because of the presence of two neutrinos with different 
origin. In this channel charged Higgs discovery will be determined by the observation of an excess of 
such events over SM expectations through a simple counting experiment. In the case of hadronic decays 
bbH^W^ — > bbriyjj the transverse mass can instead be reconstructed since all neutrinos are arising 
from the charged Higgs boson decay. This allows for an efficient separation of the signal and the main 
ti — > bbW^W^ — > bbrvji background (assuming m//± > m\y±). The absence of a lepton (e or /i) 
provides a less clean environment but the use of the transverse mass makes it possible to reach the same 
mass discovery region as in the previous case and also to extract the charged Higgs boson mass. Both 
these channels show that after an integrated luminosity of 30 fb~^ the discovery could be possible up to 
a mass of 150 GeV for all tan/3 values in both ATLAS and CMS [314, 341, 342]. 

If the charged Higgs is heavier than the top quark, the dominant decay channels are TV and 

tb depending on tan /3. They have both been studied by ATLAS and CMS [343, 344, 345, 346]. 
The charged Higgs bosons are produced in the pp —>■ tbH^ channel. For the — > tb decay, a chai^ged 
Higgs boson can be discovered up to high masses (mjj± ~ 400 GeV) in the case of very large tan f3 
values and this reach cannot be much improved because of the large multi-jet environment. For the 
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— > TV decay mode this reach is larger due to a cleaner signal despite a lower BR. In this case the 5a 
reach ranges from tan /3 = 20 for mu± = 200 GeV to tan /3 = 30 for = 400 GeV. 

For the LHC, signal and background events have been simulated in the same way as for the Teva- 
tron as described before, however, without implying any rescaling factor to match a measured ti cross 
section. Table 14.3.21 lists the resulting cross sections before (a**^) and after (a) applying the basic cuts 
'l?^^^ > 20 GeV and the hard cut p^^^^ > 100 GeV. The LHC rates allow for the discovery to be less 
challenging than at the Tevatron in the region m^± ~ m^y ± , yet the separation of signal events from 
background remains crucial for the measurement of the charged Higgs mass. 





qq.gg tt 


99 - 


^ tbH^ 


mH± (GeV) 


80 
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cr*"^ (pb) 


45.5 


72.6 


52.0 


24.5 


9.8 


a (pb) for pj''*' > 20 GeV 


17.3 


33.9 


25.7 


12.2 


3.8 


a (pb) for > (20, 100) GeV 


4.6 


6.0 


4.8 


2.9 


1.2 



Table 14.3.21: LHC cross sections of background qq,gg tt and signal qq. gg tbH^ for tan /3 — 30 and mjj± — 
80, 100, 130 and 150 GeV into the final state 26 + 2j + rjct + pT'^^ before (cr*) and after (a) the basic cuts (pt > 20 GeV for 
all jets) and the hard cut (pj"'"" > 100 GeV). 

The kinematic distributions are shown in Figs. 14.3.71 to 14.3.78 for ^/s = 14 TeV The choice of 
variables is identical to the one for the Tevatron and allows for a one-to-one comparison, the differences 
being due to a change in CM energy (and, to a somewhat lesser extent, due to the leading partonic 
mode of the production process^^). The main differences with respect to Figs. 14.2.60-14.2.67 are that 
the various transverse momenta and invariant masses have longer high energy tails. In particular, it 
should be noted that the effect of the spin differences between and events can be explored very 
effectively also at LHC energies, e.g. the ratio /p^^"' is shown in Fig. 14.3.73 which is very sensitive 
to the spin effects. These observations lead to the conclusion that the same method using spin differences 
can be used to separate signal from background at both the Tevatron and the LHC. 

The distributions of the IDA output variables for the study at ^/s = 14 TeV for two steps with 
90% efficiency in the first step are shown in Figs. 14.3.79 and 14.3.80. These distributions are qualita- 
tively similar to those for the Tevatron (Figs. 14.2.68 and 14.2.69) and the final achievable purity for a 
given efficiency is shown in Fig. 14.3.80. As for the Tevatron energy a good separation of signal and 
background events can be achieved with the spin sensitive variables and the IDA method even in case 
mH± ~ m^y±. For heavier masses the separation of signal and background events increases due to 
the kinematic differences of the event topology. 



^^As the latter is dominated by qq annihilation at the Tevatron and gg fusion at the LHC. 
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Fig. 14.3.71: pt distributions of the r jet for the tbH^ signal and the tt background for ^/s — 14 TeV (left) and the respective 
differences between signal and background (right). The lower plots show distributions without spin effects in the r decays. 
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Fig. 14.3.72: pt distributions of the leading tt^ from the r decay for the tbH^ signal and the background for ^/s — 14 TeV 
(left) and the respective differences between signal and background (right). The lower plots show distributions without spin 
effects in the r decays. 
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Fig. 14.3.73: Distributions of the ratio pi /Pt'°^ for the tbH^ signal and the background for = 14 TeV (left) and the 
respective differences between signal and background (right). The lower plots show distributions without spin effects in the r 
decays. 
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Fig. 14.3.74: pt distributions of the second (least energetic) b quark jet for the tbH^ signal and the tt background for = 
14 TeV (left) and the respective differences between signal and background (right). 



131 




- = 80 GeV 
. m^, = 100 GeV 
. = 130 GeV 
. = 150 GeV 



20 40 60 80 100 120 140 160 180 200 220 

m, (GeV) 




Sig - Bl<g 









= 80 GeV 




ni„. 


= 100 GeV 




m„. 


= 130 GeV 




"Ih- 


= 150 GeV 





20 40 60 80 100 120 140 160 180 200 220 

m, (GeV) 



Fig. 14.3.75: Transverse mass mt = i^J 2pt-'°'p""^^[l ~ cos{A(f>)] distributions of the rjct system (A(f> is the azimuthal 

angle between and pT^^^) for the tbH^ signal and the tt background for — 14 TeV (left) and the respective differences 
between signal and background (right). 
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Fig. 14.3.76: Invariant mass distributions of the two light quark jets and the second (least energetic) b quark jet for the tbH^ 
signal and the background for ^/s — 14 TeV (left) and the respective differences between signal and background (right). 
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Fig. 14.3.77: Spatial distance AR{r, 62) = \/ {Aff))^ + i^v)^ distributions (where is the azimuthal angle in rad between 
the r and 6 jet) for the tbH^ signal and the background for y/s = 14 TeV (left) and the respective differences between signal 
and background (right). 
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Fig. 14.3.78: Distributions of the total transverse momentum of all quark jets, -ffjcts ~ Pt 1 + pf + + Pt^ for the 
signal and the ti background for y's = 14 TeV (left) and the respective differences between signal and background (right). 
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Fig. 14.3.79; Distributions of the IDA output variable in the first IDA step for the tbH signal (solid, red) and the tt background 
(dashed, blue) for = 14 TeV. 
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I m^^, = 80 GeV, hot spin effects 
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Fig. 14.3.80; Upper row, middle row and lower left figure: distributions of the IDA output variable in the second IDA step for 
90% efficiency in the first IDA step (corresponding to a cut at in Fig. 14.3.79) for the tbH^ signal (solid, red) and the tt 
background (dashed, blue). Lower right figure: efficiency as a function of the purity when not taking the spin effects in the 
r decay into account for m^± — 80 GeV (dashed) and with spin effects in the r decay for m^± — 80, 100, 130, 150 GeV 
(solid, from left to right). Results are for the LHC. 
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14.4 Conclusions 

The discovery of charged Higgs bosons would be a clear sign of physics beyond the SM. In this case 
study we have investigated charged Higgs boson topologies produced at the current Tevatron and LHC 
energies and compared them against the irreducible SM background due to top-antitop production and 
decay. While sizable differences between signal and background are expected whenever mf^± / myy±, 
near the current mass limit of about mf^± ^ 80 GeV the kinematic spectra are very similar between 
SM decays and those involving charged Higgs bosons. In this case, spin information will significantly 
distinguish between signal and irreducible SM background. In fact, we have considered hadronic tv^- 
decays of charged Higgs bosons, wherein the r polarization induced by a decaying (pseudo)scalar object 
is significantly different from those emerging in the vector (W^) decays onsetting in the top-antitop 
case. For a realistic analysis which is not specific for a particular detector, a dedicated Monte Carlo 
event generation and a simplified multipurpose detector response approximation have been applied. The 
identification of a hadronic tau-lepton will be an experimental challenge in an environment with typically 
four jets being present. We have demonstrated how an IDA method can be an applied to separate signal 
and background when the differences between the signal and background distributions ai^e small. Our 
results show that the IDA method will be equally effective at both the Tevatron and LHC. While only 
the dominant in^educible tt background has been dealt with in detail, we have also specifically addressed 
the QCD background. A suitably hard missing transverse momentum cut has been applied to reject 
such jet activity and we have demonstrated that although the discriminative power is reduced by such 
a cut, the reduction is small compared to the gain from including the r polarization effects. Using the 
differences in r polarization between the signal and the dominant SM in^educible tt background is crucial 
for disentangling the former from the latter. 
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15 Energy scale for b jets in D0 
Contributed by: J. Cammin 

This section describes the determination of the energy response of b jets in the D0 calorimeter. 
Since this measurement is work in progress, and no final results are available yet, we discuss only the 
concept of the measurement. 

15.1 Introduction 

In D0, jet energies measured in the calorimeter [347], are corrected for offset, response, and shower- 
ing [348]. The response is the largest single correction factor and is measured from the energy balance in 
7+jets events using the "Missing Et Projection Fraction Method" (MPF) [348]. The response is mea- 
sured in bins of an energy estimator, E' = Et^ cosh(?7jct), which later is mapped to the raw measured 
jet energy, E-^^'^ in order to get an energy-dependent measurement of the jet response. 

The above mentioned corrections are derived for light jets and do not take into account peculiar ities 
of heavy flavor jets, such as different fragmentation and hadronization and the presence of semileptonic 
decays. The latter leads to an energy response considerably smaller than that of light jets because of the 
neutrino involved in the decay. It is therefore crucial to derive special energy connections for b jets, so 
that particle masses measured in decays containing b quarks, such as Z — > bb, t bW, or H — > bb, are 
reconstructed at the correct energy scale. 

In D0 special b jet corrections exist only for semi-muonic decays, which are applied if the jet 
contains a "soft" muon. The following section describes the response measurement to compensate for 
the remaining effects (mainly due to semi-electronic decays, for which no dedicated corrections exist). 

15.2 The concept of the measurement 

The response for b jets that do not have a soft muon tag is measured in b-tagged 7+jets events using the 
same method (MPF) as for the light jet response described above. However, b-tagged jets are a mixture of 
true b and c jets and mistagged light jets which have different energy responses. The situation is sketched 
in Figure 15.2.81, where the measured mean value of the response distribution in a particular energy bin 
is the weighted sum of the mean values of response distributions for light, b, and c jets. The response 
of b jets thus needs to be disentangled from the response measured in a tagged 7+jets sample. Since 
there are three unknows, the mean values Rj, Rb, and Rc of the responses, this can be accomplished by 
algebraically solving a system of three equations, where each equation corresponds to a measurement of 
the response in a sample with different flavor composition of the jets. In D0, the first measurement is 
performed on an untagged sample. With good approximation, the measured response is that of light jets, 
Rj, Equation (15.2.89). The second measurement is obtained from a b-tagged 7+jets sample using the 
"counting signed impact pai^ameter" (CsiP) algorithm^^ [349], Equation (15.2.90). The third equation is 
taken from a sample with a tighter b-tagging criterion: The jets must be tagged with the CsiP algorithm 
and have a "track mass" mtrk above 1.93 GeV, Equation (15.2.91). rratrk is the invariant mass of the 
tracks that tag the jet, i.e., tracks with impact parameter significance greater than two. 

^^With this algorithm a jet qualifies as a b jet if it has at least two tracks with impact parameter significance a above three, 
or at least three tracks with a above two. 
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Fig. 15.2.81: Sketch of response distribution in a particular 
energy bin. 




2 2.5 
AM:„„ GeV 



Fig. 15.2.82; "Invariant track mass" minus an offset of 



untagged : Rut ~ Ri (15.2.89) 
tagged: Rt = fiRi + fbRb + fcRc (15.2.90) 
tagged: R„^t = flRi + fl,Rb + f^Rc (15.2.91) 

The distribution of the invariant mass mtrk is depicted in Figure 15.2.82 for light, b, and c jets. 
The energy response of b jets is then obtained by solving the system of Equations (15.2.89)- (15.2.91) 
and measuring Rut, Rt, and Rmt in various energy bins: 



Rb 



fb fjb 



J. fc 



Rt - Rut ( fl — ^f'l ] — Rmtyj 



(15.2.92) 



f'\ ('\ (i\ 

The flavor fractions ^ fl ^ fc can be obtained from fits of the mass templates to the data 
distribution as shown in Figure 15.2.82, or from a similar distribution that discriminates between the jet 
flavors. The flavor composition is also a function of the energy and must be measured separately in each 
energy bin. 

The advantage of this method is that the energy response of b jets is measured directly in data 
and relies only very little on Monte Cai^lo simulations (template distributions for the fit to the flavor 
fractions). However, an inclusive 7+jets sample contains only a few percent of 7+b events, hence a 
large data sample is needed in order to keep the statistical uncertainties at a reasonable level. 

Since jets in D0 are already corrected for the light jet energy scale, the resulting b response will be 
provided as a residual scale factor Rj/Rb. This scale factor derived from tagged 7+jets data can also be 
compared to Monte Carlo simulations of 7+jets and 7+b events, a preliminary result of which is shown 
in Figure 15.2.83. The study suggests that b jets need additional energy corrections of as much as 10% 
at energies around 20 GeV and about 5% at energies of 150 GeV. 
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16 Insights into ^ 77 from CDF Searches 
Contributed by: S.-W. Lee 

We describe how all of the diphoton measurements at CDF provide important insights into the 
Higgs search at the LHC, H 77. A brief review of diphoton physics at CDF is also given here. 

16.1 Introduction 

The study of photon production at a hadron collider is important for many reasons. As the photon 
energy is well-measured, compared to jets, it can be a good tool to further our understanding of Quantum 
Chromodynamics (QCD). One of the important reason to study photons at a hadron colUder, as well as 
all QCD measurements, is that they are the backgrounds to new physics. The most famous of these is 
the Higgs search at the LHC, where diphoton backgrounds are the most serious experimental difficulty. 
In the next section we will illustrate how all of the diphoton measurements at CDF provide important 
insights into backgrounds for new physics, specifically the diphoton backgrounds to the Higgs search at 
the LHC, H 77. 

In addition, there are a large number of important and well-motivated theoretical models which 
make a strong case for looking for new physics in events with two photons in the final state. These 
theories include Supersymmetry (SUSY), Extra Dimensions (ED), Grand Unified Theories, Composite 
models of quarks and leptons, and Technicolor models. Therefore it is important to understand diphoton 
production at Tevatron experiments in order to reliably search for the Standard Model Higgs and new 
physics at LHC. 

The aim of this talk is to present the recent measurement of diphoton production at Tevatron 
experiment, CDF, to lead us to a deeper understanding of new physics signatures at LHC experiments. 

16.2 Diphoton physics at CDF 

A brief review of physics with diphoton final states using the CDF detector at the Tevatron is given 
here. These include searches for supersymmetry, extra dimensions and bosophilic Higgs, as well as 
QCD diphoton cross section measurement. Recent results from CDF Run II experiment are presented, 
but some result from Run I is also reviewed. 

Diphoton final states are a signature of many interesting processes. For example, at the LHC, 
one of the main discovery channels for the Higgs boson seaixh is the 77 final state. An excess of 77 
production at high invariant mass could be a signature of large extra dimensions, and in many theories 
involved physics beyond the standard model, cascade decays of heavy new particles generate a 77 sig- 
nature in the final state. However, the QCD production rate is large compared to most new physics, so an 
understanding of the QCD production mechanism is a prerequisite to searching reliably for new physics 
in this channel. 

CDF has good analysis tool to identify the photon signal from the mixture of photons and a neutral 
meson background. For the CDF measurement the fraction of photon candidate events that have an 
observed conversion in the materials just in front of the calorimeter is used, along with the transverse 
shower shape measured in a proportional chamber at shower maximum in the calorimeter itself. In 
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the end one of the two methods is used to evaluate point-by-point the fraction of photons in the data 
sample. [350]. 

Diphoton Cross Section 

Recently CDF has performed pure QCD test with prompt diphotons using a data sample of 207 pb^^ in 
Run II [351]. The analysis required two photon candidates with Et > 14 GeV (13 GeV) for the leading 
(softer) photon candidate in the event. The background from non-prompt photon sources is determined 
using a statistical method based on differences in the electromagnetic showers. 




Fig. 16.2.84: The 77 qr distribution, along with predictions from DIPHOX (solid), RESBOS (dashed), and PYTHIA (dotted). 
The PYTHIA predictions have been scaled by a factor of 2. Also shown, at larger qr, are the DIPHOX prediction (dot-dashed) 
and the CDF data (open squares) for the configuration where the two photons are required to have Aff) < ■k/2. 

CDF has measured the cross section for prompt diphoton production as a function of three kine- 
matic variables - diphoton mass, the transverse momentum of the diphoton system (qr), and the az- 
imuthal angle between the two photons, A0. Comparisons have been made with predictions from 
DIPHOX, RESBOS and PYTHIA. The data are in good agreement with the predictions for the mass 
distribution. At low to moderate qx and Ac/) greater than ir/2, where the effect of soft gluon emissions 
are important, the data agree better with RESBOS than DIPHOX. By contrast, in the regions where the 
2— >3 fragmentation contribution becomes important, large qx, A.(p less than 7r/2 and low diphoton mass, 
the data agree better with DIPHOX. The qT distribution is shown in Fig. 1. 

This result would appear to indicate a need to have a full theoretical calculation of diphoton pro- 
duction; a resummed full NLO calculation will be necessary. Again, an understanding of the QCD 
diphoton production mechanism is a prerequisite to searching reliably for new physics in this channel. 



Search for Supersymmetry 

Among various SUSY models, two SUSY breaking mechanism are interesting, which predict photons in 
the final states. Supergravity models can produce events which decay down to the second lightest neu- 
tralino via a loop into the lightest neutralino (xj) and a photon, where x? is the lightest supersymmetric 
particle (LSP). Gauge-Mediated SUSY breaking models (GMSB) with the x5 decaying into a photon 
and gravitino can produce a final state of two photons and large missing transverse energy (^r)- 
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Fig. 16.2.85: The NLO cross section and combined experimental limits as a function of chargino and neutralino mass in GMSB 
model. 

is often used as a pointer to possible SUSY signals because it indicates the escape of a non-interacting 
SUSY particle from the detector. The LSP signals are of particular interest as they provide a natural 
explanation for the dark matter. 

CDF has searches 202 pb~^ of inclusive diphoton events of Run 11 data for anomalous production 

of as evidence of new physics. Events are selected as having two photon candidates with Et > 13 

GeV in the central. CDF observe no candidate events, with an expected standard model background of 

0.27 lb 0.07 (stat) it 0.10 (syst) events. Using these results, CDF has set limit on the lightest chargino 

M~± > 167 GeV/c^, and the lightest neutrahno M~o > 93 GeV/c^ at 95% C.L. in a GMSB scenario 
XT '6 

with a light gravitino [352]. 

Fig. 2 shows the combined CDF and DO result for the observed cross section as a function of xf 
and Xi along with the theoretical LO and NLO production cross sections. The final mass limit for the 
lightest chargino is 209 GeV/c^ which translates to a mass limit of 1 14 GeV/c^ on the lightest neutralino 
and a limit of 84.6 TeV on A. This result significantly extends the individual experimental limits [353]. 

Search for Extra Dimensions 

Recent theories postulate the existence of new space-time dimensions. Such extra dimensions might 
be found by looking for graviton exchange processes in the diphoton final state. For example in the 
Randall-Sandrum model with a warped extra dimension, diphoton resonances can be produced via the 
graviton. 

CDF has searches for diphoton mass resonance with a data sample of 345 pb~^. Two isolated 
photons , each with Et > 25 GeV, ai^e required in the analysis. The main background comes from 
standai^d model diphoton production which accounts for 30% of the events, and from jets which fake 
photons. No deviation from standard model expectations is observed, and set upper limit on the cross 
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section times branching ratio of the Randall-Sandrum graviton production and decay to diphotons. The 
lower mass bounds obtained for the first excited states of the Randall-Sandrum graviton are 690 and 
220 GeV/c^ for coupling, k/Mpi = 0.1 and 0.01, respectively. [354] 
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Fig. 16.2.86: Combined 95% C.L. RS gaviton mass limit of the diphoton and dilepton searches 



Fig. 3 shows the combined 95% C.L. RS graviton mass limit of the diphoton and dilepton searches 
in the graviton mass versus k/Mpi plane. 

Search for Basophilic Higgs 

The signature of high mass photon pairs is attractive for seai^ches for new physics as the photon is the 
lightest gauge boson, and hence might be more easily produced in decays of new physics. There are 
models in which a Higgs boson could decay into two photons with a branching ratio much lai^ger than 
predicted in the standai^d model; bosophilic Higgs boson. 

In Run I CDF has searched for departures from standard model expectations for inclusive high 
mass diphoton production in association with a W or Z boson [355]. This analysis is complimentary 
to the diphoton cross section analysis, in which very strict photon selection requirements are used to 
reduce the large jet fake backgrounds maximizing signal significance, but which become progressively 
less efficient with Et for high energy photons. It is also complimentary to the recent diphoton+X search 
analysis which was focused on non-resonant diphoton signatures such as GMSB SUSY. 

CDF found no evidence for a resonant structure and set an upper limit on the cross section times 
branching ratio for pp — > — > 77 between 60 and 200 GeV/c^ (see Fig. 4). A 95% C.L. lower limit 
on the mass of a bosophilic Higgs boson (one which couples only to 7, W and Z with standai^d model 
couplings) is set at 82 GeV/c^. 

16.3 Conclusion 

In this article we summarize the current CDF experimental results of diphoton physics, test of standard 
model and searches for new physics in final states containing energetic photons, at Tevatron. we also 
describe how all of the diphoton measurements at CDF provide important insights into backgrounds for 
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Fig. 16.2.87; Upper limit at 95% C.L. on the 77 + W jZ cross section as a function of 77 mass. The dashed curve shows the 
prediction for cross section times branching fraction for a bosophilic ^ 77 with branching fraction from reference [1] and 
the cross section for associa ted Higgs production is a standard model NLO calculation from reference [356]. 

new physics, specifically the diphoton backgrounds to the Higgs search at the LHC. As we learned it 
is important to understand diphoton production in order to reliably search for the standard model Higgs 
and new physics. There may be an interesting connection between Tevatron and LHC. 
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